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On Projective Planes Of Order Nine 


By Marshall Hall, Jr., J. Dean Swift and Raymond Killgrove 


1. Introduction. There are four known projective planes of order nine: the de- 
sarguesian plane; the plane coordinatizable by a Veblen-Wedderburn system which 
is not a field; its dual; and a self-dual non-desarguesian plane. The last was orig- 
inally found by Veblen and Wedderburn [4]; its structure has been generalized by 
Hughes [3]. 

Each of these planes was known to have an elementary abelian addition in an 
appropriate ternary ring. A natural question was whether or not there were 
any other planes with this property. A search was made partly by hand, partly by 
machine for such planes. The result of the search is that no further planes exist 
with this property but that the three non-desarguesian planes all may be coordina- 
tized in more than one way with an elementary abelian addition. For notation and 
terminology used in this paper the reader is referred to [1]. 


2. The additive pencil and lines consistent with it. In a plane of order 9 with an 
elementary abelian addition we may assume initially three pencils of finite lines 
zr=c,y=c,andy = x+c,c = 0,1, --- , 8, where the addition in the last pencil 
is that of the elementary abelian group. Specifically this last pencil may be repre- 
sented by the latin square: 
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where, as is customary, in each row of the square: yoy:YoY/sysYsYey7ys 1S to be read as 
the second row of the permutation: 


4 = a Oe 
(2.2) (° 1 3 Ss a * 
Yo Yr Yo Ys Ya Ys Yo Yr Ys 


and represents a line whose nine finite points are (7, y;) 7 = 0,1, --- , 8. 

In our search it is sufficient to find the affine planes of order nine; the extension 
to projective planes by adding the points at infinity is trivial. 

It turns out that there are 2241 lines consistent with the pencil (2.1), these being 
represented by permutations (2.2) which do not have as many as two points (7, y;) 
in common with any one of the lines of (2.1). In practice the lines for which yy = 0 
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are first found. These are then subjected to the permutations given by (2.1) to 
obtain the full list. If Z is a consistent line and A, is a permutation of (2.1) then 
LA, is also and conversely. Indeed if L; and L2 are any two linesand P any permuta- 
tion, the number of intersections (coincidences) of ZL; and Lz is obviously equal to 
the number for L,P and L2P. In our case the group goes into itself under multiplica- 
tion by A, so that the consistency is maintained. 

The discovery of the 249 lines through (0, 0) may be done by hand or by simple 
mechanical sorting procedures. To check and to facilitate the consideration of other 
possible structures a program for the automatic computer SWAC has been written 
which will determine the lines consistent with any square which is the table of a loop. 

The problem of the search is now to find sets of 63 additional lines from the 2241 
such that each resulting total of 90 lines contains exactly one line joining any two 
distinct points, (a, b) and (c, d). For a = c or b = d we have the line x = a = 
or y = b = din one of the pre-determined pencils. 

The search described would be far too extensive in practice and a procedure 
must be found to reduce the total number of sets to be looked for and the total 
number of isomorphic planes to be found. To do this we consider automorphisms 
which preserve the 27 fixed lines of the three basic pencils. The group G which takes 
the three’ pencils into themselves consists of 81-48 automorphisms. Of these 81 
translations T have the form: 


(2.3) (x,y) > (x +a,y +b). 


The subgroup H fixing (0, 0) has 48 elements and is essentially the automorphism 
group of the elementary abelian group of order 9. These automorphisms are of the 
form: 


(2.4) (x,y) > (a(x), a(y)) - 
where a@ is a permutation of 0, 1, --- , 8 which leaves 0 fixed. This group contains 
a subgroup U of order 6. The permutations a of U are: 
a= I 
a2 = (345) (687) 
a; = (354)(678) 
wee ay = (36) (47) (58) 
as = (37)(48) (56) 
as = (38)(46) (57) 


The group H may now be given by 
(2.6) H = U + UVa; + Ucn; + Uary + Uerrs + Ua, + Uorg:r + Ue 


where 

a; = (12)(36)(48)(57) 

ay3 = (13472685) 

aig = (1428)(3765) 
(2.7) as, = (15862743) 

a3, = (16452387) 

a37 = (17832546) 

ou; = (1824) (3567). 
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To apply the automorphisms 7’ of 2.3 to a line L the product A, “LA, is com- 
puted where A, and A, are the permutations of (2.1) containing taking 0 to a and b 
respectively. To apply an automorphism a of H to L, a ‘La is calculated. As these 
automorphisms leave the lines (2.1) invariant they take consistent lines to con- 
sistent lines and this fact may be used in the construction and checking of the con- 
sistent lines. The discussion of the method by which the 2241 lines are obtained from 
249 may be simply phrased in terms of translations. 

There are still other automorphisms which exchange the pencils x = c (through 
(©) on L.), y = c (through (0) on L..) and y = x + c (through (1) on L,.) among 
themselves. The reflection R about y = z is given by 
(2.10) (x,y) = (y, 2). 


This is applied to L by computing L™. 
There is also the automorphism C: 


(2.11) (x, y) > (24,2 — y). 


C and R together generate a group K of order 6 which permutes (= ), (0) and (1) 
in all possible ways. 


3. Classification of consistent lines. The 249 lines which pass through (0, 0) 
and are consistent with (2.1) are interchanged by the group H. These lines fall into 
12 equivalence classes with respect to this group. 

We list here one representative of each class and the number of lines in each 
class. 


Number of 
Class Line Conjugates 

1. © 22.8.8) %..4: 6.4 1 

20216874 8 i S& 16 

3. 8% 224 7 #8. 84,806 8 

4. @:'23.48. 0B B.6)i4.28 8 

SS  #@.& 38 2 68156. 4:2 48 

(3.1) 6.28: &-8, 2%, 1,446 4 48 

fie! Gok, B23. 4, 4% 48 

8 03625814 7 6 

S . &@..4,.3. 7.4, &.8.4 6 

.. O13. @: & 1, 8. 4:21.72 48 

i). Ma & 4-4-2 4. 34 6 

i: 0.2: @..¢ ib, £7). 3 6 

Total 249 

The distribution of these lines may be represented as follows: 

Form of line Class 
1 pene BG 7ST = a hCUe 
(3.2) 02--- 1wwztesteeeeset ee 
Os --- (¢ = 3, ---, 8) . A 2 ote. ee oe oe be 


The classes tabulated in (3.1) and (3.2) are equivalences under H alone. With 
respect to G, obviously each of the 2241 lines is equivalent to one of the lines of 
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(3.1) but some of the classes may be further amalgamated. For example, let L, be 
the representative of class 2 in (3.1); if we apply to L» first the translation (x, y) 3 
(2 + 3, y + 8) and then the automorphism a; = (12)(36)(48)(57) of H, the repre 
sentative of class 3 results: 


As 'L2As = Agl2As = (063) (174) (285) - (12) (364857) (084) (165) (273) 
(12) (38) (47) (56) 


I 


(12) (36) (48) (57) (12) (38) (47) (56) (12) (36) (48) (57) = (12) (37) (46) (58) = Zs. 


Ina similar manner the classes 5, 6, 7 amalgamate under G. Specifically if we de- 
note the elements of H: (12)(34)(68) and (12) (274658), by ai and ag respectively 
and follow the rule of denoting the line (permutation) representing the ith class 
of (3.1) by L; we have: 


—1 
ay2Ae D;A iar mad Le 
and 
—1,-1 
ag Aj LsA garg = L, ° 


Further, classes 9 and 10 amalgamate under G but classes 1, 4, 8, 11, 12 remain 
distinct. However, using the group K which permutes the basic pencils, we find that 
classes 8 and 12 are equivalent under the transformation C while 11 and 12 are 
equivalent under R. Thus classes 8, 11, 12 amalgamate under the group K. These 
are the only new equivalences provided by K and there finally remain 6 classes: 1; 
(2, 3); 4; (5, 6, 7); (8, 11, 12); (9, 10). 

Although many of the procedures discussed in this section are essentially me- 
chanical in nature, programming is at least as difficult as the hand computations 
and the latter method was utilized for all these results. 


4. Construction of pencils through (0, 0). A plane 7 containing the additive 
pencil (2.1) must contain a line through (0, 0) and (1, 2), that is a line whose 
representation is of the form 02 --- . From (3.2) it is clear that such a line must 
belong to one of the classes 1-6. But, using the equivalences deduced under G, we 
may assume that it belongs to one of the classes 1, 3, 4, 6 since a line of class 2 is 
equivalent to one of class 3 and a line of class 5 will be transformed to one of class 6. 

Among these four classes we shall establish a preference order or hierarchy in 
the order 1, 6, 4, 3. That is, we shall eliminate all lines of class 1 from consideration 
when making a search using a class 6 line; again when using a class 4 line we shall 
eliminate any line of class 1, 5, 6, or 7; finally using a class 3 line for a start any 
line of classes, 1, 4, 5, 6, and 7 is excluded. This preference order helps to eliminate 
duplications since each plane will appear only with its most preferred line. 

The first stage of the search consists in finding, for each start, the possible pencils 
through (0,0). When Apo of (2.1) is eliminated and we ignore the first column, such 
a pencil is represented by a latin square of order 8 on the numbers 1, 2, --- , 8 whose 
first line is the designated start and whose subsequent lines are selected from the 
249 permissible lines through (0, 0) with due respect for preference order. 

A considerable portion of this search was carried out by hand. The complete 
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search was made independently by SWAC and the results compared for the dupli- 
cated parts. 

The mechanical procedure first finds for each of the four starts, the lines through 
(0,0), (1,7),7 = 3, --- , 8 consistent with the start; a further rejection is then made 
on the basis of preference. Next, a particular line through (0, 0), (1, 3) is selected 
and all remaining lines consistent with it are used in the high speed memory of 
SWAC to form by a continuation of these processes all possible latin squares. Then 
the next line through (1, 3) is taken and so on to completion. 

With the set of possible pencils before us, we may proceed in either of two ways 
(or a combination of them). We may eliminate equivalent pencils on the basis of the 
automorphism groups H and K which leave (0, 0) fixed, and attempt to complete 
the surviving representatives to planes or we may eliminate the pencils which can- 
not be completed and then discuss possible equivalences among the survivors. The 
first choice is better adapted to hand methods, the second to machines. 

In the cases which produced a large number of possible pencils it was desirable 
to introduce a secondary preference order. In particular, the class 1 line, L, , goes 
into itself under all substitutions of H; moreover lines through (0, 0) in classes 2, 
3, 4, 5 and 6 are not consistent with L, . Hence any line through (0, 0) consistent 
with L; may be carried by elements of H into one of L; , Lio , Ls , Lo , Lun , Ine . These, 
in this order, were adopted as secondary preferences and the total output of squares 
of type 1 was thus reduced to 50. 

Under the automorphisms of H-K these pencils fall into 17 classes. Similarly, 
the class 6 start produced 57 pencils In 11 classes; the class 4 start 12 pencils in 3 
classes. No pencils are found in class 3. 

In the mechanical procedure, the 249 possible lines through a point (2, 0) are 
first found by translating the 249 lines through (0, 0). Each of these lines is then 
compared with a square and those lines which intersect no line of the square in more 
than one point are stored and eventually tabulated with the square. It frequently 
occurred that in the tabulation for a particular square no line through (0, y) was 
present for at least one y. This square could be at once rejected as impossible to 
complete. 

In all routines used on SWAC the major operation was one of sorting. To facili- 
tate this sorting, positional coding was used. Each element y; of a line was assigned 
eight consecutive spaces in a register (since SWAC uses 36 bit words, double pre- 
cision was required to store a full line). The number 7 is represented by a 1 in the 
ith space and 0’s in the other 7; 0 is represented by 0’s in all eight places. Using 
complementation and extraction, the number of intersections was representable as 
the number of 1 bits in a resulting word pair. 


5. Representative pencils which are completable to a plane. Most of the squares 
found cannot be completed; further, most of the squares which do lead to planes 
do so uniquely. After a few additional pencils are added from the tabulations de- 
scribed above the remainder of the lines are essentially decided and it is a simple 
matter to write them down. 

When, in either order, the equivalences of pencils through (0,0) have been taken 
into account and completability determined, we find 12 pencils remaining, 8 of class 
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= 0, Class 4 representative pencils which can be completed 
posi- eeg31f¢86 5 @ 3 021768 5 43 
e9s6éfo FT 2 4 8 G Sseetztis 2 & 
. 23 6 ee 7 e9 ea ees7is8 ¢€ 
(2) 0 6 8 C24 SE 44 @868786848686 4 
4 9¢€sadgress fa Sse saeaesegazaai ez 
: Sere ges ze 4 9756518326 4 
5 e8eeGeitis2f 08 45167 38 2 
1 o212768664 82 
9 03658271 4 
6 O.426 4.82.8,23 8 & 
2 (445) 058273846 1 
e 6 34 t. 7.6.6 3 
4 0758 241 383 6 
7 Gus & &. 3.2.3. 7 
é Each of the 12 has a unique completion to a full plane except 1.44 which has nine 
6 completions. But, using our preference order, the class 6 plane may be removed 
2 from further consideration. In the completion of the 6.17 square the class 1 line, 
+ 354021687, occurs unequivocally. 
” The tabulations yield 8 possible pencils through (2, 0) consistent with 1.44. 
4 With one exception each of these is uniquely completable; the first has two com- 
. pletions. These pencils are listed here, omitting the lines x = 2, y = Oand y = x + 1. 
1 (2, 0) pencils for 1.44 
3 en eoeaF 65843 210876543 
5) 36074232861 360517428 
6 4802653217 470238651 
2 (1441) 570383281 6 4 (1.445) 580642137 
63058 147 2 63072658: 1 4 
4 750614238 740381265 
‘ 840137625 $5016437 2 
6 2108 76 5 4 8 7 ie eis FT S 
7 7s 9 et v2 O's ey ‘O'S 2S"? '@'e 
3 a S°O'O'3 2 3 T's see 'O 7 FF o's 
2 (1.4423) 5 70 1°@ 48 82 (146) F's O'R £ O'R TT 
6 3 'O 7 2 8's 2°46 o's @ 2 BT FS Be 
78 8.848 Fe’ ft 7's O'S & F' ee'a2 
S44 °'O" Ss 'e'O''! 8° 7 S's 8's bh 4 '8°F 3 
2108765 4 3 2108 76 5 4 8 
370285461 37028 5 46 1 
2 et ee ee oe ee ee 48063127 8 
(1443) 5 8 0642 13 7 (144.7) 5609724831 8 
6-4 ces 18-5 - 2c’ 6EeeRaA £2 Be 2 
730518 6 2 4 730651862 4 
6.8 653 6 4.80840 840362157 
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(2, 0) pencils for 1.44—continued 


2108 76 5 4 8 7s 2S TY OS SS 
»pewree sti Ty 4 sweUewe & 2 tt eS 

er @ezspees& i 6068 2%76ii1%5 

(1.44.4) 56072 43 1 8 (1448) 5 701 8 4 6 38 2 
65601478 383 2 Stereet &€ 8S 

40383 81 2 6 5 7502 48 3 61 

Pun eS 16483 7 Se2@Cegisazg?® 


To give a final separation of the completable planes it suffices to indicate the 
pencils through (3, 0) for 1.44.1. Again lines x = 3, y = 0, y = x + 6 are omitted: 


(3, 0) pencils for 1.44.1 
r4v7F¢68 56 2.6 32 18 503 6 2 7 4 
23604713 6 24706 3 1 5 8 
ono 4.8.2 16.8 7 35 40 2 16 8 7 
(1.44.11) 4 2 60 3 8 5 7 1 (1441.2) 4 1605 7 8 3 2 
es i OV Ss 4:2 83 oe 2°98 -6° 4 7 4 
7 oe 2°65 6 8 14 6S | 2 SF 8 2a 
So Se 6 ao & 2 Ss 2S or oS 2 Oy 


6. Equivalence of completed planes. It remains to recognize the distinct planes 
among the 19 completed planes (excluding that of class 6). For the purposes of 
tabulation we distinguish the four previously known planes as follows: 

A. The desarguesian plane—coordinatizable by GF(9); this appears as 1.35, 
this pencil through (0, 0) being readily recognizable as giving the multiplication 
table of the field. 

B. The Veblen-Wedderburn plane—coordinatizable by a near field. Again this 
multiplication may be recognized in the square given for 1.36. 

C. The dual Veblen-Wedderburn plane—This may be exhibited as the plane 
constructable from 1.44.2. 

D. The plane whose coordinatization in terms of an elementary abelian addi- 
tion was indicated by Hughes [3]. This appears in 1.44.1.! 

Of these, planes A, B, C have lines of the form y = xm + 6 in their standard form 
and, using the additive pencil (2.1) and the multiplicative pencil through (0, 0) 
tabulated in the last section, the remaining lines may easily be added by the reader 
who wishes a complete tabulation. For B the distributive law (2; + x2)m = xym + 
xem holds; the multiplications are automorphisms of the additive group. For C the 
distributive law x(m, + m2) = xm, + 2xmMz is valid. 

For D we give the complete tabulation of 1.44.1.1 below. In order that the reader 
may follow the details of the arguments on equivalence of planes we also list the 
completions of 1.44.1.2, 1.44.6, 1.4 and 4.3. 

We now maintain that the nineteen planes may be classified as follows: 

A. 1.35 

B. 1.36, 1.4, 1.44.6, 4.4 

C. 1.44.2, 1.44.1.2, 1.44.3, 1.44.4, 1.27, 1.28 

D. 1.44.1.1, 1.44.5, 1.44.7, 1.44.8, 1.26, 1.30, 4.3, 4.5. 

To begin with class B, 4.4 is the Veblen-Wedderburn plane coordinatized by 
the system in which —1 is not in the center. This system is designated by the letter 





du 
ta! 
an 


wom =] > 


So Or 
“Ie CO bo Ore Ww 


bo 


e the 
tted: 


on 


wd ns CO 


noe Ww Ow 


6 
6 1 
lanes 


es of 


1.35, 
ation 


1 this 
plane 
addi- 
form 
0, 0) 
pader 
im + 


S the 


pader 
t the 


d by 
letter 





ON PROJECTIVE PLANES OF ORDER NINE 241 


U in [1]; as shown in this reference (p. 274) the same plane results as that pro- 
duced by the nearfield. The translation L — A; ‘LA: or (2, y) > (x + 1, y + 2) 
takes 1.4 into 1.44.6. It is somewhat more difficult to show the equivalence of 1.44.6 
and 1.36. In the projective completion of 1.44.6 we make the following replacements. 


(0,0) — (#) 
(2) — (0, 0) 
(0) — (0) 

(1, 1) — (1, 1) 


(#)(¢ bd ot 22-828 3 yo 23 4 5 6 :) 
e~8 S23 6S" 6"T @ 0 123 4 5 6 8 
where the last replacement indicates that the line x = 1 is replaced by the line 
x = yin the manner shown. With the new coordinates 1.44.6 is transformed to 1.36. 
This establishes the equivalences under B. 

For C, the translation L — A; ‘LA; takes 1.27 to 1.44.3 and 1.28 into 1.44.4. 
The transformation L — L™ interchanges 1.44.2 and 1.44.1.2. The automorphism 
C + (a, y) > (a, x — y) interchanges 1.44.4 and 1.44.1.2. Finally the mapping 
L — L (12)(36) (48) (57) interchanges 1.44.3 and 1.44.4. This combination shows 
the equivalence of the planes in class C. 

In class D, L — A; ‘LA; takes 1.26 to 1.44.8 and 1.30 to 1.44.7. The mapping 
L — (06)(17)(28) £(012) (375648) takes 4.5 into 4.3. We introduce new coordi- 
nates into the projective completion of 4.3 by making the replacements: 


~I =] 


(2) > (@) 
(6) — (0) 
(0,0) — (0, 0) 
(8, 4) — (1, 1) 
(7)(9 123 45 6 7 8 a)(® 4268 3 1 
“Se 76186%22¢ 9° "© 426 8 3 1)° 
Under this change the plane 4.3 is found to be identical with plane 1.44.1.1. The 
transformation L — (354)(678) L(345)(687) takes 1.44.8 into 1.44.5. Also L - 


(345) (687) L(354) (678) takes 1.44.7 into 1.44.5. Finally, L — (13472685) L(1586- 
2743) takes 1.44.5 into 1.44.1.1. This completes the tables of equivalences. 


~I 
or or 


Plane 1.44.1.1 


eh or Oo vo te By Se Fe ge "921.68 7 83 4 
wee ae ee ee ee ee oe oe. Th Bl Bo Se 
ae ee Te Ok ML ek See ~ 2 BONES FG Sears 
’2°34.66 7844 '3 °* 2 6 4:90 3°4 82 86:97 
"4 87846617 '3° * €.66 4°00 2:2 @ Ss 
°>£4 4 868 7 23 9 1 ’ 6 ¢€@.3 2.8168 97 6 
? ' O.4 3 Oia ie 8-44 64° 7 3°32 0:2 % 
"78 e12320604 63 "7 £2.42 43 $1 .9-3 
"8672015 3 4 " £28 6 €& 2 2 i @ 


Note: Starred lines are common to 1.44.1.1 and 1.44.1.2. 
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Plane 1.44.1.2—continued 
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Plane 1.4—continued 
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Plane 4.3—continued 
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7. Concluding remarks. The failure to find, in this extensive search, any planes 
beyond those known for many years strengthens the conjecture that the four planes 
are a complete set for order 9. This number is the smallest order still in doubt as 
to the list of possible planes and also is the smallest order with more than a single 
plane [2]. It is still not possible within the limits of present mechanical equipment 
to make a complete search for order 9 as was done for 8. Yor one thing, no list 
of latin squares of order 8 yet exists. It is planned to continue the search with 
other likely additive pencils beginning with the cyclic group. 

The reader will have noted that, whereas the procedures used in the basic search 
were mechanized or easily mechanizable, the methods of section 6 were more 
involved and clearly ad hoc. No satisfactory mechanical way to identify two iso- 
morphic planes exists whether they be presented by a coordinate system or by an 
incidence matrix. The preparation of such a method is an interesting question. 
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Functional Approximations and Dynamic 
Programming 


By Richard Bellman and Stuart Dreyfus 


1. Introduction. In this note we wish to indicate some ways in which the theory 
of approximation can be used to increase the range of present day computers. 
Although we are primarily interested in applying these techniques to the functional 
equations occurring in the theory of dynamic programming, [1], it should be 
noted that these same methods are applicable, and even more readily, to the classi- 
cal functional equations of mathematical physics. 

What we wish to do is to trade additional computing time, which is expensive, 
for additional memory capacity, which does not exist. 


2. Dimensionality Difficulties. A typical problem arising in the theory of con- 
trol processes is that of maximizing a functional of the form 


T 
(1) J(y) -/ g(x, y) dt, 
where x and y are N-dimensional vectors related by the differential equation 
‘ dx 2 
(2) =" h(z, y), 2(0) =«. 


As we have discussed elsewhere at some length, [2], [3], questions of this nature, 
although nominally within the domain of the calculus of variations, in actuality 
cannot be reduced to the point of numerical solution by means of classical tech- 
niques. 
Writing 
(3) Max J(y) = f(e, T), 
v 


the theory of dynamic programming replaces the foregoing variational problem 
by that of solving the nonlinear partial differential equation 


(4) x = Max [g(c, v) + f(h(c, v), af/ae)], 
f(c, 0) = 0, 
where 
fo _ (Ff Ff ... of 
(5) %. (2,2, 2). 


For computational purposes, it is often convenient to use the approximate differ- 
ence equation 


(6) f(c, T + A) = Max [g(c, v)A + f(e + he, v)A, T)). 
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In terms of the capacities of modern computers, we have an extremely efficient 
algorithm if N = 1, a scalar problem, and a feasible algorithm if N = 2. If N = 3 
or more, we face fast memory difficulties if we attempt to proceed in a routine 
fashion. 

The reason for this is the following. To store a function of N variables in the 
usual way, we tabulate the values of the function at a set of lattice points within 
the domain of interest. If there are M different possible values of c: , of c2, and 
so on, the total number of grid points will be MM”. For M = 100, and N = 3, this 
yields a quantity outside of present capabilities. 

In various papers, [4], [5], we have indicated some methods which enable us 
to circumvent these difficulties. These methods combine analytic devices with 
the method of successive approximation. In this paper, we wish to present a new 
method, based upon approximation techniques, which appears to have wide ap- 
plicability. 


3. One-dimensional Case. In order to illustrate the application of the method 
in its simplest form, let us consider the problem of determining the sequence of 
functions {f,(c)}, n = 1, 2, --- , given by the recurrence relation 


file) 


Max g(c, v), 


(1) 
fnsi(e) = Max [g(c, v) + fa(hi(e, v))]. 


Let us suppose that c takes values only over [—1, 1] and that the function h(c, v) 
similarly assumes values over this interval for all c-values and all permissible 
v-values. 

The standard approach involves a grid of values in the interval [—1, 1] where 
the number of grid-points depends upon the accuracy that we desire. Let us pro- 
ceed in a different manner. In place of considering that the function is determined 
by the set of grid-points, we shall consider the function to be determined by a 
Fourier expansion in terms of a suitable orthonormal set. For the interval [—1, 1], 
a convenient set is the set of normalized Legendre polynomials. Thus, for some 
fixed value R, we write 


R 
(2) falc) = D> aden Pr(c). 
k=0 


This is, of course, an approximation comparable with that of using a finite number 
of grid-points. 

We could envisage using a power series expansion rather than an expansion 
in terms of orthogonal functions. Generally speaking, an orthogonal expansion 
is to be preferred, both on the grounds of accuracy and the grounds of ease of de- 
termination of the coefficients. As far as the calculation of the values of P;(c) are 
concerned, the simple three-term recurrence relations connecting the successive 
members of the sequence of Legendre polynomials make the computation of these 
values not much more difficult than that of the calculation of the powers c’. 

The function f,(c) is now replaced, as far as storage in the computer is con- 
cerned, by the set of coefficients 


(3) A a= [do,n Oy ee Ar n\. 
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Whenever we wish a value of f,(c), we compute it by means of the approximate 
relation in (2). Consequently, in computing the values of f,4:(c) by way of (1), 
we require only the storage of the set of values given by A, . 

There is, however, a difficulty in this approach. The coefficients are determined 
by the relations 


(4) ain = | falc)Pa(e) de. 


How are we going to calculate them? If we replace the integral by a simple Rie- 
mann sum, we are led back to the necessity for tabulating the values of f,,(c) at a 
set of grid-points. 

It is here that we invoke the theory of mechanical quadrature. In place of 
evaluating the integrals in (4) by means of Riemann sums, we use an interpola- 
tion formula of the following form: 


1 Ss 
(5) [, t(OPale) de = & asf.lo) Palos, 


j=l 
where the a; and c; are carefully chosen. 

If the quantities c; are chosen to be the zeros of the Legendre polynomial of 
degree S, and the coefficients a; are chosen to be the Christoffel numbers, the 
formula in (5) is exact if the integrand f,,(c)P:(c) is a polynomial of degree 2S — 1. 

Since the quantities a; and c; are tabulated, and in any case, could easily be 
determined in advance, the amount of effort entailed depends upon the choice of 
S. Since our original approximation, contained in (2), is equivalent to that of 
assuming that f,(c) is a polynomial of degree R, it would be reasonable to take S 
equal to R. However, there is no necessity for doing this. 

If we use the equation in (5) to compute the values of the coefficient set A, , 
it is clear that at each stage we require only the values f,(c;), 7 = 1, 2,---, R. 

The computation then proceeds in the following fashion. Given the values of 
fn(c;), we compute the coefficient set A, . Using these coefficients we can determine 
the values of f,(h(c, v)) which occur in (1) in the course of computing the new 
values fn4i(¢;). 


4. Discussion. Apart from the fixed set of instructions, and the values such as 
a;P,(c;) which are determined at the beginning of the process, we require the 
set An = [on , Gin, *** , Gr,n| to be retained in the fast memory at the nth stage. 
This is a set of (R + 1) values. 

On the other hand, a grid size of 6 over [—1, 1] would require 1/6 values. In 
one dimension, the difference between R + 1 and 1/4 is not particularly important, 
and the great amount of additional computation required by the method described 
above can more than outweigh this advantage. 

Consider, however, the two-dimensional case. The straight-forward approach 
based upon a grid size of 6 in the c; and c intervals requires (1/6 )* values of f, (1 , C2). 
On the other hand, if we set 


R 
(1) fn(Cy , C2) _ ) AijnPi(a)P (ce), 
i,j=0 


and proceed as above, we require only R(R + 1)/2 values, the coefficients 
Qijn,t,J ead 0, 1, “i , R. 
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Proceeding to three dimensions, we compare (1/5)* and R(R + 1)(R + 2)/6, 
the number of coefficients a;;x.,,. Let us use some typical values of 6 and R and 
compare the values, say 6 = .01 and R = 5, 10. 





Dimension (1/8) R=5 | R=10 
l 100 6 il 
2 } ae 4 55 
3 10° 56 | 286 
4 108 126 | 1001 
5 10° 252 =| 3003 





We see that variational problems involving four and five state variables which 
are completely untouchable by direct methods are within the scope of the method 
we have outlined. Combining this method with the Lagrange multiplier technique, 
and the method of successive approximations, we have a way of attacking pre- 
viously impregnable problems. 


5. A Numerical Example. Consider the problem of determining the sequence 
{v,} so as to minimize the function 


(1) Pr = Duk +a dos, 
k=1 k=1 
where 
(@) Ung = 2Un — Un +n, Uo = ©, - 
(b) v, must be chosen so as to keep un+: within the interval [—1, 1). 
Introducing for N = 1 and —1 < ¢ S 1, the function 


(3) fw(c) = Min Fy, 
{v} 


we derive the recurrence relation 


fx(c) = Min [c’ + Wo” + fya(2c — & + v))I, N 


IV 
bo 


(4) 
fic) =c. 


This yields a simple computational determination of the sequence {fy(c)}. 

As a test of the foregoing method, this sequence was first determined by the 
usual method based upon a grid of values over [—1, 1], and then following the 
procedure described above, using the approximation 


(5) f(c) = 2 ax. Pi(c), 
k=0 


and a value of S = 14 in (3.5). 
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Let {fi(c)} denote the sequence determined using the grid, and {f,*(c)} the 
sequence obtained via Legendre polynomials. A comparison of values is given 
below for k = 6. 





c | wO | Ko 
1.0 | 1.782 | 1.77 
8 | 1.370 | 1.36 
2 | 1538 | «145 
0.0 | .006 | 0.0 
—.2 | .202 | .20 
—.8 | 4.876 | 4.89 
—1.0 | 8.666 | 8.67 


As can be seen, the agreement is quite good. 


The RAND Corporation, 
Santa Monica, California 
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The Influence of Critical Boundary Conditions 
on Finite Difference Solutions of Two 
Point Boundary Value Problems 


By A. R. Mitchell 


1. Introduction. It is often assumed in the solution of two point boundary 
problems by finite difference methods that an accurate answer is obtained pro- 
vided the mesh length is small and the numerical calculations are carried out with 
sufficient accuracy. The object of the present note is to illustrate by means of 
worked out examples the extent of the disagreement which often exists, even 
for small mesh lengths, between corresponding exact solutions of ordinary differen- 
tial and difference equations. Although the differential equation chosen for purposes 
of illustration in the present paper is a particularly simple one, the difficulties out- 
lined are quite general and arise in the numerical solution of more complicated 
equations subject to two point boundary conditions. 

The equation considered is 


(1) TY + p/w + hy = 0, (k > 0) 
dx? dx 
with the associated second order central difference replacement 


(2) (1 + $ph)yrus — (2 — h’)y + (1 — Sph)ya = 0, 


where p, k are constants, y;41, Yr, Yr are the values of y at x = (r + 1)H, rH, 
(r — 1)H respectively, and h = +/kH (where H is the mesh length). In all cases 
the range of the problem is from x = 0 to x = L where 


(3) L= (N+ 1)H, 
and N is the number of internal nodes. Corresponding exact solutions of (1) and 
(2) are compared for two types of boundary condition. 

2. Function given at both ends of the range. Consider the boundary conditions 


y = Oatx = Oand y = Y at x = L. The solution of (1) becomes 


(4a) ink sin $4 — pr/ka (nil-Vie) y 
sin 31/4 — pl 





if p’ < 4, and 
(4b) seals = 4/p — 4/ kx gin-Vie) y 
sinh 44/p? — 41 
if p> > 4, where 1 = +/K L. No finite solution exists for general x when 








- 2rr mae ; 
(5) wer on (r = = my ***2 
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The solution of (2) becomes 





; 2 +. a si 
(6a) Yr = ( *) ge me y 
2 — ph sin (N + 1)0 
where 
2-—h 
cos 6 = Ea 
V4 — ph 
ifp +h < 4, and 
; N—yp . ’ r 
(6b) Yr pe 7 } (ry = 1,2,--- N) 
where 
—— 2-h -~i & inte © v— 4 
= 9 fe ph 
if p +h > 4. No finite solution exists at general nodes for ranges satisfying 
[x - i . Kr ’ - 
7) = ~/2(N 1)} 1 20S ———— }.. = 1,2,---N 
(7) L= SY2(N + | -4/1- ws 5 COs | (K : V) 


Ranges for which either the differential or the difference equation has no finite 
solution are said to be critical. 

It can be shown that for a prescribed value of K, a critical range of the difference 
equation given by (7) either approaches a critical range of the differential equation 
given by (5) or tends to infinity as N tends to infinity. Due to the presence of 
the mesh length, however small, in a difference calculation, there is never exact 
correspondence between critical ranges of differential and difference equations, 
and an attempt is now made to determine the effect of this lack of correspondence 
in critical range on the agreement between exact solutions of differential and as- 
sociated difference equations. 

Calculations are carried out for p = 0, ~/2, 1.9, 2.0, +/4.99 and for h = 0.1, 
0.2. Values of y are obtained from (4) at « = L/2, and from (6) at the 
node r = (N + 1)/2. The difference between corresponding exact solutions of 
the differential and difference equations is expressed as a percentage of the solu- 
tion of the differential equation, and the resulting values are quoted in Table | 
It should be pointed out that when x = L/2, (5) requires modification, and in 
fact no finite solution of the differential equation exists for 

t. (ry = 1,3,5, ---) 
Vi-p 
Also the sign in front of the percentage differences is plus or minus depending on 
whether the difference equation solution is greater or less than the differential 
equation solution. 

There are the following points to be noted in the table: 

(1) The values of | for which the differential equation has no finite solution 
are indicated, and it is clear that poor agreement exists between the exact solu- 
tions of the differential and difference equations in the vicinity of these critical 
values. 
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(2) As p approaches the value 2 from below, agreement between the exact 
solutions deteriorates and when p = 1.9, for example, the error is in excess of 2% 
when h = 0.1 for any 1 > 8.2. 

(3) Once p has attained the value 2, the critical ranges of the differential 
equation are at infinity and the percentage errors increase monotonically with 1. 

(4) The P-condition numbers are given for p = 1.9, 2.0, +/4.99 when h = 0.1. 
The magnitude of the P number indicates the order of difficulty of obtaining an 
accurate numerical solution of the appropriate difference equation. The numbers 
are large in the vicinity of critical ranges. 


3. Function and derivative given, one at either end of the range. Next con- 
sider the boundary conditions y = 0 at x = 0 and y’ + ay = B at x = L. The 
solution of (1) becomes at x = L, 


; 28 
(8a) yVk = — $$ 
5— pt V4 — pcot V4 — p hl 
if p < 4, and 
= 28 

(8b) q k Se OE 

Vv 6 — p+ Vp — 4 coth Vp? — 441 
if p > 4, where 1 = +/k L and 6 = 2a/+/k. The solution ceases to be finite when 


(9a) l= 2(rx — 9) 


where 


if p> < 4, and when 





. — log. ~—t VP ah. 
VP = 4 t=94+- VF +4 
if p’ > 4. It should be noted that the derivative condition at x = L causes critical 


ranges to exist for p’ > 4 as well as for p’ < 4. 
If the derivative condition at « = L is replaced by 


(9b) l= 


Qyw — [(2 — ht) + (2 + ph)B0h] yuss + (2 + phyh © = 0, 


Vk 


the solution of (2) becomes at the node r = N + 1, 





h(2 h 
me ee Bh(2 + ph) 


2 Sets o({2+ #) __ sin NO © 
IE RE el Path (; — ph) sin (N + 1)6 


(10a) 
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where 


+ (ap = 2008 + 30 


cos 20 = — =a 








l— pr 
if p’ + h’ < 4, and 
a Bh(2 + ph) 
(10b) y — P . ‘ Tw -- ue 
where 


ip ato thvp +h Rs 
” 2+ ph , 
if p +h’ > 4. No finite solution exists when 
‘ PP Kr —of} ‘ : 
lla) 1=2(N+1)]1— 4/1 ——=— —~ ws oe ae 
(lla) (N+ | / iW +17 | (K = 1,2,---N) 
where 
[(4 — hk’) — p} 
te = —___—_ aa ; 
me = 81 = Fhe) + Ip? — 2) 
if p +h’ < 4, and when 


2 oe aT ee N+1 
am [Ranga 2] 





h + hrv/p? +h? — 4-2 


_ (1 — yph’)h + 3nlh’ — 2)h—- Vp? +h? —4 
6(1 — ap*h*)h + gph? — 2)h + ~/p? + h? — 4 


if p> + h? > 4. Again, it can be shown that the critical ranges of the difference 
equation given by (11) approach the corresponding critical ranges of the differ- 
ential equation given by (9) or tend to infinity as N tends to infinity. 

Calculations are carried out for p = 1.9, 2, 3 with h = 0.1. A range of values 
of 6 is covered for each value of p. The values of y+/k/8 for the differential and 
difference equations are obtained from (8a) and (10a) respectively when p = 1.9, 
and from (8b) and (10b) respectively when p = 2, 3. In Table 2, the percentage 
differences between corresponding exact solutions of the differential and difference 
equations are quoted for p = 2, 3, whereas the actual differences between cor- 
respuading exact solutions are quoted in units of .0001 for p = 1.9. This is because 
the exact solution of the differential equation for p = 1.9 is sometimes so small 
that the percentage differences are misleadingly high. 

There are several points arising from Table 2. 

(1) An entry © in the table means that the differential equation has no finite 
solution for the appropriate values of / and 6. As in Table 1, poor agreement exists 
between exact solutions in the vicinity of these critical values. 

(2) When p = 1,9, critical values occur periodically for each value of 6, but 
when p = 2, 3 there is at most one critical range for each value of 6. For example 
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when p = 2, as 6 increases from 0 to 2, the value of the critical range increases 
from 1 to «. For any value of 6 > 2, there is no critical range. 

(3) The lack of agreement between exact solutions in the vicinity of critical 
ranges seems to be less pronounced than was the case in Table 1 


4. The effect of reducing the mesh length on the accuracy of finite difference 
solutions. So far, exact solutions of the difference equation have been obtained for 
mesh lengths of 0.1 or 0.2. Calculations are now carried out to show the effect of 
reducing the mesh length on the accuracy of finite difference solutions. 

Considered the differential equation 


- + pk — ky =0 (k > 0) 


subject to the boundary conditions y = 0 at x = 0 andy’ + ay = Batz = L. 
The values of y+/k/8 are calculated from the conventional second order central 


difference replacement for a range of values of § and N when »/kL = 4 and p = —6, 
and are shown in Table 3. 
The differential equation itself has no finite solution when 6 = — 12.32, whilst 


the difference equation has no finite solution at the following values of 6 and N 


nS 10 12 15 31 63 127 10 
6 +1.95 +66.98 —85.28 -—28.61 -—14.40 -—12.79 -—12.44 —12.32 


It is the lack of correspondence between critical conditions of the differential and 
difference equations which causes poor agreement to exist between the exact so- 
lutions. For example, when »/kL = 4, 6 = —15, the difference solution at x = L 
is —.1180 for N = 3, and as N tends to 29, the solution increases monotonically 
to +. For N > 29, the solution increases from — * reaching the value — .9034 
when N = 63. Thus although the mesh length h( =+/kH) is as small as js, the 
exact solution of the difference equation is approximately 20% less than the exact 
solution of the differential equation. Of course, for values of 5 less than —15, the 
infinite step in the solution of the difference equation occurs at a value of N less 
than 29, (see Table 3) and so the solution of the difference equation for N = 63 
is correspondingly more accurate. When 6 = —30, for example, the error is re- 
duced to 2.65%, which nevertheless is still appreciably large. 


TABLE 3. Values of xy 














| ra 
hee 2 a 
| = 12 —12.32 | 35 —18 mH a ~i.. | =o =—7 | =<# 
| | } 
a rR Resi’ VND oBt 3 ‘ 
3 | —.1826 | —.1433 | = 1180 | — .1002 a —.0871 | —.0771 | —.0691 | —.0626 
10 | —.0263 | —.0253 | —.0244 | —.0235 | —.0227 | —.0220 | —.0213 | —.0206 
12 | +.0262 | +.0273 | | + .0284 | +0297 | +0311 | +..0326 | 4.0343 | +.,0362 
+ @ (14) 
15 | +.1020 | +.1204 | | +.1469 | +.1885 | +.2628 | +.4338 |+1.2420 |—1.4391 
| | | + © (29)| +e = (2) + © (20)) + © (18)] + = (16) 
31 | +.3706 | +.8343 | |-3.3183 | —.5551 | —.3029 | —.2083 | —.1587 | —.1282 
63 | +.5282 +2.5442 | | —.9034 | — “3036 | — 2435 | —.1784 | —.1407 | —.1162 
@ | +.6016 |+6.1623 | + —.7476 | — - 


.9024 | — .2305 1713 | —.1363 | —.1132 
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Fic. 2.—Comparison of finite difference solution with solution of differential equation 


for varying 6. 


Two figures are now giver to illustrate some of the results contained in Table 3. 
Fig. 1 illustrates the manner in which the exact solution of the difference equation 
converges towards the exact solution of the differential equation as the number of 
internal nodes is increased. The outstanding feature is the infinite step in the 
solution. Reduction of the mesh length only improves the solution of the difference 
equation once the infinite step has been passed. Fig. 2 compares the exact solutions 
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of the difference equation when h = %, 7g respectively with the exact solution of 
the differential equation, for a range of values of 6. The solutions are infinite at 
6 = —14.40, —12.79, and —12.32 respectively, and the diagram clearly illustrates 
how this difference in the critical value of 6 causes poor agreement to exist between 
the exact solutions over a wide range of 6. 


5. Concluding remarks. Although in the present note, no attempt has been 
made at an exhaustive study of the accuracy of finite difference solutions of two 
point boundary problems, sufficient examples have been given to show the lack of 
agreement which can exist between exact solutions of a differential equation and 
its associated central difference replacement. It is also unlikely that conditions are 
improved by using alternative difference schemes either of the same or of higher 
order, since critical conditions of the differential equation have their counterparts 
in all difference replacements. In fact the use of higher order replacements may well 
complicate the situation further by introducing additional critical conditions into 
the system. This was certainly found to be the case in step-by-step difference solu- 
tions of initial value problems by Todd [1]. 

An important point illustrated by the results in Fig. 1 is that in order 
to guarantee a reasonable answer using difference methods it is necessary but by 
no means sufficient for the exact solution of the difference equation to tend to the 
exact solution of the differential equation as the mesh length tends to zero. In fact, 
the computor is well advised to obtain solutions of the difference equation for a 
range of mesh lengths in order to detect any infinite step in the solution of the 
difference equation. For example, in Fig. 1, convergence occurs only for h < +%. 

Finally, with regard to the accuracy of the numerical results presented in the 
present paper, it should be stated that all calculations were carried out on a 10 x 10 
electric desk machine, the latter being used to its full capacity in all calculations. 
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Expansion of the Confluent Hypergeometric > 
Function in Series of Bessel Functions 


By Yudell L. Luke 


Abstract. An expansion of the confluent hypergeometric function #(a, c, z) in 
series of functions of the same kind has been given by Buchholz [1]. By specializa- 
tion of some quantities, there is obtained an expansion in series of modified Bessel 
functions of the first kind, I,(z), where v depends on the parameter a. Tricomi [2, 3] 
has developed two expansions of similar type where both the order and argument 
of the Bessel functions depend on the parameters a and c. In the present paper, 
we derive an expansion in series of Bessel functions of integral order whose argu- 
ment is independent of a and c. Our expansion is advantageous for many purposes 
of computation since the parameters and variable of (a, c, z) appear in separated 
form. Also, for desk calculation, extensive tables of J,,(z) are available, while for 
automatic computation Bessel functions are easy to generate [4]. 

Special cases of the confluent function, such as the incomplete gamma function, 
are also studied. For the class of transcendents known as the error functions, 
including the Fresnel integrals, it is shown that our expansion coincides with that 
of Buchholz [1]. By specializing a parameter and passing to a limit, we derive ex- 
pansions for the exponentiai integral and related functions. Other expansions for 
the error and exponential integrals are derived on altogether different bases. Finally, 
some numerical examples are presented to manifest the efficiency of our formulas. 


1. Representation of the Confluent Hypergeometric Function in Series of Bessel 
Functions of Integral Order. The confluent hypergeometric function is defined by 
the series 
. (a),2" 


O(a,c,z) = iF ,(a;¢;z) = 4 
a y k=0 (c).k! 





(1.1) 
a ad l(a + k) 
k (a) ’ 
where I'(z) is the gamma function. We also have the integral representation 
I'(c) [ zt,a—1 c—a—1 

a = —_ t 1-—t t (a) 
(12) &(a,c,z) T@re —a) + e ( ) dt, R(c) > R(a) > 0, 
and Kummer’s formula 
(1.3) @(a, c,z) = e@(c — a,c, —2Z).. 


For these and other properties of the confluent function, see [5). 
Let ¢ = sin *@ in the integrand of (1,2). Using the Jacobi expansion [6] 


(1.4) e(!200s2 _ > « cos 2k0l,(z/2), o = 1, « =2, k >0, 
k=0 
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(1.2) becomes 


Sa &(a,c,z) = Wat aj > ns I,(z/2) 
. [cos 2k0(sin 0)** "(cos 6)" de, 
Since 
\* k 4 ‘ 2m 
(1.6) sos wid we (~1)? > (—1) (*) (k)m(cos 6) 





m=0 (1/2) m P 


4 


>. , : , 
where e is the binomial coefficient, and 


ieney [ (sin 0)** "(cos 6)” do = 1 P(z)P(y) 
0 


2T(a+y) ’ R(z)>0, Ry) >0, 





we find that 
(1.8) @(a, c,z) = ce”? > (—1)*aRi(a, c)1,(2/2), 
k=0 


where R;(a, c) is conveniently expressed in hypergeometric form [7] as 
(1.9) Ri(a, c) = sF2(—k, k, a; c, 3; 1). 


An alternative expression for R,(a, c) is readily reduced from a result in [8]. We 
have - 


Me) TK = ©) p.(a,c) = G(a,e) + (—1)'Gy(e — a,c), 
I'(c) 
k . =s ae ed aaa 
(1.10) idiad x (—1) cos (¢ a)a0(h c + 1)T (2c. 2a) 





3-2IT(k + 1 fe — 2a) 
h(t — 2a,2c — 2a;k+1+e¢e— 2a; 5). 


If a and ¢ are fixed and k is sufficiently large, then the »/;’s in (1.10) are of 
order unity and 








Lae — k on _ ~ ee 2a—2c 
T(a)T(e a) R, (a,c) = (—1)* cos (ec e)s0'(3e 2a)k*~ 
T'(c) Dre 2a—1 
(437) P 
cos arI'(2a)k “ 
+ 7 
Here we have used the fact that for fixed a and 8, 
, T(k +. a) a—s 
1.12 ——____ ry kh", k large. 
1.42) T(e + B) . 





> 0, 


We 


2a 











EXPANSION OF THE CONFLUENT HYPERGEOMETRIC FUNCTION 
Again, if z is fixed and k is large, then 


(1.13) I,(z) & 





(z/2)* 
i 
It follows that (1.8) converges like 

> (Ak*~ + Bk) z* 


Dk je | ’ 





(1.14) 


k=1 
A and B are constants. A similar argument shows that (1.1) converges like 


rs ko ck 


k=1 k! ‘ 


(1.15) 


and so (1.8) converges more rapidly than (1.1). 
Some general properties of R,(a, c) are next of interest. The combination of 
(1.3) and (1.8) shows that 


(1.16) R.(a,c) = (—1)*R,(e — a,c); R,(a, a) = (—1)*. 
Further relations follow by applying known contiguous relations of ®(a, c, z) to 
(1.8). For example, since 
(1.17) $ Hae) = 844+ 1,c + 1,2), 
dz c 
it follows that 


(118) R(a+1,¢+1) = —— [Riss(a, ¢) — 2R,(a, c) + R,-(a, c)], 


and 

(1.19) (@a—c+1)R,(a,c) — aR,.(a + 1,c) + (ce — 1)Ri(a,ce — 1) = 0, 
results from a contiguous formula of the form (1.19) with R,(a, c) replaced by 
(a, c, z). To obtain a pure recursion formula for R,(a, c), apply the differential 
equation satisfied by (a, c, z) to (1.8). Employ the difference-differential prop- 
erties of J,(z) and equate to zero the coefficient of J,(z). We get 

(1.20) (k+ c)Risi(a,c) = 2(c — 2a)R,(a, c) + (k — c)Ry4(a, c). 


For certain values of the parameters, the coefficients R,(a, c) can be expressed 
in simple form. Thus from (1.10), 
[T(c)} 
(ec + k)I(e — k)’ 





(121) R,(4,c) = c ~ 0, —1, —2,---, R(c) > 0. 


This also follows from (1.9) and the known expression for a »F; with unit argument. 
If c = 2a, then from (1.9) and (1.20), R,(a, 2a) is zero if k is odd, and 
(3 — a). 


(3 + a),’ 


1 3 5 
ax =—§, “aah se te 





(1.22) Rx.(a, 2a) = 


Since 


(1.23) (vy + 3, 2v + 1, —z) = M(v + 1)(2/4) "1, (2/2), 
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it follows that 


(—1)*e, In (2) 
—-k+1)lot+kt+ 1) 


provided » is not a negative integer The latter is a special case of a known result [9]. 


(1.24) I,(z) = (z2/2)’T(v + 1) oR 





2. The Incomplete Gamma Function. If c = a + 1, we write 


(2.1) y(a,z) = a ze @(1,1+a,z) = a '26(a,a+ 1, —2), 


(2.2) y(a,z) = [et “2 dt = T(a) — é et dt, R(a) > 0. 


(2.3) y(a,z) =a 2c" > & R,(a)1i(z/2) 
k=0 


R,(a) = R,(a, a + 1) sF'2(—k, k, a,a + a, 1) 


(—1)';F2(—k, k, l;a + 1, 3; 1) 
—(-1)*a 


= ao = 9-}- — 9.1 
HE — ak — at 1p eit — 2% Ak — @ + 254) 


I 


Il 











(2.4) 
4 ao anI'(2a)kT'(k — a) 
P-Tk+at+i1) ’ 
and 
(2.5) (k+a+ 1)Rey(a) = 2(1 — a)Ri(a) + (k_— a — 1)Riur(a), 
1 
(26) (a +1) = —SEY ira) — 2h (a) + Rea(a)) 


If a = 3, we deduce the following expressions for the error function and related 
functions. 





s e —z2/2 
27 — feta a 2 (oe /2) 
(2.7) Erf z l e dt 2 ai 
(2.8) ec“ Eric =" [ the «> ee T12/2) | 
, 0 k= 4k? — 1 


29) (2n)'C(2) + i8(e)] = [tet at = ayes Soaeae 


A more general expansion for the error function can be found as follows. Now 


zsin@ 8 

(2.10) [ edt = ze” [ e'**/) 00838 aos ¢ dt, 
0 0 

and using (1.4), we obtain 


, _2e sin (2k +1)6 , sin(2k —1)@| —.20,, . 
(2.11) Erf (zsin@) = sda 1 - pm \ T,(2°/2), 








; [9]. 
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Nie 
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EXPANSION 





OF THE CONFLUENT HYPERGEOMETRIC FUNCTION 


and from the latter follows the interesting representation 


(2.12) [ Erf (z sin 6) sin (2k + 1) 0 d0 = 
0 


Tz 


eal 21a) 
4Qk +1)° (T,(2’/2) + Ts1(2°/2)). 


3. The Exponential Integral. This is a limiting form of the incomplete gamma 
function. From (2.2), we have 


(3.1) 


| te dt = lim {[(a) — y(a,2z)}. 


a+0 


Another useful form is 


(3.2) 


“0 


[« 


(1 — e") dt = lim {a 2’ — y(a,z)} 
a+0 


Using (1.4) and (2.3), we write 


> ieee 19 ieee Te 
(3.3) cg | eu" du = tn + 1) 1\ fF + 25541 — Blo) 1,{z/2). 
z ) k=1 | 





a a 


Application of L’Hospital’s theorem gives 


(3.4) 


lim {[T(a + 1)z2°% — 1\ 





a0 


= —(y+I1nz) 


a 


where y is Euler’s constant. Define 


(3.5) 


Then 


(3.6) 


he = 


=. ee 
1 — R,(a) k (—1) ¥ (k) mn 
~ Sere Sat 





—FEi(—z) + (y + Inz) = / te * dt + (y + Inz) 


In view of (2.5) and 


(3.7) 


- [ (1 —eé*)d =2> > fe" I, (2/2). 
0 k=1 
(3.5), we have the recurrence formula 


(k + 1) fess = fi + (k La 1) fra 7 4, 


and by induction or otherwise, we can prove that 


(3.8) 


< 





ba as — fi + 


7 -(-DH +440 +(-1)}, fr=0,f, = 2, 


1 1 
ae +5"), 


fonss so fox + 


Employing the theory of the ¥(z) function, 


2 


% +1 ? Forse li Fon+s +=; 


Ts 


the logarithmic derivative of the 
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gamma function [10], we find that 

(3.9) lim (fi — 2In 2k) = 2y, 

and so (3.6) converges like 

(3.10) } Uz’, tu = (Ink)/2™k! 
An alternative form of (3.6) is 


(feu-oe 
{ 0 
| 


II 


2[1 — e*"Iy(2/2)] + 4 > Gis2€ "Th 42(2/2), 


(3.11) 

(ge = 3h — 1 
Now 
(3.12) Fi(z) = [ te! dt = Bi(ze**) + ix, 
and so 

e te’ dt = &*(y + Inz) + Qe?" Iy(z/2) — &*] 
(3.13) ba 


— 4 a (-1 ) gorse e?T40(z/2) 


k=0 
We also have 


[ t*(1 — &") dt = (y + Inz) — Cifz) + iSi(z) 
(3.14) 


= 2? S* f, 7 J.(2/2). 
k=1 


4. A Second Expansion and Special Cases. Buchholz [1] has given the expansion 
(a, ¢,2z) = e""(2/4)*“T(a — 3) 


(4.1) $e (k+a-HQa- 1, 
k=0 k\(e), 





(2a =e 7, (2/2), -R(a) > 0, 


which coalesces with (1.8) when a = }. If a, c and z are fixed, then an analysis 
similar to (1.11)—(1.14) shows that (4.1) converges like 


oo | saliiendis 2k 


a ’ 
tm Ss 2**K! 


(4.2) 


and from this point of view, the convergence of (1.8) and (4.1) are nearly alike 
when R(c — a) > 0. 

We now deduce from (4.1) series expansions for the exponential integral and 
related functions. The method of proof is akin to that in the previous section. 
If in (4.1), we put a = c, there follows a representation for e”” since ®(a, a, z) = e’. 
Use this and (4.1) with c = a + 1 and z replaced by —z. Then after the manner 
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of (3.3), we have 
a a 
fas) /2 —l1 —t 1/2 ‘ 
(4.3) é / te dt = (xz) » hy. (a) In4a1/2(2/2), 
z k=0 


and employing the duplication formula for gamma functions, the elements h,(a) 
take the form 


2k 4+ @ — 1/2)P (§ SS a) r (f + a) 
(4.4) h(a) = Tk! 





{a Ss. 
Tr(a — 1k + ak +a — 1){° 


For k = 0 and 1, application of L’Hospital’s theorem gives 





(4.5) h(0) = —3(7¥ + Inz), hi(O) = ho(O) + 1. 
Otherwise, 
(4.6) h(0) = 2% Sy pho k 
k(k — 1) 
Thus 


—Ei(—z) + (y+ Inz) = | te‘ dt + (y+Inz) = [ (i —eé“)dt 
z 0 
(4.7) en 
Ms 1/2 —2/27 1/2 —_ Ok+3) _ 
= (42) 6" Tya(2/2) + (xz) > EL DEPD 
and forms for the related functions follow readily enough after the manner of 
(3.12)-(3.14). By a familiar argument, (4.7) converges like 


—2/2 7 
e 'Iissn(z/2), 





(4.8) > Dp 2", “N= [2"*kT'(k + >) 


k=1 


and so (4.7) converges more rapidly than (3.6). However, the difference is not 
great since u;/v, — k* In k (see (3.10) ). 


5. Further Expansions. The following representation is due to Tricomi [2, 3]. 


(5.1) (a,c, z) = T(c)e” 8 A,(a, c, h)2*Ec44-1(—az), h = 0, 
k=0 
(5.2) E,(z) = 2° J,(22), E,(—z) = z”"I,(22), 
where the coefficients A;,(a, c, h) are defined by the generating function 
(5.3) > Ar(a, c, h)e* = & [1 + (h — 1)zy*(1 + hz)”, z| <1, 
k=0 


and satisfy the recurrence system 
(5.4) (k + 1)Agus = [(1 — 2h)k — helA; 
+ [a(1 — 2h) — h(h — 1)(0 +k — 1)JAeu — h(h — 1)Ar-, 
Ag=1, Ar:=-—he, As = }W'e(c +1) + a(} — A). 
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A second expansion given by Tricomi [2, 3] is 


(5.5) (a,c,2z) = T(c)e”” >> A,*(K, c/2)(2/2)"Beu (Kz), K = ° 
k=0 


=<-a, 

where 

(5.6) > AM(K, c/2)2z° = &*(1 — z)* "71 + 2), z| <1, 
(k + 1)Aty = (k +c — 1)Any — 2KAb oe, 

(5.7) 


Both (5.1) and (5.5) are useful to study the behavior of when the parameters 
are large [11], but are not too suitable for many purposes of computation, since 
the order and argument of the Bessel functions depend on the parameters a and ec. 

In (5.1), put a = 1, replace c by a + 1 and set h = 0. Then 


(58) y(a,z) = (ae “2? D ex(—1)2"Fe4a(22"”), e(—1) = YS 
k=0 m=) L: 


follows from (2.1) and (5.3). The method of proof surrounding (1.11)-(1.14) 
shows that (5.8) converges like > ro z'/T(k + a + 1) which is the same as (1.1) 
for the above selection of parameters. 


6. Further Expansions for the Error and Exponential Integrals. We start with 
the representation 


2*T (: a hay ? 
[ I(t) dt = — 2 
0 


fy—ptil 
Crary 














(6.1) 
: _ 1 
ss (y+ 2k + 1)r (5 ~+k 
- >) (-1)* —+__= f Tymii(z), Re +v+1)>0 
dX ( ) “fy t+uts3 +241 (2) (se tet 
I ‘ + k 
2 
which is readily verified by differentiation. Put 1» = 0, and consider the formulas 
when vy = 3 and vy = —3. Then by subtraction and addition, we get 
(6.2) Erf z = (/2)' > (-1)"""7,_4(2). 
k=1 
i ~ 4 9 
(6.3) Erfi z = (/2)' >> (-1)"" yaa(2’). 
k=0 


Here [k/2] is the largest integer, including zero, contained in k/2. The latter two 
equations were also verified by Tricomi [3] by means of the Laplace transform*. 





* There is a typographical error in [12] for (6.2). There the summation should start with 
n=1,notn=0. 
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Another representation which follows from (5.8) is 
1 2 = 
. . z k ‘ 
(6.4) Erf z = 3(xz)’e* > ex( —1)2*In44(2z). 
k=O 


By a familiar argument, the convergence of (2.7) is superior to that of (6.2)—(6.4). 

To obtain another expansion for functions related to the exponential integral, 
it is convenient to again use (6.1) with z replaced by iz. The ensuing formula is 
of the same type as (6.1) with J,(z) replaced by J,(z) and the factor (—1)* omitted 
behind the summation sign. Now put 4» = —43 and v = 3. Then for the sin in- 
tegral, we have 


ne sin ¢ r\'" = (4k + 3)k! 
(65) [ ot at = (x) > 


=~ (Be 











Similarly, with » = —4, we can form an expression for 


[ t’*?(1 — cost) dt. Let wp — — 3; then 


“0 


1) 
soe) Pe (1 — cost) dt = 3(Z) yet DG) 


A(k + 1) — WL) + Wk + 4) — W4)} Joeaalz). 


Our previous analysis shows that for z fixed and k large, the ratio of the kth term 
in the expansion (6.6) to the kth term of the Taylor series expansion for the cosine 
integral is proportional to k(In k)/2™. Similarly for (6.5), the ratio is proportional 
to k/2™. If (3.6) is compared with its corresponding Taylor series representation 
in like fashion, we again obtain the former ratio. Thus the proper combination of 
(6.5) and (6.6) gives about the same convergence as (3.6) and (4.7). 

As a remark aside, using (6.1), we can produce 





lim { [1 — Jo(t)] at 
po—l 0 


[ eu- nwa = +z) 4 [ roa 
o 0 


ll 


22>) (2k + 3)[v(k + 2) — (1) JJx4s(2) 


k=O 
=1—2 (2) + 2°D (2k +5 [Wik + 3) — WL) — 1] Jxe4s(z). 


For tables of this integral when z is real, see [13]. 
We now derive yet another expression for the exponential integral. It is known 
[14, 15, 16]* that 


* We take this occasion to point out some typographical errors in [15] and [16]. The formula 
(dJ,(z)/dv),.. should read 


aJ,(z)| **Jn(z) , 
WA) — Gals) + Ul > ‘S , Gre) = —5 Yale). 
Ov |vak =o ni(k — ad 
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(6.8) EFi(—22) = (ne/2)"e" | sts) 4 91,(2) .: 
O, = |vm/2 Oy |ra—1/2 


where Ei(—z) is given by (3.6). To evaluate the partial derivatives in (6.8), 
we use the expansion [17] 
(2/2)’ & (—1)*2"1k(z) 


(6.9) Ie) = “Tey 2 GE BEL? 





and find that 


(6.10) [ (1 — &") dt = 2e°*” cosh z/2 + 2g;(z) —2g2(z), 
0 


(6.11) e- [ t(e' — 1) dt = —2e*” cosh 2/2 + zg,(z) + 2go(z) , 
0 


where 


> ( 1)"(z/4)*e* "1, (2/2) 
6.12 lm = eee oe de € _ 2 
se On(2) to) ie (2K + 3 — 2m)? . waitin 


Of all the exponential integral expressions given, this converges the best. 


7. Numerical Examples.* 
EXAMPLE |. Using (1.8) and (1.9), we illustrate the computation of 


(7.1) e *b(m + 1 — ia, 2m + 2, 2ixr) = Do aR Ji(2), 
k=0 
(7.2) R, = 3F2(—k, k, m + 1 + ta; 2m + 2, 4;1). 


Suppose a = 1, m = 2. Then values of 7R, are easily generated using (7.2) and 
(1.20). We have Ry = 1, iR, = 4, 7R, =4%, etc. Using standard tables of Bessel 
functions, if « = 1, seven terms of (7.1) give the six-decimal value 1.451140. If 
x = 2, nine terms give 1.293748. For tables of (7.1), see [18]. If « = 2, about 15 
terms of the Taylor series are required to achieve six-decimal accuracy. 

EXAMPLE 2. We employ (2.7) and (2.8) to compute the error integrals for 
z = 2. Using nine terms, we find Erf 2 = 0.88208 14 and e “Erfi 2 = 0.30134 04. 
About 19 terms of the Taylor series expansion are needed for the same accuracy. 

EXAMPLE 3. To illustrate computation of the exponential integral, it is suffi- 
cient to consider (6.10). If 2 = 4 and we use six terms of each series in (6.12), 
then the value of the integral (6.10) is 1.96728 94. Similarly, if z = 8, only eight 
terms of each series in (6.12) are needed to yield the value 2.65669 49. For z = 8, 
about 32 terms of the Taylor series expansion would be required to achieve the 
same accuracy. 


* The author acknowledges the assistance of Sue Chapman and Betty Kahn for the de- 
velopment of numerical examples. 
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Quadratic Residues and the Distribution 
of Primes 


By Daniel Shanks 


1. Introduction. Chebyshev stated, [1], that there were “many more” primes 
of the form 4m — 1 than of the form 4m + 1 and he indicated, [2], at least three 
senses in which this assertion was to be understood. If the numbers of primes 
in these two classes which are <n are designated as r_(n) and (nm) respectively, 
and if 


(1) A(n) = m_(n) — r4(n), 
and 
A(n) 
(2) t(n) = ————_, 
/n/log n 
then Chebyshev stated that for a certain subsequence, n; , of the integers, 
(3) Lim t(n;) = +1. 


Since+/n/log n — ~, x(n) could thus be made to exceed r4(n) by any amount 
by an appropriate choice of n. 

This theorem (which was one of the above-mentioned three senses) was proven 
by Phragmén, [3], and later by Landau, [4], [5]. The other two “senses’’ have never 
been proven, [5, p. 647]. Since the theorem, (3), is entirely noncommittal as to 
the behaviour of A(n), or of t(n), for values of n other than those of the subse- 
quence n; , it is however not a very convincing argument in favor of Chebyshev’s 
“many more’’ primes assertion. Littlewood, in fact, showed, [6], that there is a 
constant K such that for two other sequences n; and n, we could have 


(4) t(n;) > K log log log n; 
and 
(5) t(m) < —K log log log n, . 


Since log log log n — ~ (with great dignity) these inequalities imply that the 
right side of (3) may be changed from +1 to any real number, positive or negative, 
for some sequence, n; . Because of this Cramér wrote, [7]: 

‘‘Chebyshev’s general assertion—‘there are many more primes of the form 
4n + 3 than of the form 4n + 1’—can therefore be true only in a rather limited 
region.”’ 

The purpose of this paper, [8], is to determine the sense in which Chebyshev’s 
general assertion is (nonetheless) probably true, and to identify the nature of 
the cause which tends to diminish the class of 4n + 1 primes relative to the 4n — 1 
primes. The problem is generalized to primes, and also to composites, in any arith- 
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metic progression. Finally, there are some comments concerning the Littlewood 
phenomena, r(n) > Li(n). 


2. The Function A(n). Tables of x(n) and of 2,(n), [9], [10], [11], [12], pub- 
lished prior to the short note of John Leech, [13], were not of sufficient completeness 
(or, in one ease, of sufficient accuracy) to correctly identify Cramér’s ‘‘rather 
limited region.’’ Leech computed (on EDSAC) A(n) for n S 3,000,000 and de- 
termined that 26,861 was the smallest* n for which r,(n) > 2_(n): 


(6) A(26,861) = —1. 


The author, unaware of Leech’s note, carried out a similar computation (on an 
IBM 704). The following description of A(n) in this range, 1 to 3-10°, agrees 
(where it overlaps) with the shorter account of Leech. In Table 1 we divide this 
range into six regions and in each region give the maximum and minimum values 
of A(n), the number of intervals during which A(n) remains zero, positive, or 
negative, and the total number of n’s for which A(n) is zero, positive, or negative. 
The final figures show that A(n) > 0 for 99.84% of the n < 3- 10°. Table 2 identifies 
some extrema (and some zeros) of A(n). 

This detailed description makes it highly plausible that the predominantly 
positive character of A(n) in this range of n is not merely a passing fancy of the 
integers (as is almost implied by Cramér’s remark) but a permanent phenomenon 


TABLE | 
Description of A (n) for 1 = n S 3,000,000 


Regions of » Max Min a + Int.'— Int. 0 n’s + n’s —n's 
1 - 462 +6 0 5 4 0 10 452 0 
463 — 26 ,832 +31) +1 0 ] 0 0 26 ,370 0 
26,833 — 26 , 926 +1); -—1 6 4 l 60 32 2 
26,927 — 616,768 +105) +1 0 l 0 0 589 ,842 0 
616,769 — 633 , 882 +12) —8 | 101 | 48 | 52 | 1,282 12,428 3,404 
l 0 0 2,366,118 0 


633,883 — 3,000 000 +|+256 +1 0 


~ 


1,352 (2,995,242) 3,406 
00.05% 99.84% (00.11% 


TABLE 2 
Some Special Values 


A(227) = +6 
A(461) = A(462) = 0 
A(17,971) = +31 


A(26 ,861) = A(26,862) = —1 
A(359 ,327) = +105 
A(623 ,681) = —8 
A(627 ,859) = --- = A(627,900) = 0 


A(2,951 ,071) = +256 


* This first axis crossing of A(n) was discovered independently by J. W. Wrench, Jr. 
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for which a sufficient number-theoretic cause should be assigned and of which a 
more precise formulation is desirable. 


3. The Function 7(n). We therefore return to ¢(n) and seek a better insight into 
this function. We do not compute ¢(n) itself, but instead, the somewhat simpler: 
(7) r(n) = A(n)-~/n/x(n) 


where 2(7) is the total number of primes <n. From the prime number theorem, 
a(n) ~ n/log n, we have 


(8) tT(n) ~ t(n). 

For example, for n = 2,000,000 we have 
a— (n) = 74,516; wi(n) = 74,416; a(n) = 148,933; 
A(n) = 100; t(n) = 0.9496; and t(n) = 1.0259. 


The function 7(1000 k) was computed for the 2000 values k = 1, 2, --- , 2000. 


9 =) 


The minimum, the mean, and the maximum of these 2000 values are 


7(629,000) = —0.0464, 
2000 


> 7(1000k) = 1.0613, 


c oe 
(9) 2000 “F 


7(127,000) = 2.0961, 


respectively. Only slightly greater extremes would have been obtained had we 
computed r(n) for all the 2-10° n in this range. The distribution of these 2000 
values of 7 between the extremes was determined by counting the k’s for which 


(10) — < 7(1000k) < ———__ (m = —1,0,1,--- , 33). 


These counts, v(m), are tabulated in Table 3 and plotted in a bar graph in Fig. 1. 

The following comments are now in order: 

a.) The distribution is roughly normal with a mean of (nearly) +1. 

b.) The rare cases of A(n) < 0 (i.e., 7(n) < 0) are now to be thought of as no 
more unusual than the equally rare other extreme: A(n) > 2x(n)/+/n (ie., 
t(n) > 2). 

c.) The implication of Littlewood’s inequalities, (4) and (5), is that the distri- 
bution function has tails of infinite extent. Of course, occurrences of 7 far out in a 
tail will be very rare. 








TABLE 3 
awh! £1 £} 2in8).4) 6] 6| 7] 8] O] » 
v(m) | 4) 7) 9] 7] 10) 4) 38] 14] 11] 48/| 55) 80 
Ps ne oh ee salt Se TE re Ue” ee ee eee eee eee eee 
m 11| 12| 13] 14 | i5| 16! 17| 18| 19} 20| 21 
| 


1 
v(m) 76 | 114 | 139 | 15 


m 23 | 24 | 
v(m) 52 | 51 


| 151 | 152 | 136 | 129 | 109 , 107 





25| 26| 27| 28] 29| 30| 31] 32] 33 
299| 21| 19] 17| 12| 5| O| 2 











eee ~~ | 





QUADRATIC RESIDUES AND THE DISTRIBUTION OF PRIMES 275 


~8 C) 8 16 24 32 40 
250 — ——— 


+ ‘ 
v(m) is the number of k's, such 
thet | <= k = 2000 and 


M™<« k m+i 
| 16 <= r(i00o0k) < ie 
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Fig. 1—The distribution of 2000 values of 7(1000 k) 


d.) Chebyshev’s function +/n/log n (or its equivalent r(n)/+/n) is seen to 
be an appropriate normalizing factor—it reduces the function A(n) to a function of 
n with relatively small variation. For example, throughout most of the second 
million, (1,015,000 to 2,000,000), 7(1000k) has the following very modest be- 
haviour: 


(11) 7(1,811,000) = 0.426 <= 7(1000k) < 1.610 = 7(1,521,000). 
e.) The function 
1 = 
? oe 1000k (gs = 1,2,--- ,40 
T 50s i 7( ) (s ) 


provides a running history of the mean value of 7, and is found to change very 
little: 


(12) M5 = 1.0222 S 7, S 1.1799 = 7; (s = 1,2, ---, 4). 
4. A Conjecture. The above discussion suggests that the proper formulation 
of Chebyshev’s ‘‘many more”’ primes assertion is the following 


CoNnJECTURE. The mean value of 7(n), 


N 





1 
wai 7): 
has a limit as N — ~ and this limit equals +1; 
N 
(13) Lim > r(n) = +1. 
N+o N —_ 1 2 


The conjectured limit, +1, may seem a little rash in view of the limited data, 
(12). However, there is other evidence. The remark made after (5) may suggest 
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that the +1 right side of (3) is of no significance since it may be replaced by any 
real number. Nonetheless, from Landau’s proof, [5], of (3) it is readily seen that 
+1 has a very special rele. For every a the Dirichlet series 


(14) F.(s) = > A(n) — av/n/log n 


neti 





n=2 


defines a function which is regular in the half plane ®(s) > 1 and for the real 
values } < s S 1, but is singular for s = }. However, only for a = 1, that is, only 
for F,(s), does the function have a limit as s — 4+. Using this fact, and r(n) ~ 
i(n), what is to be proved is that the mean value of the a, in 


i) 


An 
14a) F,(s) = —_— 
(14a) (8) = i; 
has a limit as n — © which is equal to zero. 
The conjecture has not been proved, either as its stands, or in one of the two 
following weaker forms: 


ConJeEcTURE. If no zero of 


1 1 ] 
ee 


& 7 & 


L(s) =1-5+ 


has a real part >4, then (13) is true. 


ConJEcCTURE. If 
N 


Lim Z r(n) 


N+o N — 1 2 





exists, it equals +1. 

However, even the original conjecture seems sufficiently likely to be true to 
merit recording. The closest thing to it, in the literature known to the author, 
is Theorem 2.351 of Hardy and Littlewood, [6, p. 151]. This, together with their 
Theorem 2.34 and a third result due to Landau, [7], may be combined and entirely 
rewritten to read: 

THeEoreEM. If no zero of L(s) has a real part >4, then there is a positive K 
such that 


| 


} > A(n) > K VN 
+, opty log N 


=] 


for all sufficiently large N, and conversely. 

From this theorem, Abel’s lemma on partial summation, [14], and (8), we 
can easily prove the following 

Coro.uary. If no zero of L(s) has a real part >4, then there is a positive 
K such that 
1 
N-1 





¥ en) >K 


for all sufficiently large N. 
We will present further evidence for the conjecture, in section 9 below, after 
we have determined the number-theoretic cause of the 4m + 1 prime deficiency. 
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5. Localization of the Deficiency. If the classes 4m + 1 and 4m — 1 are split, 
modulo 8, we have 


4m+1 4m — 1 
8m + 1 8m +5 8m +3 8m +7 








and we now inquire whether the deficiency exhibited by the 4m + 1 primes is 
shared equally by the 8m + 1 and the 8m + 5 primes. Landau’s generalization of 
(3) for any modulus, [5, p. 704-711], and other simpler considerations given below, 
suggest, on the contrary, that the deficiency is confined entirely to 8m + 1. 

Let 72(m) be the number of primes of the form am + 6b which are <n. A 
table of 23,(1000k) was computed for b = 1, 3, 5, and 7 and k = 1, 2, --- , 1000. 
Examination of this table shows that while the numbers of primes of the forms 
8m + 3, 8m + 5, and 8m + 7 are continually jockeying for first place, the number 
of 8m + 1 primes is always in last place. Let f,(n) be the fraction of the tabular 
values, with argument <n, for which am + b is in first place. (Two-way or three- 
way ties for first place are pro-rated 4 or 3 of a point respectively.) Similarly let 
h(n) be the fraction of the tabular values for which am + 6 is in last place. In 
the present case, since 8m + 1 is in last place throughout the entire range of the 
table, we have for every k = 1, 2, --- , 1000: 


(15) f:(1000k) = 1,(1000k) = 1,(1000k) = 1,(1000k) = 0. 


A brief summary of zs, and the four other f, and J, functions is shown in Table 4. 
No significant deviations from equality, [15], are to be noted among 8m + 3, 
8m + 5, and 8m + 7,— the slightly low f,(1,000,000) is compensated by the 
slightly high f;(750,000). We conclude that the deficiency of the 4m + 1 primes 
resides solely in the 8m + 1 primes. On those rare occasions (regions 3 and 5 of 
Table 1) when 24(n) > 2_(n), we invariably find that 8m + 5 is in first place 
with a sufficient lead over 8m + 3 and 8m + 7 to overcompensate for the leads 
which these latter classes hold over 8m + 1. It is not known whether 8m + 1 
itself can ever take the lead but in view of what is written below about 24m + 1, 
this is probable — for sufficiently large n. 


TABLE 4 
Distribution of the Primes Modulo 8 





i 
ts i 





n 78,1 rs, 3 | 78,5 8,7 h | js 
} | | 2 
250 ,000 5473 5525 5523 | 5522 | 1.000 | 0.446 | 0.322 | 0.232 
500 ,000 10334 | 10418 | 10397 | 10388 | 1.000 | 0.337 | 0.344 0.319 
750 ,000 14998 | 15087 | 15078 | 15074 | 1.000 | 0.380 | 0.382 | 0.237 
1,000 ,000 | 19552 | 19653 | 19623 | 19669 | 1.000 | 0.401 0.287 | 0.312 





6. Distinctions between Modulo 12 and Modulo 10. Alternatively, we could 
have split 4m + 1 and 4m — 1 modulo 12: 


4m + 1 4m — 1 
12m + 1 12m + 5 12m + 7 12m + 11 
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TABLE 5 
Distribution of the Primes Modulo 12 

















n wi2, 1 | wie, 5 | wie, 7 m2, u h ts fr fu 
relic | “a = ert E wk —- 
252 ,000 5504 5567 | 5564 | 5566 | 1.000 | 0.352 | 0.420 | 0.228 
504 ,000 10404 | 10476 | 10480 | 10472 | 1.000 | 0.231 | 0.274 | 0.495 
756,000 | 15100 | 15196 | 15204 | 15186 | 1.000 | 0.396 | 0.258 | 0.346 
1,008 ,000 19715 19771 | 19812 | 19797 1.000 | 0.338 | 0.402 | 0.260 





TABLE 6 
Distribution of the Primes Modulo 10 








| 








n | 10, 1 10, 3 10, 7 | 10, 9 | h | fs ir lo 
— } . a a a saa a — ——— 
250 ,000 | 5495 5520 | 5517 5510 | 0.492 | 0.556 | 0.444 | 0.508 
500,000 | 10386 | 10382 | 10403 | 10365 | 0.555 | 0.439 | 0.561 | 0.445 
750 ,000 15027 | 15084 | 15073 | 15052 | 0.533 | 0.445 | 0.555 | 0.467 
0.502 | 0 0.418 | 0.498 


1 ,000 ,000 19617 19665 | 19621 | 19593 | . 582 





A table of 22,,(1008k) was computed for b = 1, 5, 7, and 11 andk = 1, 2, --- , 1000. 
(An interval of 1008 was chosen, since its divisibility by 12 simplified the computer 
program and eliminated even slight inequalities.) The results are similar: the 
deficiency of 4m + 1 resides solely in 12m + 1, and when A(n) < 0 we now find 
12m + 5 in the lead. Again, 


(16) fR=k=,=l, =0 


throughout the range computed. Again, the three stronger classes are all equally 
strong—see Table 5. 

In contrast now consider the four classes of primes modulo 10, 10m + 1, +3, 
+7, and +9. This time we find 


(17) fi:(1000k) = 1;(1000k) = 1,(1000k) = f,(1000k) = 0 


for k = 1, 2, --- , 1000 so that for these arguments neither 10m + 1 nor 10m + 9 
is ever in first place.* We find that 10m + 3 take turns in first place while 10m + 1 
take turns in last place. From Table 6 we conclude that this time we have two 
strong classes (which are equally strong), 10m + 3, and two weak classes (which 
are equally weak), 10m + 1. 

These striking distinctions between modulo 8 and 12 on the one hand and mod- 
ulo 10 on the other are easily explained in terms of the corresponding modulo 
multiplication groups, and this explanation provides us with a simple, number— 
theoretic, sufficient cause. 


7. Group Multiplication and Fluctuations. Consider the multiplication table 
modulo 8 of the four residue classes 1, 3, 5, and 7 


* At some intermediate arguments, e.g., n = 135852, 10m + 1 is tied for first place and simi- 
larly, for n = 969240, 10m + 7 is tied for last place. 
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and assume that in an interval centered around n there are fewer primes of the 
form 8m + 1 than of the other three types. Consider the excess primes of the other 
three classes. Their products with themselves and with each other are composites 
(of order n?) whose residue classes are contained in the lower right 3 x 3 box. Now 
note that these products are distributed into the four residue classes in proportions 
3:2:2:2—there is an extra product congruent to 1. Assume on the contrary an 
excess number of primes of the form 8m + 1. The composites from these primes lie 
in the upper left 1 x 1 box. Again we find an excess of composites congruent to 1. 
Similarly, too few or too many primes congruent to 3, 5, or 7 will also lead to an 
excess number of composites congruent to 1. Briefly, any fluctuation from equality 
in the distribution of the primes into the four residue classes will create an excess number 
of composites and therefore a diminished number of primes congruent to 1 modulo 8. 
This idea, that the phenomenon in question is essentially a fluctuation phenomenon 
(with analogies to the fluctuation phenomena of physics), arose during a discussion 
of this problem between the author and T. S. Walton. The idea is thus partly due 
to him. 

Since the multiplication group modulo 12 is isomorphic to that modulo 8, the 
behaviour found modulo 12 is similar. But for 10 we do not have the “vier” group 
but the cyclic group 

















3 9 1 7 

7 1 | 9 | 3 | 
| 
| 9 7 | 3 | 1 
| 








and this time we find that fluctuations in 10m + 3 or 10m + 7 diminish the 10m + 9 
primes while fluctuations in 10m + 1 or 10m + 9 diminish the 10m + 1 primes. 
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8. Quadratic Residues and the Distribution of Primes. In general, it is readily 
seen, those residue classes which occur on the principal diagonal are the ones whose 
numbers of primes are diminished. These are the quadratic residues. We expect then, 
for any modulus, that the prime race will separate (in the mean) into two races— 
that between the non-residues, up front, and that between the residues, in the rear. 
In particular, since +1 is the only quadratic residue of 24, the deficiency of 4m + 1, 
which was found to reside solely in 8m + 1 or 12m + 1, can be further localized 
to their intersection, 24m + 1. 

For most larger moduli, unlike the unusual 24, there will be many quadratic 
residues and the deficiency, ~/n/log n, will be shared among all of them. Thus the 
separation between the residues and the non-residues will not be as sharp, and 
extensive interplay is to be expected. 


9. Propagation of the Deficiency and Higher Order Effects. The quadratic 
fluctuation effect just discussed is not merely strong enough to maintain the mean 
deficiency in the 4m + 1 primes, but, as we shall see presently, as n increases it 
becomes too strong and we must examine the compensating cubic and higher order 
effects in order to obtain an accurate picture. We confine ourselves to the modulus 
4, although the generalization offers little difficulty. We have seen, (9), that the 
mean value of + up to n = 2-10° is nearly one and we wish to show how a mean 
value equal to one can propagate itself to larger values of n. Our computation, 
however, is only approximate, and while the result gives further evidence for the 
truth of (13), it is not a proof of that conjecture. 

Let a be a positive integer and let +.(n) be the number of positive integers 
of the form 4m + 1 which are <n and which are the products of a (not necessarily 


distinct) primes. Our previous 7,(n) is now 7{?(n). With a similar definition for 


a ( 


(n), let 
(18) A®(n) = x(n) — w{(n), 
(19) = (n) = w(n) + 2{(n), 
and 
(20) en) = A) Vp. 


>@(n) 


Note the slight difference between + and r“’—the new denominator, =“ (n) = 
a (n) — 1, omits the count of the prime 2. 

We wish to compute the mean value of r°’(n) and we assume that the mean 
values of +(x) for x = n“* (a = 2, 3, 4, ---) are all equal to 1. It follows that 
of all the odd primes of order n““, a fraction equal to (1 + n~’/**) is of the form 
4m — 1 anda fraction equal to}(1 — n~*/**) is of theform 4m + 1. The composites 
of order n, which are the products of a primes, have as prime factors primes whose 
(geometric) mean order is n'*. With reference to the multiplication table (mod 4): 


(1) 


we oe 
+1) —-1| 


—1 | +1 
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we now find, by induction on a, that among all possible a—fold products of these 
primes, a fraction equal to (1 + (—1)**'n’*) is of the form 4m — 1 and a fraction 
equal to $(1 — (-—1 )***n) is of the form 4m + 1. In other words, on the average, 


we will have, for a = 2,3,---, 


(21) #(n) = (-1)""". 


Thus we find a simple and interesting generalization of the Chebyshev phenomenon. 
There will be, in the mean, an excess of 4m — 1 composites of odd degree a, and 
an excess of 4m + 1 composites of even degree. Further, the mean fractional excess 
will be n™, independent of the degree a. 
Now, for any n, 
1 2) 2 (3)/ 
0 = A'’'(n) + A” (n) + A? (n) + --- 
with an error of 0 or 1, and therefore, with the same error, 
> (1) 2) 3 
(22) AY’ (n) = —A”(n) — A™(n) — --- 
As long as n is not too large, say 1000, most of the composites are of second degree, 
(1) Ms. wes : , : a 

and we have A*’(n) + —A”’(n). This, essentially, was our picture in sec. 7 above. 
But as n increases, the other even degree composites would tend to increase the 
deficiency while the odd degree composites tend to diminish it. To determine the 
balance we use the known fact, [5, p. 627], that in any arithmetic progression the 
number of integers with a odd is asymptotically equal to those with a even. Specifi- 


cally, for the progression 2m + 1, we have 


(1 + e) 2S” (n) = Z™ (wn) — ZS (wn) + SB (en) — --- 


, 
where « > 0 asn — ~. Therefore from (20), (21), and (22) we have, in the mean, 


A? (n) = (1 + &) B® (n)/Vn. 
Thus 
(23) t(n) =1t+. 


and a mean value of one propagates itself. 
There are two remarks of interest concerning equations (21) and (22). 
1.) The sequences x(n) and A(n) are defined as follows: 


A(n) = (—1)° 
where a is the number of prime factors of n, and x(n) = 0, 1, 0, and —1 for num- 
bers of the form 4m, 4m + 1, 4m + 2, and 4m + 3 respectively. Let 
(24) b, = A(n)x(n). 
A generalized Chebyshev assertion now reads: there are ‘‘many more’’ integers 


m for which b,, = +1 than for which b,, = —1 and the mean excess is of order 
44/n. Since it is known, [5, p. 674], that 


f(s) = Oban = (1 — 2)¢(28)/L(s), 
n=l 


it is clear that a function-theoretic formulation of this assertion is concerned with 
f(s) and with the zeros of L(s). 
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2.) Let n be sufficiently large such that the number of primes is small com- 
pared with the number of composites. Then the terms of like sign in the right 
side of (22) could be combined so that the right side would be the difference be- 
tween two nearly equal numbers. We therefore have an analogue of a difficulty 
in numerical analysis—subtractive loss of significance. In fact, the analogy is 
very good, since it is known, [15, p. 342], that the normal order of the number of 
prime divisors of n is log log n and that the number of composites of degree a is 
given asymptotically by 


x(n) ~ n(log log n)**/(a — 1)! log n. 


This means that the numbers of composites of degree a have an approximate Poisson 
distribution around a mean of log log n and that the equation (22) is comparable 
to the following equation, which is numerically sensitive when n is large: 
—1 l 2 
1—e “" = log log n — 4(loglogn)° + ---. 

We should therefore expect disturbances in the balance, (23), and more or less 
random oscillations around this mean. The oscillations should grow as log log n 
increases and may be expected to show strength when the condition =“ « 

2 3) : : oe 
= ~ >” is met. This should occur for log log n = 2.5+ or log log log n ~ 1. 

With these remarks we are led to our final topic. 


10. The Littlewood Phenomena. Analogous to Chebyshev’s assertion is the 
erroneous inequality: 
(24) a(n) < li(n) 


which was thought to be correct both by Gauss and by_Riemann, [7, p. 791, 795]. 
Again we have the Littlewood counterexample 





(25) a(n) — Li(n) > K vn log log log n. 
log n 


And again we obtain, this time from Riemann’s prime formula, [7, p. 795], 
(26) Li(n) = x(n) + 3x(/n) + ---, 

and from #(4/n) ~ 2 +/n/log n, the suggestion to define 

Li(n) — x(n) 


(27) rank o/n/log n 

Or similarly let us define 

’ _ Li(n) — x(n) - 
(28) p(n) = “a Vn. 


Analogy with (13) now suggests the possibility that 


(29) Lim : 


Now fi i 1 





z 
» p(n) = +1 


is the proper formulation of the erroneous (24). 
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If, indeed, the differences Li(n) — x(n) and x(n) — x,(n) are usually of 
the same sign and order of magnitude we should find that 


(30) a*(n) = x(n) + a(n) — w4(n) = 2x_(n) + 1 


agrees with Li(n) much better than x(n) does. This is, in fact, the case as is seen 
in Table 7, [16]. We therefore regard (29) as plausible. It should be investigated. 

A final remark concerning the Littlewood phenomena may be of value. It is 
sometimes said or implied, [17], that (24) is true ‘‘im Bereich der Tabellen’’, where 
the tables in question include such entries as n = 10° and n = 10°. But x(n) — Li(n) 
has not been evaluated for all intermediate values of n (say for 10° < n < 10°) 
and consequently, it is not certain that (24) holds there. In the analogous case 
of A(n) discussed above, had we computed A(10,000k) up to k = 300 we would 
have always found A > 0 since all the axis crossings of A(n) up to n = 3.10° occur 
at intermediate values of n. Further, since the first crossing of x(n) — Li(n) may, 
like Leech’s 26861 above, be of very short duration, it could well be missed unless 
the computing interval in n were rather small. The author knows of no compelling 








TABLE 7 

” x(n) x*(n) a(n)/Li(m) x*(n)/Li(m) 

10° 168 175 0.9459 0.9853 
2-108 303 311 0.9625 0.9879 
3-108 430 437 0.9712 0.9870 
4-108 550 561 0.9728 0.9923 
5-108 669 679 | 0.9777 0.9923 
6- 10° 783 799 | 0.9782 0.9982 
7-108 900 | 915 | 0.9843 1.0007 
8-10 | 1007 | 1015 0.9811 0.9889 
9-10° | 1117 i125 | 0.9825 0.9895 

10° | 1229 1239 | 0.9862 0.9943 
2-10 | 2262 2273 0.9884 0.9932 
3-10° | 3245 3267 0.9903 0.9970 
4-10¢ | 4203 4235 0.9929 1.0005 
5-10* | 5133 5167 0.9935 | 1.0001 
6-104 6057 6077 0.9958 | 0.9990 
7-108 6935 6971 0.9928 | 0.9980 
8-10 | 7837 | 7867 0.9950 0.9988 
9-10° | 713 8729 0.9949 0.9968 

10° | 9592 9617 0.9961 0.9987 
2-108 | 17984 18013 0.9971 0.9986 
3-108 | 25997 26033 | 0.9965 0.9979 
4-10° | 33860 33919 | 0.9981 0.9999 
5-108 | 41538 41613 | 0.9983 1.0001 
6-105 49098 49151 0.9985 0.9995 
7-105 56543 56589 | 0.9982 0.9990 
8-10° 63951 64071 | 0.9986 1.0005 
9-108 71274 | 71379 | 0.9988 1.0002 

10° 78498 | 78645 0.9983 1.0002 
2 0 


108 148933 149033 0.9992 


.9998 














284 DANIEL SHANKS 


reason why the first crossing should occur at anywhere near as large an n as Skewes’ 
. 1034 
fantastic 10” , [18]. 
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Tables of Abscissas and Weights for Numerical 
Evaluation of Integrals of the Form I e*x"f(x) dx 


By Philip Rabinowitz and George Weiss 


There are many physical problems [1, 2] whose solutions involve integrals of 
the form 


(1) [ven dz, 


hence it is desirable to have a table of the weights and abscissas suitable for use 
in the numerical integration formula of Gaussian type 


(2) [ xe “f(x) dx = > din f(tiw) + E; 
=) 
where 
NIN ae 
E= N\(N +n)! rie . 


(2N)! 


It is known that if Ly(x) is the generalized Laguerre polynomial defined by 


N r 
=o... = (N+ a\(—2)" 
(3) Ly(x) = yy S i ") a 
then [3] the abscissas xj appearing in Eq. (2) are the roots of Ly(x) = 0, while 
the coefficients ajy are given by 


1 1 

Lkw [Ly (xiw)}” 
where L’ denotes the derivative of L. 

The most extensive tables of aty and xy have been given by Salzer and Zucker 
[4] who tabulated the weights and roots to at most 14 places for N = 2(1)15. A 
very short table of the ajy and zy to four places has been given by Burnett [5] 
for n = 2, 3,4 and N = 2, 3. The present tables contain the roots, weights, and 
the weights multiplied by exp(ziy) for use in the numerical integration formula, 
Ey. (2). The tables are good to 18 significant figures for n = 0(1)5. For n = 0 
we give the appropriate numbers for N = 4(4)32 and for n > 0 we give the num- 
bers for N = 4(4)16. 

As an illustration of a possible use of these tables we cite the calculation of 
such lattice sums as [6] 


(4) din = 


“ eo rT 2 
n'm 

(5 T +224 = ——— , 

oar 2 mai (m* + n?)! 
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By using the identity 


—z t-1 
e “udu 


2 
0 


z= {r()}7 | 


we may convert the expressions for 7... into an integral form 


(6) Tint = =o u'a(u)o,(u) du 
where 
(7) o(u) = >one™™ 


n=l 


which has an exponential behavior. Previous tabulations of lattice sums have relied 
on the use of the Poisson summation formula to find the behavior [6] of o,(«) for 
small u or on direct summation of a large number of terms of Eq. (5) by electronic 
computer [7]. The present tables enable one to evaluate the integral in Eq. (6) by 
direct numerical integration. To facilitate this method we are now computing 
tables of theta functions evaluated at the zeros of Laguerre polynomials. 

The values in the tables were computed on WEIZAC using double precision 
floating operations which are accurate to approximately 21 significant figures. The 
zeros of Ly (x) were computed using Newton’s method. If xis" is the ith approxi- 
mation to the kth root of Ly(z), then 


i+1 ) i) 7 dy (i 
(8) tin = apn — Ly(akn’)/Ly (atx) = tin” — A 
where Ly(x) is calculated from the recurrence relation [3] 


NLIx(x) — (2N —1+ 7 — x)Ly-(x) + (N — 14+ n)Ly_2(x) = 0, 


(9) 
iA ae: | Be) = 1, 
and 
(10) Lt’ (2) = ~ (NL¥(2) — (N + n)L3_4(2)). 


The iteration was terminated when | A| < 2™”. For the initial guesses for the 
iteration, various sources were used. For n = 0, N = 4, 8, 12 the values given by 
Salzer and Zucker [4] were used. For n = 0, N = 16, 20, the values given by Head 
and Wilson [8] were used. In all other cases, we made use of the following inequality 
given by Szego [9] 


ae 
4N + 2(n + 1) 


_ (Qk + n+1) + (k+ n+ 12k +n4+ 1) +3 —- ni} 
i 2N+n+1 


< Xin <~ R 





(11) 








Ww 
fi 
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where jx, is the kth zero of J,(x) and is tabulated in Watson [10] for n = 0(1)5, 
k = 1(1)40. The way this was done is as follows. First we determined 
iw = (tin — A)/(B — A) forn = 0, N = 15, 16 using the values given by Head 
and Wilson [8]. From the two graphs of Aj.s and Ais we estimated that 
4(k/(2N — 1))°* was a good approximation to zy , independent of n. Hence our 
initial guess for xiy was 

2 k ¥ 2 2 

Jin + 8 (sy*5) (2k +n+1)° — 47%, + (2k +n+ 1) 

3) _ ‘{(k+n+1)?+4—n}'] 
i = 7 - . 
4N + 2n+2 


(12) 





This choice of zis” was a rather good approximation to xfy , differing from it 
by at most 10% and usually by less than 1%, and always being closer to it than 
to any other root of Ly(x). Convergence to zjy occurred in almost every case 
in at most 10 iterations, and usually in much less. However, in some cases in which 
k = N, Newton’s method favored convergence to ry_:,~ . This caused very little 
trouble for an approximate value for xy,y was computed from the relation derived 
from Eq. (3) 

N 


(13) > thy = N(N +n) 


k=1 


and used as an initial value in the iteration. In general, this relation together with 
a second relation derived from Eq. (3) 


(14) I] atv = (N +n)! 


' 


k=l nN: 


were used as checks on the computation. They always checked to 18 significant 
figures. In addition the following quantities were computed by the machine and 
these also checked to 18 significant figures (all summations for k = 1, --- , N). 


> atv = n! 
Daivtiy = (n + 1)! 
> aty(xtw)? = (n+ 7)! 


Diaiw(ziw)” = (nm + 15)! N ¥4 
Dcaiw(xiw)” = (n + 23)! N ¥4,8 
{ Deatw( xen)” = (n+ 31)! n= 0,1,2;N = 16 
| Satw(xkv)*? = (n + 27)! n = 3,4,5;N = 16. 
The values for n = 0, N =’32 were not checked by the machine and were only 


checked by hand using Eq. (13), while for n = 0, N = 28, only Eqs. (13) and 
(14) were checked by the computer. This modification was necessary because the 
values a%s.os , @31.32 , ANd ay2,32 were out of the range of the double precision floating 
routine and had to be computed by hand. 








TABLE 








Zeros of Generalized Laguerre 
Polynomials 


Weights 


Weight times Exponential of 


=a" -) 
Zero aun oP (iy 





3 .22547689619392312 (—1) 
1.74576110115834658 (0) 
4 .53662029692112798 (0) 
9.39507091230113313 (0) 


.70279632305101000 (— 1) 
-03701776799379912 (—1) 
-25108662986613069 (0) 
-26670017028765879 (0) 
-04590540239346570 (0) 
-07585160101809952 (1 ) 
-574067864 12780046 (1 ) 
-28631317368892641 (1) 


Nee AN Oe 


.15722117358020675 (—1) 
.11757484515130665 (—1) 
.51261026977641879 (0) 
.83375133774350723 (0 ) 
.59922763941834848 (0) 
.84452545311517735 (0) 
-62131684245686704 (0) 
.30060549933063477 (1) 
-71168551874622557 (1) 
.21510903793970057 (1 ) 
.84879672509840003 (1) 
.70991210444669203 (1) 


WONNK RK OCOrPNR 


.76494104789278403 (— 2) 
.626963289 15080832 (—1) 
.14105777483122686 (0) 
. 12928364509838062 (0 ) 
.43708663389320665 (0 ) 
.07801861454976791 (0) 
.070338535048234 13 (0) 
.43831433639193878 (0) 
.22142233688661587 (1) 
.54415273687816171 (1) 
-91801568567531349 (1) 
.35159056939919685 (1 ) 
.85787297428821404 (1) 
-45833987022866258 (1 ) 
. 19404526476883326 (1) 
.17011603395433184 (1) 


OrPRWNN RR eH OAC WN eS OO 


7 .05398896919887534 (—2) 
3.72126818001611444 (—1) 
9. 16582102483273565 (— 1) 
1.70730653102834388 (0) 
2.74919925530943213 (0) 
4 .04892531385088692 (0) 
5.61517497086161651 (0) 
7 .45901745367 106331 (0) 
9 .59439286958109677 (0) 
1. 20388025469643163 (1) 
1 .48142934426307400 (1 ) 
1.79488955205193760 (1 ) 
2.14787882402850110 (1) 
2.54517027931869055 (1) 
2.99325546317006120 (1) 
3.50134342404790000 (1 ) 
4.083305705672857 11 (1) 
4.76199940473465021 (1) 
5 .58107957500638989 (1 ) 
6 .65244165256157538 (1) 


n=0,N =4 
.03154104341633602 (— 1) 
.57418692437799687 (—1) 
.88879085150053843 (— 2) 
.39294705561327450 (—4) 
n=0,N=8 
.69188589341637530 (— 1) 
.18786780814342956 (— 1) 
. 75794986637 171806 (— 1) 
.33434922612156515 (—2) 
.79453623522567252 (—3) 
9 .07650877335821310 (—5) 
8 .4857467 1627253154 (—7) 
1.0480011748715' °8(—9) 
n=0,N = 12 
2.64731371055443190 (— 1) 
3.77759275873137982 (—1) 
44082011319877564 (— 1) 
044922221 16809307 (—2) 
.01023811546340965 (— 2) 
66397354186531588 (— 3) 
.03231592662999392 (— 4) 
36505585681979875 (— 6) 
66849387654091026 (—7 ) 
34239103051500415 (—9) 
.06160163503502078 (— 12) 
14807746742624168 (— 16) 
n=0,N = 16 
.061517 14957800994 (— 1) 
31057854950884166 (— 1) 
.65795777644214153 (—1) 
36296934296377540 (—1) 
.73289286941252190 (—2) 
12999000803394532 (—2) 
.84907094352631086 (—3) 
04271915308278460 (—4) 
48445868739812988 (—5) 
. 82831933087 119956 (—7) 
88102484107967321 (—8) 
86235024297388162 (— 10) 
12707903322410297 (— 12) 
).29796700251786779 (— 15) 
.05047370003551282 (— 18) 
.16146237037285519 (— 22) 
n=0,N = 20 
68746801851113862 (— 1) 
.91254362006068282 (— 1) 
66686102867001289 (— 1) 
.66002453269506840 (— 1) 
.48260646687923705 (—2) 
.49644173092832211 (—2) 
.20255084457223685 (—3) 
.14496238647690824 (—3) 
.55741773027811975 (—4) 
.54014408652249157 (—5) 
.08648636651798235 (—6) 
. 3301209095567 1475 (— 8) 
.75798117905058200 (—9 ) 
.72550240251232087 (— i1) 
.76752925157819052 (— 13) 
.37284424336243841 (— 15) 
.15501433950039883 (— 17) 
.53952214058234355 (— 20) 
. 286442725569 15783 (— 24) 
.65645661249902330 (— 28 ) 


Owwa 


CO We WDNNNM Oh Nowe OO 


NWMWRKORr NR KH FR to wd 


OI 


Se Nl oll el el ee eee 


Or PWNNN HR Ooh 


8. 
2. 
3. 
6. 


4. 
1. 
1. 
2. 
3. 
4. 
5. 
8. 


2. 
6. 


eNO OWN Ree 


= 00 > 


MONO ROWN DNDN RR HO er 


32739123837889247 (—1) 
04810243845429682 (0) 
63114630582151786 (0) 
48714508440766227 (0) 


37723410492911373 (—1) 
033869347665597 64 (0) 
66970976565877575 (0) 
37692470175859948 (0) 
2085409 1334792628 (0) 
26857551082513220 (0) 
81808336867 192193 (0) 
90622621529221140 (0) 


97209636044410798 (—1) 
96462980430597231 (—1) 
10778139461575215 (0) 
53846423904282912 (0) 
99832760627424276 (0) 
50074576910086687 (0) 


.06532151828238624 (0) 


7232891 1078277160 (0) 


.52981402998173506 (0) 


59725846183532109 (0) 
21299546092587700 (0) 


.05438374619100811 (1) 


. 250363 14864247252 (—1) 
. 25836052762342454 (—1) 


31961391687087088 (— 1) 
14609924096375170 (0) 
47175131696680859 (0 ) 
81313468738134816 (0) 
17551751969460745 (0) 
56576275016502921 (0) 


99321508637 137516 (0) 


47123448310209029 (0) 


.02004408644466886 (0) 


67251660773285426 (0) 
48742065798615247 (0) 


).58536123328921366 (0) 
.27635798436423448 (0) 


18242775516584348 (1) 


$1080062418989255 (— 1) 
22556767878563975 (—1) 
66909546701848151 (—1) 
15352372783073673 (— 1) 
169539707 19554597 (0) 
43135498592820599 (0) 
702981 13798502272 (0 ) 
98701589079274721 (0) 
28663578125343079 (0) 


.60583472755383333 (0) 


94978373421395086 (0) 
3253957820093 1955 (0) 
7422554705898 1092 (0) 
2142367 1025188042 (0) 


. 76251846 149020929 (0) 
-42172604424557430 (0) 
.25401235693242129 (0) 
.38731438905443455 (0) 
. 15132873098747960 (0) 
.28933886459399966 (1 ) 
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Zeros of Generalized Laguerre 
olynomials 


Weights 


Weight times Exponential! of 
Zero = “om exp iy 





5.90198521815079770 (—2) 
3 .11239146198483727 (—1) 
7 .66096905545936646 (— 1) 
1 .42559759080361309 (0) 
2 .29256205863219029 (0) 
3 .37077426420899772 (0) 
4 665083703467 17079 (0) 
6. 18153511873676541 (0) 
7 .92753924717215218 (0) 
9 .91209801507770602 (0) 
1 .21461027117297656 (1) 
1 .46427322895966743 (1) 
1 .74179926465089787 (1) 
2 .04914600826164247 (1) 
2 .38873298481697332 (1 ) 
2 .76359371743327174 (1) 
3 .17760413523747233 (1) 
3 .63584058016516217 (1) 
4 .14517204848707670 (1) 
4.71531064451563230 (1) 
5 .36085745446950698 (1) 
6. 10585314472187616 (1) 
6 .99622400351050304 (1 ) 
8. 14982792339488854 (1 ) 


5 .07346248498738876 (—2) 
2.67487268640741084 (—1) 
.58136628354791519 (—1) 
. 22397 180838490772 (0) 
. 9667676 1247377770 (0) 
.88888332603032189 (0) 
.99331165925011414 (0) 
.28373606284344256 (0) 
.76460340424350515 (0) 
.44121632827 132449 (0) 
.03198504629932601 (1) 
.240790341446067 17 (1) 
.47140851641357488 (1) 
. 72486634 156080563 (1 ) 
.00237833299517127 (1) 
.30538901350302960 (1 ) 
.63562973744013176 (1) 
-99519668335961821 (1) 
38666055 165844592 (1 ) 
.81322544101946468 (1 ) 
.27896723707725763 (1) 
. 789207 16336227 437 (1) 
.35112979596642942 (1) 
.974879608464 12408 (1) 
.67569772839064696 (1 ) 
.47867781523391618 (1) 
.4317837 1072270431 (1) 
.65824206275273191 (1) 


o>) 


CWONOMINS SWWNNNN RH RH HOO OWN Ree 


.44893658332670184 (—2) 
.34526109519618537 (—1) 
.76884629301886426 (—1) 
.07244875381781763 (0) 
.72240877644464544 (0) 
.52833670642579488 (0) 
.49221327302199449 (0) 
.61645676974976739 (0) 
.90395850417424395 (0) 
.35812673318624111 (0) 
.98294092421259610 (0) 
.07830186325399721 (1) 
.27636979867427251 (1) 


CONT OT ROOD RON 


n=0,N = 24 
1 .42811973334781851 (—1) 
2.5877410751 7423903 (—1) 
2 .58806707272869802 (—1) 
1.83322688977778025 (—1) 
9 .81662726299188922 (—2) 
4.07324781514086460 (— 2) 
1.32260194051201567 (—2) 
3 .36934905847 830355 (—3) 
6 .72162564093547890 (—4) 
1 .04461214659275180 (—4) 
1 .25447219779933332 (—5) 
1. 15131581273727992 (—6) 
7 .96081295913363026 (—8) 
4 .07285898754999966 (—9) 
1 .50700822629258492 (— 10) 
3 .91773651505845138 (— 12) 
6 .89418105295808569 (— 14) 
7 .81980038245944847 (— 16) 
5 .35018881301003760 (— 18) 
2 .01051746455550347 (— 20) 
3 .60576586455295904 (— 23 ) 
2.45181884587840269 (— 26 ) 
4 .08830159368065782 (— 30) 
5 .57534578832835675 (— 34) 

n=0,N = 28 
1 .23778843954286428 (—1) 
2.32279276900901161 (—1) 
2 .47511896036477212 (—1) 
1 .92307113132382827 (—1) 
1. 16405361721 130006 (—1) 
5 .63459053644773065 (—2) 
2 .20663643262588079 (—2) 
7 .02588763558386773 (—3) 
1 .82060789269585487 (—3) 
3 .83344303857123177 (—4) 
6 .53508708069439831 (—5) 
8 .97136205341076834 (—6) 
9 .84701225624928887 (—7) 
8 .56407585267304245 (—8) 
5 .836838763 13834429 (—9) 
3 .07563887784230228 (— 10) 
1 .23259095272442282 (—11) 
3.68217367410831200 (— 13) 
7 .99879057596890965 (— 15) 
1 .22492250032408341 (— 16) 
1.27112429503067374 (— 18) 
8 .48859336768654320 (—21 ) 
3 .40245537942551185 (— 23) 
7 .420156588867485 13 (— 26 ) 
7 .60041320580173769 (— 29) 
2.87391031794039581 (—32) 
2.5418229038893 1800 (— 36) 
1 .66137587802903396 (— 41) 

n=0,N = 32 
1.0921834195238497 1 (—1) 
2.10443107938813234 (—1) 
2.35213229669848005 (—1) 
1.95903335972881043 (—1) 
1. 2998378628607 1761 (—1) 
7 .05786238657174415 (—2) 
3.17609125091750703 (—2) 
1.19182148348385571 (—2) 
3.73881629461 152479 (—3) 
9 .80803306614955132 (—4) 
2.14864918801364188 (—4) 
3.92034196798794720 (—5) 
5.93454161286863288 (—6) 


1.51494412859509452 (—1) 
3 .53256582529923847 (—1) 
5 .56784563288152601 (—1) 
7 .62685317697309109 (—1) 
9 .71872632246547580 (— 1) 
1 . 18535789303780109 (0) 
1 . 4042656272844 1853 (0) 
1 .6298686 1575704148 (0) 
1 .86363505533207295 (0) 
2. 10729115108148019 (0) 
2 .36290589 104193500 (0) 
2 .633008753 16385675 (0) 
2 .92075757972772467 (0) 
3 .23018513349235362 (0) 
3 5665733773687 5657 (0) 
3 .93704375545516030 (0 ) 
4 .35153118886351197 (0) 
4 .82448185489803577 (0) 
5.37802207978918232 (0) 
6 .04841781261995523 (0) 
6 .90089835218049651 (0) 
8.069965 156 14695609 (0 ) 
9 .90279331948422546 (0) 
1 .38205320947920057 (1 ) 


. 302207 49260603611 (—1) 
.03513987299129915 (— 1) 
77992610899299887 (—1) 
.53972597062647 847 (—1) 
.3201 1455708276465 (—1) 
.01271781273780542 (0) 
.19675156048109225 (0) 
.38483686392073 190 (0 ) 
.577779655787 12887 (0) 
.7764896694435896 1 (0 ) 
.98200839256 100634 (0 ) 
19554512993098825 (0 ) 
41852439028379250 (0 ) 
.65264936115550918 (0 ) 
89998868599490885 (0 ) 
. 16309775591259801 (0 ) 
44519246556 198147 (0) 
750405 15240709831 (0 ) 
.08417387436021254 (0) 
.453857 16433002825 (0 ) 
4 .86974939776545707 (0) 
5 .34685205218439829 (0) 
5 .90818195463222266 (0) 
6 .591519030727657 14 (0) 
7 .46490124681590107 (0) 
8 .66878601010575288 (0 ) 
.0566 1475342792109 (1) 
.464469151725447 14 (1) 


BP ePWWWNNNN Hee HHH OOS We 


— ht 


.14187105768104849 (—1) 
.66065216897615217 (—1) 
18793137324852994 (—1) 
.72532846499804707 (—1) 
.27648788380971311 (—1) 
. 845367 19340249717 (—1) 
.04361887589207697 (0) 
. 20534927415235258 (0) 
.37022133852178119 (0) 
.53877725646864475 (0) 
.71161935268645726 (0 ) 
.889424063449484 10 (0) 
2 .07295934024653367 (0) 


et et et et OO I OT oe 
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Zeros of Generalized Laguerre 
‘olynomia! 


Weights 


Weight times Exponential of 


Zero = a” ) 


fon 
kN exp (x 


kN 





1.49311397555225573 (1) 
1.72924543367153148 (1) 
1.98558609403360547 (1) 
2.26308890131967745 (1) 
2 .56286360224592478 (1) 
2 .88621018163234747 (1) 
3 .23466291539647370 (1) 
3.61004948057519738 (1) 
4.01457197715394415 (1) 
4 .45092079957549380 (1) 
4 .92243949873086392 (1) 
5.43337213333969073 (1) 
5.98925091621340182 (1) 
6.59753772879350528 (1) 
7 .26876280906627086 (1 ) 
8.01874469779135231 (1) 
8.87353404178923987 (1) 
9 .88295428682839726 (1 ) 
1.11751398097937695 (2 ) 


7 .43291927981431435 (— 1) 
2.57163500764627847 (0) 
5.73117875168909963 (0) 
1 .09538943126831905 (1) 


4 .09383573203185153 (— 1) 
1.38496318480313988 (0) 
2 .95625455616886207 (0) 
5.18194310104007138 (0) 
8.16170968814581733 (0) 
1 .20700551268371548 (1) 
1.72497355261489875 (1) 
2.45859552436527819 (1) 


2.82858348239914112 (—1) 
.52326041364613503 (— 1) 
.01649213857768073 (0) 
-49235406977802504 (0) 
.40549102001572284 (0) 
.79281394041241498 (0) 
1 .07073886889890908 (1 ) 
1 .42271523637899807 (1) 
1.84719966342299991 (1) 
2 .36417837524008953 (1 ) 
3 
3 


“IC1Ww nO 


.01200586261063443 (1) 
.88892843760953185 (1) 


2.16140305239452255 (— 1) 
7 .26388243251803954 (—1) 
1.53359316037354132 (0) 
2.64497099861191096 (0) 
4 .07097816088019065 (0) 
5.82585551510560455 (0) 
7 .92850418530666709 (0) 
1.04038082899510393 (1) 
1 .32846610707070382 (1 ) 
1.66151732168666126 (1 ) 
2 .04560060200272169 (1) 
2 .48938470253519108 (1) 
3 .00598629202025761 (1 ) 
3 .61706945436791780 (1) 
4 .36403651841768370 (1) 
5.35291511602684204 (1 ) 


41640457866755222 (—7) 
60456787912078148 (— 8) 
35060222662580674 (—9) 
. 28138297 104092888 (— 10) 
30589949189133608 (— 11) 
79937928872709406 (— 13) 
23780165772926646 (— 14) 
17182344342071943 (— 16) 
54213383339382337 (— 17) 
11979229016361861 (—19) 
.05442967378804543 (— 21) 
34698258663739516 (— 23) 
.66129413039735937 (— 26) 
41856054546303691 (— 28) 
91337549445422431 (—31) 
19224876009822236 (— 34) 
67151121924013699 (—38) 
33861694210625628 (— 42) 
.51053619389897424 (— 48) 
n=1,N=4 
46870593218776310 (—1) 
77635772363868313 (—1) 
41777847310521364 (—2) 
31584968630324014 (—3) 
n=1,N=8 
7632541405723668 (— 1) 
.38985360731142368 (— 1) 
.89996070781313351 (—1) 
.51413846166973434 (—2) 
.93264664870734492 (—3) 
.08642136813304225 (— 4) 
34895820979709614 (—6) 
.72139282319322132 (—9) 
n=1,N = 12 
.01453765665592267 (—1) 
17699483362918101 (—1) 
3998627 1374452795 (—1) 
79532599439199819 (—1) 
20319390353283915 (—2) 
49161743576890271 (—3) 
67135797278597991 (—4) 
63620912038435440 (—5) 
18456700779946096 (—7) 
.32315410040323253 (—9 ) 
.83970243442130593 (— 11) 
.37762852227037045 (— 15) 
n=1,N = 16 
}.32773328795394308 (— 2) 
31090461520719732 (—1) 
.16933542163999097 (—1) 
.37894217875244607 (—1) 
10272743359247882 (— 1) 
30895883563969837 (— 2) 
.52571640779246630 (—3) 
42711436792047207 (—4) 
).99619095556818742 (—5) 
.61208604213355474 (—6) 
.10154125147424469 (—7) 
.83568251273250421 (—9) 
.48148886588507814 (— 11) 
.73484535903961320 (— 14) 
4 .08414605500225599 (— 17) 
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3.62395496247891406 (—21) 
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2.26310663399696361 (0) 
2 .46088907248823613 (0) 
66750812639711715 (0) 
88439209292204179 (0) 
11326132703958617 (0) 
35621769259580256 (0) 
61586985648426880 (0 ) 
89551304494854956 (0) 
19939410471158549 (0) 
.53311497853436176 (0) 
4.90427028761124484 (0) 
5.32350097202366611 (0) 
5.80633321423362136 (0) 
6.37661467415965254 (0) 
7 .07352658070724211 (0) 
.96769350929590065 (0) 
-20504033127818968 (0) 
.11630130907678735 (1) 
.53901804152606398 (1 ) 


bo 


— me Os] 


39700286292286231 (—1) 
25091709508283381 (0) 
28714219375495781 (1) 
523538307768847 15 (1 ) 


NINO Oo 


82553823014304850 (—1) 
75360553036282958 (0 ) 
.57541442577641522 (0) 
33773006498366668 (1 ) 
77973624568047784 (1 ) 
38782591536632858 (1 ) 
.03842391519935203 (2) 
.24706083439768244 (2 ) 


NR orl Orr bo 


_ 


.34620745484392543 (—1) 
-8 .23391640164269559 (— 1) 
2.55395234400434764 (0) 
5.90001913843717294 (0) 
15836368199363577 (1) 
.05762635029790113 (1) 
42793100498839410 (1) 
48810758737070439 (1 ) 
.61552836575302900 (1 ) 
35579615855293272 (2) 
21678124758432626 (2 ) 
16862096414656053 (2) 


eR WOW nN 


85446678289795171 (—2) 
. 778025087 13723198 (—1) 
.46892374608307634 (0) 
35027298906274298 (0) 
46355621447558663 (0) 
12157486249617326 (1) 
.81106406690387663 (1) 
779683947 48499728 (1) 
11449606829674792 (1) 
.93794445636475115 (1) 
43201366918135643 (1) 
18866087257456150 (2 ) 
.68085784287149817 (2) 
.42124561031249231 (2) 
66328790805694462 (2) 
}.40615970837448162 (2) 


WN KOON RF OW =I 








TABLE—Continued 





Zeros of Generalized Laguerre 
Polynomials 


Weights 





1 .22676326350030207 (0) 
3.41250735869694597 (0) 
6 .90269260585161340 (0) 
1 .24580367719511386 (1) 


6 .99330392297772446 (—1) 
1.89881649533754637 (0) 
3.67761476834163509 (0) 
6 .09929454816086345 (0) 
9 .26742581328238619 (0) 
1 .33607382722601 106 (1) 
1.87281386688430816 (1) 
2.62686410414766043 (1) 


.90239109177454061 (—1) 
_32377645578306233 (0) 
_54213223523224643 (0) 
16451935324397066 (0) 
'21800163077006985 (0) 
74078952492664456 (0) 
17872181154281359 (1) 
.54365948600615973 (1) 
.98104815433429478 (1) 
511111947685387 « . (1) 
17262495017079404 (1) 
_06488781934720538 (1) 


PON RR RH OOH 


.77613508344740734 (—1) 
.01749195760256570 (0) 
.94775802042424383 (0) 
. 17692724488986868 (0 ) 
.71624006979179562 (0) 
.58058826577491250 (0) 
78946527064707413 (0) 
13683230828333189 (1) 
43506267274370444 (1 ) 
77810957248416460 (1) 
172108479657 13089 (1) 
625813867511 10680 (1 ) 
15245960042758187 (1) 
.77389210025289391 (1 ) 
53185461100898426 (1) 
53325835388358 122 (1 ) 


> ee OO CO 


_ 
S 


Cre OWN Nee OO 


75552164718549145 (0) 
26560586565682345 (0) 
.05794068313800185 (0) 
. 3920931 8040196832 (1 ) 


— OO ie 


.02996168735087822 (0) 
.43991423401425973 (0) 
.4131867638385 1689 (0 ) 
.01921044285466378 (0) 
.03653585615970610 (1 ) 
.4634328491 1153273 (1) 
.01808472740399451 (1) 
.79171925451893479 (1) 


we ee 


.31333453524153573 (—1) 
. 7220058777 1537790 (0) 
.08711820886747629 (0) 
.84907183643667852 (0) 
.03654192274923744 (0) 
9 .68900419917576858 (0) 
1. 28619916279440914 (1) 


ONT Com I 


moe NT 


doe bo bOI 


Oe me eR OO Pe OO 


WWW RAI W eR WOO 


—s10e 


n=2,N=4 
.25524997698654378 (— 1) 
.06342429197919458 (0) 
.06696131028353551 (— 1) 
35457929379748887 (—3) 
n=2,N=8 
.24797751904342970 (—1) 
-99530958902891337 (—1) 
13670609157191155 (—1) 
31541145990131662 (—1) 
.91357867305206514 (—2) 
30770014757761649 (—3) 
60219339506628624 (—5) 
52333939442792311 (—8) 
n=2,N = 12 
51482251459639793 (—2) 
77196653262217159 (—1) 
.20466902576669237 (— 1) 
.93826609735938395 (— 1) 
.75757334769862195 (—1) 
.38929300594179609 (—2) 
.51882418847753131 (—3) 
.87787474819631979 (—4) 
.68668021445018989 (— 6) 
.59806360065300972 (—8 ) 
.25743502885190380 (— 10) 
.99508663346059607 (— 14) 
n=2,N = 16 
.86064094670787025 (— 2) 
.93347390190440231 (—1) 
.83219363383550989 (— 1) 
.81874148596173479 (—1) 
.38180537473790217 (—1) 
.22105963944979891 (—1) 
.81146258006637302 (—2) 
. 14314919248225670 (—3) 
.85648533767438300 (— 4) 
. 2015800563 1090712 (—5) 
.34236243815651751 (—7) 
.32646044204804137 (—8) 
. 15266648290842947 (— 10) 
.94706915124608697 (— 13) 
.63797825636053360 (— 16) 
.45457612313612400 (— 20) 
n=3,N=4 
. 860334074 14649995 (0) 
. 35689 101902892032 (0) 
.64453972843517610 (—1) 
.83209339810621145 (—2) 
n=3,N=8 
47659268688580562 (— 1) 
12957515366677356 (0) 
37008143827449573 (0) 
13728082731325425 (—1) 
32197368598833373 (—1) 
.66810091481878879 (—3) 
.04305893787 150415 (—5) 
.56535793848369390 (—7 ) 
n = 3, N = 12 
.52730239391842707 (—1) 
.07163463296299004 (0) 
.09369597825 120887 (0) 
.75927719487292856 (0) 
.38876099420103765 (— 1) 
. 6353847 4605744857 (— 1) 
.90887802863997279 (—2) 


Weight times Exponential! of 
n nn 
kN 7-7 iy) 


Zero = a 


47416635265254532 (0) 
226559603 12239398 (1 ) 
-05652340494177747 (2) 
12047669645817929 (3) 


m bo Ww bo 


-52384060238173233 (—1) 
33925681844090349 (0 ) 
82270577764238765 (1 ) 
-03161617804640828 (2 ) 
.08474054427306454 (2 ) 
29856650024944985 (2 ) 
17892467611196491 (3) 
}.46102645632067 593 (3 ) 


dO We dO Ore 


~ 


55349132168384772 (—1) 
79310697 194985889 (0) 
15478262411010895 (0) 
-17835397356631661 (1) 
-81773930663432268 (1 ) 
11926747809213372 (2) 
.62936424650790635 (2 ) 
-4993653053 1789858 (2 ) 
88125704585776709 (3 ) 
69994780357075735 (3 ) 
.55128122489510488 (3 ) 
-79932534090427829 (4) 


KWH OF NOWOR 


-09069825088867 147 (—2) 
11471686778248555 (—1) 
0900865 1869529020 (0 ) 
39492556822190547 (1 ) 
77907618869590957 (1) 
803392147 18968970 (1 ) 
8456467 1246403063 (2 ) 
-58534172771350011 (2) 
}.58551410235728057 (2) 
16135832614353863 (3 ) 
-99151272949358401 (3 ) 
. 36092665 129919429 (3 ) 
65798450641837801 (3 ) 
68455967415686438 (3 ) 
747650568202679 12 (4 ) 
707067273202597 18 (4) 


Swe Owe O*! 


worounw-- 


.07647588421845741 (1) 
39037661575137080 (2 ) 
41474077829377203 (3 ) 
.03578296138366883 (4) 


NwNe 


.25387504291512878 (0) 
44306070979065246 (1) 
-956076917 17435238 (2) 
-02145993175804680 (3 ) 
. 19605176057880639 (3 ) 
51220362170062521 (4) 
-25709678492352907 (4) 
-083993368507 16076 (5 ) 


Noe eR te 


3. 17350448025970021 (— 1) 
5.99659048258816783 (0) 
4.58809098751416806 (1 ) 
2. 24522264593643932 (2) 
8 .40432710605229296 (2 ) 
2.63937440838205532 (3 ) 
7 .356481878887 16898 (3 ) 
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TABLE—Continued 








Zeros of Generalized Laguerre 
Polynomia! 


Weights 


Weight times Exponential of 


= n nn 
Zero = a”. exp ey) 


kN 





1.66361029842912300 (1) 
2.11344882654905403 (1) 
2 .65616764973530448 (1 ) 
3330947935537 48032 (1) 
4 .23811857710775980 (1) 


5.6744345899 1574145 (—1) 
1.33290773275989334 (0) 
2.38148248007005849 (0) 
3.72382664209342697 (0) 
5.37212395216187100 (0) 
7 .34193662826135234 (0) 
9 .65333213726122905 (0) 
1 .23323014070182379 (1) 
1.54128500654077 155 (1) 
1 .89402755758602723 (1) 
2 .29765957 156018692 (1) 
2.76101814474260749 (1) 
3. 29745092032584955 (1 ) 
3 .92898232533641479 (1) 
469768962767 103126 (1) 
5.71135140237534688 (1) 


2.31915524835569955 (0 ) 
5.12867199360117916 (0) 
9 .20089134897243093 (0) 
1 .53512814090706903 (1 ) 


1 .39445874535841327 (0) 
3 .00412262031589098 (0) 
5.16118127238140954 (0) 
7 .94175644134279969 (0) 
1. 14570496332457916 (1) 
1 .589356767 14650236 (1) 
2. 16116099678724373 (1) 
2 .95362536480182341 (1) 


1 .00157108495383152 (0) 
2. 14370613511225703 (0) 
3.64933321432680176 (0) 
5.54486813974946731 (0) 
.86076259241096344 (0) 
.06377654649934495 (1) 
.39325683951155931 (1 ) 
.78270063100842237 (1 ) 
. 24457538203894399 (1 ) 
.79955513560437942 (1 ) 
.48721707324949280 (1) 
.40889427543252506 (1 ) 


WN N RRR aI 


-82339164085635910 (— 1) 
.6700718367 1874001 (0) 
.83292171947257737 (0) 
.28441836819021891 (0) 
.0378687 1096889479 (0) 
.10954933413974021 (0) 
1 .05201017048502218 (1) 
1 .32960399123563184 (1) 
1.64718887086961996 (1) 
2 .00934999730530394 (1 ) 
2 .42235299428062549 (1 ) 
2.89511495806519383 (1 ) 
3 .44109430408020216 (1) 
4 082489527 13820481 (1) 
4 .86170549734729113 (1) 
5 .88737277583532394 (1) 


DOrnrs! 


.12749487761940805 (—3) 
.07769783328946458 (— 5) 
.27387146853583194 (—7 ) 
.65351819752401242 (— 10) 
.30506360143886766 (— 13) 
n = 3,N = 16 
.50981121009449315 (—2) 
. 65273224 (-—1) 
-48989313857906502 (0) 
.86124704489653486 (0) 
.30047554848272360 (0 ) 
-478196468569 15096 (— 1) 
.43843043208105701 (— 1) 
.37498472421622760 (—2) 
.442444693506 10863 (—3) 
.52340599558850768 (— 4) 
.50126530501356105 (—6) 
.06842629643597444 (—7 ) 
.92524335897222214 (— 10) 
.61863342780076823 (— 12) 
.54373278073215328 (— 15) 
.58180355287779219 (— 19) 
r= 4, N = 4 
.57222031198722368 (0) 
.37936031844206074 (1) 
.54028469036795246 (0) 
.38918132242164933 (—2) 
n=4,N=8 
.27551104534833669 (0) 
5597 1394358976204 (0) 
.98706281698755810 (0) 
42018619472331421 (0) 
17052700563446937 (—1) 
98824021635586082 (— 2) 
89788576623485837 (— 4) 
10804739994015747 (—6) 
n=4,N = 12 
59159548108807725 (—1) 
. 26582820687 688745 (0) 
.82396142854537 168 (0) 
75650688373216868 (0) 
7373097 2528196261 (0) 
.29328933322615071 (—1) 
19937171658365840 (—1) 
73699187701495750 (—3) 
284909545 16935571 (— 4) 
61140180598097782 (—6) 
2374658387 2488555 (—9) 
.01723739243414224 (—12) 
n=4,N = 16 
29739702663052232 (—1) 
49748665808338508 (0) 
.94612614281207374 (0) 
.43026705086785796 (0) 
.05687666231872829 (0) 
.90916131357755320 (0) 
.55584651819837666 (—1) 
.56009598652813629 (—1) 
.75191359650891542 (—2) 
.18227257300007 126 (—3) 
.58506311659349061 (—5) 
9 .50616079827638761 (—7) 
9 .38309759503402079 (—9) 
3 .62281000770432509 (— 11) 
3.75119190628700090 (— 14) 
4.01670675392835967 (— 18) 


WWWO HOF DD Ol ee Om COOWwWwe 


Cw co 


Row P Oe 


WOON NO ST O OI 1 


One 


> 00 DD 


1 .89270227412433568 (4) 
4 .64320234766113947 (4) 
1. 12370201079018605 (5) 
2 .82363592457219416 (5) 
8.41579093496417164 (5) 


1. 14816936994675722 (— 1) 
2. 14354639377120874 (0) 
1 .61220335260503062 (1 ) 
7 .70976931131275901 (1) 
2.80018019446226477 (2) 
8 .45668526756508810 (2) 
2 .24015331269839424 (3 ) 
5.38903830867048069 (3 ) 
1 .20653652088452424 (4) 
2.56137300941415792 (4) 
5 .23686431947 102649 (4) 
1 .04639104770985088 (5 ) 
2 .07677073163230993 (5) 
4. 18695241024345828 (5) 
8.93856486144031606 (5) 
2 .28131757123540439 (6) 


6 .682030295337 17859 (1) 
2 .32826105813133801 (3 ) 
3 .50699000772320933 (4 ) 
4.36118522676941444 (5) 


5. 14386973811936581 (0) 
52468188096643765 (2) 
74144726764228389 (3 ) 
243087 10943862802 (4 ) 
78085582369019424 (4) 
. 18617808357226932 (5 ) 
43383544036118813 (6) 
44715365181727399 (6) 


Se WOR Ree 


77831924001506431 (—1) 
7860767817 1097753 (1 ) 
.00823651354474138 (2 ) 
.98505168139673751 (3) 
69270306643015021 (3 ) 
73399567 482394664 (4 ) 
34831503047832052 (5) 
2730956 1007931643 (5) 
27919289421063792 (6) 
75999417306024579 (6) 
14970110208989482 (7 ) 
. 23832281907530249 (7 ) 


ee OO ee OOO Ot 


2836863765601 11066 (— 1) 
7.95547 187624795022 (0) 
8 .40596269993120727 (1) 
5 .39142636980744192 (2 ) 
2 .53782571556701110 (3) 
9 .67609960197269259 (3 ) 
3. 17018863503759819 (4) 
9 .28000231797978033 (4 ) 
2 .49553149851783729 (5) 
6. 29816124368594188 (5) 
1.51876543553876336 (6 ) 
3 .55901373314319768 (6 ) 
8 .25712602731401404 (6) 
1.94568396347320976 (7 ) 
4.87851881245905145 (7 ) 
1 .48732275042273460 (8 ) 
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TABLES OF 








ABSCISSAS AND WEIGHTS 


5 


TABLE—Concluded 
Zeros of wanes A Laguerre Weights 
=5,N=4 
z 91078636936770048 (0) 2. 96843023155522145 (1) 
6 .00000000000000000 (0 7 .00000000000000000 (1 ) 
1.03341103467862578 (1 ) 1.97469504802402056 (1) 
1.67551032838460417 (1) 5. ors te ie 


1 .78786047157082745 (0) 
3 .58823885764338781 (0) 
5 .92009002612269490 (0 ) 
8.86684640118214593 (0) 
1 .25435889444281568 (1) 
1.71405288833306908 (1 ) 
2 .30233426297284257 (1) 
3. 11295037859936706 (1) 


1.29748784000685508 (0 ) 
2 .586209436344067 19 (0) 
4.22697995386950259 (0) 
6 .25087547037452539 (0) 
8 .69028125595059872 (0) 
1. 15872437894981120 (1) 
1 .49995922575396278 (1) 
1 .90103582323823253 (1) 
2.37456904925278601 (1) 
2.9414476128886308+ (1) 
3.64163482053466644 (1 ) 
4 .57744569372735530 (1) 


1.01975621006729237 (0) 
2. ponrnnees ne cmmensg ad 
3 .30047946506096635 (0 ) 
4.85757645713809794 (0) 
6.71275559761595590 (0 ) 
8.88305877614446853 (0) 
1.13897106617410895 (1) 
1 .42597409258898484 (1 ) 
1.75281787111426743 (1) 
2.12414122811650920 (1) 
2.54627 162726269377 (1) 
3 .02820392701258409 (1) 
3.58350465381012185 (1) 
4 .23454277633176611 (1) 
5 .02404476674638389 (1) 
6 .06147860634799663 (1) 


=5,N = 

+ saaieamaineeens (0) 
3.3761326797 1511042 (1) 
5.12313725300345871 (1) 
2.54081265226598562 (1) 
4.537588981 16486575 (0) 

2 .74338324361631772 (—1) 
4.37283100034351893 (—3) 
8 .82759451633260886 (—6) 

n=5,N= 

1. 13334726805127857 (0) 
1.27102752975832523 (1) 

3 .59951248234067119 (1) 

4 .09672458660868638 (1 ) 
2.21781320782776474 (1) 
6.09556014277879247 (0) 

8 .58547769806400297 (—1) 
5.98481941810913418 (—2) 
1 .89538249456284785 (—3) 
2 .30998625486114025 (—5) 
7 .74406887601347796 (—8) 
3.01616108514055821 (— 11) 

= 5, N => 5 

3.51651938860095021 (—1) 
5.12310712750064588 (0) 
2.04008497899245158 (1) 
3.57964501907962966 (1) 

3 .32890537 160876807 (1) 
1.79120868308842187 (1) 
5.81751150687781879 (0) 

1. 15803031442409197 (0) 
1406413041 18632739 (— 1) 
1 .01865573672515247 (—2) 
4.21355482124525534 (—4) 
9 .27055945773008095 (—6) 
9 .67230928638172189 (—8) 
3 .93647895812592163 (— 10) 
4.29130683419748559 (— 13) 
4.84789862852962661 (—17) 


moh 


bo em ST et 


NOK WHE NOH NN 


mb OO ee OT OT be Or © 





Weight times Sumeuhdia of 


Zero = Cin exp (iy 


.45339079699888520 (2 ) 
-82400155444914586 (4) 
.07501194316356380 (5 ) 
-.07522522732126438 (7 ) 


.85855186820191489 (1) 
-22115795065788180 (3 ) 
-90898809485486296 (4 ) 
.80216730854013967 (5) 
-27185211458041413 (6) 
.626470463 14642504 (6 ) 
36187637191007579 (7 ) 
.91886557731016559 (8 ) 


14815342846979275 (0) 
.68784065874957923 (2 ) 
.46602657834341949 (3 ) 
12401457309494012 (4) 
31845757933745043 (5) 
56583037728656770 (5) 

.80546343478958589 (6 ) 
-07930634493047940 (7 ) 
-89331206067960144 (7 ) 
37452950475636027 (8 ) 
-06277099810570424 (8 ) 
28585454857795443 (9 ) 


7496 1600786600682 (—1) 
.88857731425075226 (1) 
.53386139750131562 (2 ) 
.60742858155172212 (3 ) 
7391395859309 1 106 (4) 
29124780054447953 (5) 
14310915694565138 (5 ) 
805704688963309 16 (6 ) 
761083593 19635732 (6) 
71023416056506905 (7 ) 
81909868948838928 (7 ) 
31349522839 105173 (8 ) 

. 535844404298 19067 (8 ) 
.§7149353332157447 (8) 
.82968933922063744 (9) 
.02381299844288556 (10) 
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Numerical Integration Over the Triangle 


By Gerald E. Bartholomew 


1. Introduction. Integration formulae for numerical integration over the simplex 
in n-space have been given by Hammer, Marlowe, and Stroud [1]. This paper 
describes a method for the evaluation of integrals over a triangle, based on a 
formula developed in [1]. The mid-points of the sides of the triangle are used as 
the points of evaluation in this method. 

Other possibilities may be used; for example, the trisection points of the medians 
that are not the centroid could be used as the points of evaluation. These points 
are symmetrically spaced and the weights, w; , again are equal. The centroid could 
also be used as a single point of evaluation with w; equal to the area of the tri- 
angle. As another modification to this method, one could calculate the appropriate 
mid-points of the reference triangle for a reasonable number of subdivisions, then 
retain them as constants rather than calculate each point as it is used. Recently, 
(1958) Mr. Walter Leffin carried out computations by the use of these methods. 


2. Mathematical Formulation. Consider the formula 
(1) [ [ 40 v) du dv = > w; fu; v;) 
T + 


where T' is a triangle, w; = A/3 for all 7, A being the area of T, and the points 
(u;, v;) are the mid-points of the sides of T. Formula (1) then gives a first ap- 
proximation to the integral. Now divide 7 into four triangles by connecting the 
mid-points of the sides. Evaluate (1) over each subtriangle T’ and sum. This 
gives a second approximation of the original integral. Iterate this process until the 
relation 


(2) | In naa < 6, 


is satisfied for some prescribed &, where J,, is the nth approximation to the integral. 

To find the desired mid-points of the sides of each triangle for successive sub- 
division, it was found convenient to subdivide a reference triangle R. Then the 
affine transformation 


U = Ayr + apy + ds 


(3) 

Vv = Anx + Any + des 
will give the proper points of T for the evaluation of (1). R is defined by the vertices 
(0, 0), (1, —1), (1, 1). The transformation (3) must then map R — T in the 
following manner 


(0,0) — (m, %) 
(4) (1, —1) — (t#, v2) 
(1, 1) — (us, v3) 





Received August 16, 1957; revised June 9, 1959. The machine calculation was supported 
by Army Ordnance. 
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where (1, :), (U2, V2), (Us, V3) are the vertices of 7. Substituting (4) into (3) 
and solving for the A;;, we obtain 





U2 + Uz — Qu Ve + v3 — 2r, Uz — Us 
a ee ada = etnies teat a 
(5) 
Vs — Ve 
a > a= uM a3 = UY. 


Since the area of FR is one unit, the absolute value of the Jacobian of (3) is the 
area of 7’. Note, also, that the area of each 7” is one-fourth of the area of 7. Hence 
the w; are easily determined for each approximation. 

The mid-points of the triangles in the previous (n — 1) subdivisions together 
with their vertices constitute the vertices of the triangles of the nth subdivision. 
Hence one could find the mid-points for the nth subdivision by use of the formulae 
of analytic geometry. This is undesirable because of the large amount of informa- 
tion that must be retained. For example, the 5th approximation has 768 mid- 
points, 408 of them are distinct. Therefore, a procedure is used that greatly reduces 
the amount of information that must be retained after each subdivision. This 
method uses three incremental values, A; , A2, A; , and the previously calculated 
mid-point, along with certain features of the triangle R to determine the next 
mid-point. 

Since R is symmetrical about the z-axis, by reflection across the x-axis one 
generates a point p,* = (z;, —y:) for every point p; = (2;, yi), where y ¥ 0. 
Also every mid-point not lying on the boundary of RF is a mid-point of two tri- 
angles. Hence one must check whether (a)a = |y| or (b)x = 1 to determine 
whether it is an interior point with respect to R. The increments A; , As , A; have 
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Fia. 1.—The sequence of point generation for first and second subdivisions. 
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the initial values 3, 1, 0 respectively. A; and A, are divided by a factor, two, for 
each successive subdivision, while A; is altered during each subdivision and is re- 
set to zero at the start of the next subdivision. 

The flow diagram in section 3 will show specifically how the increments are 
used to generate the p; . Figure 1 shows the order in which the points are generated 
for the first and second subdivisions. The points p; are actually calculated along 
straight lines originating at the z-axis with slope one. Hence the obvious check 
must be made to see whether the next point should be calculated along the same 
line, or whether the next point should initiate the next line. 
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3. Flow diagram. Figure 2 contains the flow diagram. The notation used in the 
flow diagram should be explained. The expression a — b means that a replaces b. 
For example, 1 — A; means that A; is set equal to one. | BE is the partial calcu- 
lation of J, that has already been made, or it is J, if the calculation is completed 
for the nth approximation. Dow-_vf is then Ip+. 


4. Results of Trial Calculations. A code was devised for the IBM 650 and 
several trial calculations were run. Table 1 gives the results that were obtained. 


It might be pointed out that formula (1) is exact in the first approximation for 


a quadratic function. This problem was suggested by Professor Preston C. Hammer. 
University of Wisconsin, 
Madison 6, Wisconsin 


1. P. C. Hamner, O. J. Martowe & A. H. Stroup, ‘Numerical integration over simplexes 
and cones,’’ MTAC, v. 10, 1956, p. 130-137. 








TECHNICAL NOTES AND SHORT PAPERS 


Efficient Method for Solving Atomic 
Schroedinger’s Equation 


By Sherwood Skillman 


1. Introduction. One of the basic numerical problems in atomic quantum 
theory is the solution of the Schroedinger’s wave equation for a spherically sym- 
metric potential. In practice, one is usually concerned with such potentials ob- 
tained by the Hartree-Fock self-consistent fields [1] or by the Thomas-Fermi-Dirac 
statistical field methods [2, 3]. 

This paper describes a highly efficient and rapidly convergent technique for 
solving the radial Schroedinger’s equation for an arbitrary atomic-like potential. 
The method has been programmed for the IBM 650 computer and the numerical 
results obtained are in good agreement both with pertinent experimental results 
(x-ray term levels) and with previous theoretical work [4]. 


2. Problem Definition. The Schroedinger’s wave equation in atomic units 


+ 


(r in Bohr radii, E in Rydbergs) is given by 








2 
d P(r) si q(r)P(r) for r = 0 
dr? 
where 
. ’ lil 1 
qr) = V(r) - E+ as ‘. 


In this equation P(r) is the radial wave function, V(r) is the potential energy, 
and F is the desired eigenvalue satisfying the boundary conditions P(0) = 0 and 
P(«) = 0. The parameter /, the angular-momentum number, is an integral num- 
ber 0, 1, 2, 3, ete. These numerical values of / are usually designated as s, p, f, 
d, ete. For each I, there is a set of eigenvalues satisfying the boundary values and 
labeled by an integer n. The lowest state for each / value is taken as n = | + 1, 
the principal quantum number, and for every value of n there is a distinct energy 
value EH. The various radial wave functions and associated eigenvalues are specified 
using the subscript nl. 


3. Integration Method. The integration mesh was made flexible for many dif- 
ferent type atoms by choosing it as a function of the atomic structure. The inte- 
grating variable « was made equal to r/u, where 


9 





al oT 
e* 4(6x2)8 
and z = the atomic number. It was found that Ar = .0025 was a satisfactory 


starting step and the mesh could be doubled every 40 points. The potential values 
were calculated for the same mesh using the Fermi-Thomas-Dirac equation and 
the desired boundary conditions. 


Received March 20, 1959; revised June 19, 1959. 
299 








FLOW CHART FOR SOLVING SCHROEDINGER’S EQUATION 
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METHOD FOR SOLVING SCHROEDINGER’S EQUATION 


CONVERGENCE OF TRIAL 
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QUANTUM STATE: Is 
TRIAL Es OE, 
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CONVERGENCE OF TRIAL 


EIGENVALUES 





ATOM: Germanium 
QUANTUM STATE: 2s 





“a t 


























TRIAL to. BE .. 
i - 150.000 | + 57.042 
2 |- 92.958 |- 7.956 
3 |- 100.914 |- 1.150 
4 |- 102.064 |- .015 
5 |- 102.079 |-— .000 0o2 
Fia. 3 
CONVERGENCE OF TRIAL EIGENVALUES 
P56 
ATOM: Cerium Laie, a 
QUANTUM STATE: 5d IV 
TRIAL San DE sg 
-2.22 000 |+.28 567 
2 - 1.93 433 |-.03 142 
3 -1.96 575 |-.00 126 
4 -1.96 701 |-.00 0003 
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A method applicable only to a linear second-order differential equation which 
lacks the first derivative was used to solve Schroedinger’s equation. This method 
introduced by Numerov [5] is approximately four times faster than the Runge- 
Kutta-Gill procedure and the truncation error is proportional to h°. 


4. Matching Procedure. Schroedinger’s equation is integrated outward to a 
point where P(r) changes from a traveling-wave behavior to a damped-wave 
behavior. The choice of this matching radius R,; is made dependent upon the value 
of q(r), the kinetic energy, becoming negative. This choice makes the program 


Ri 
quite general as it will handle any value of / and EH. The integral [ P*(r) drand 
0 


the logarithmic derivative of P(r) at r = R,, (P’(r)/P(r))out are stored. 

Then an inward integration is performed starting at R,. = (5 + 1)R,. At R, 
the final normalized wave function is essentially zero, thus satisfying the second 
boundary condition P(«) = 0. Since q in the neighborhood of R, is essentially 
a constant, a good approximation for P(R.) = ¢ **V2(®2)- This value is then 
used to start the inward integration, which is continued until R, is reached. The 
values of / P*(r) dr and (P’(R,)/P(R:))in are calculated and stored. 


“Ry 


5. Calculating Improved Eigenvalue. The difference at R, between the values 
of (P’/P)out and (P’/P) in is used to estimate what change in E, (AE), is neces- 
sary to obtain a better match. For each trial eigenvalue an outward and an inward 
integration is performed. A formula developed by Hartree (1) is then used to 
obtain AE as follows: 


{f° P*(r) aod, P*(r) dr| sn |(F), +(F) 
L PP(R1)oue PR) ) m . 
The value AF is then added to the trial EK and the process repeated 


until | AE/E | < 10°°. The radial wave functions are then made continuous at R;, 
and normalized. 





6. Flow-Chart of Procedure. The main features of the program are outlined 
in Fig. 1. 


7. Convergence of Trial Eigenvalues. The program has been used in applica- 
tions pertaining to the germanium, silicon, and cerium atoms and also the zinc- 
sulfide and gallium-arsenide crystals. It has been used for the seven quantum 
states 1s, 2s, 3s, 2p, 3p, 5d, and 4f. Some typical eigenvalues convergence rates are 
shown in Fig. 2, 3, and 4. 


Radio Corporation of America, 
David Sarnoff Research Center, 
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a A Hartree, The Calculation of Atomic Structure, John Wiley & Sons, New York, 
1957, 'p. 
e. , Encyclopedia of Physics, Vol. 36, Springer-Verlag, Berlin, 1956, p. 109. 
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4. R. Latrer, Physics Review 99, 1955, p . 510. 
5. M. G. Satvaport & M. L. BARON, ‘Womerieed Methods in Engineering, Prentice-Hall, 
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Explicit Formulae for 25 of the Associated 
Legendre Functions of the Second Kind 


By Dietrich Suschowk 


Since the present writer knows of no place where to find explicit expressions 
for some more of Legendres associated functions of the second kind than are given 
for instance in [1] or [2], he hopes the small collection of formulae at the end of 
this note might be of some use to others. 

The associated Legendre functions Q,*(z) of the second kind for |x| > 1 and 
nonnegative integer values of n, k can be defined by [1] 


k/2 a 
g,*(2) =< @ =D) | — 1)" log + +t) 








2°n! = dx t# 1 


id z+1 d,2 ‘ 

= me 8 it ‘— | , 

~ Dae ;| 1082 =) ae | 
These functions can be written as follows: 


Ga |sa Anu (a ) log = 24 + Ane (x | (n2 k) 








: (n — k)! 
(1) Q,(z) = (-1)! 
= 4 2) ‘ 
(@? — 1)*? 4 A nk (x) (n< k) 
where A‘'?(x) (i = 1, 2) is a polynomial and 


n+k fornZzk 


degree of Ank (z) = {n, degree for n < k ( ASe (x) = 0) 


n+k—-—1 fornZzk 


p a7, 
degree of Anz (2) = . —n-—1 forn<k 


(2) If the degree of A (x) is even (uneven), A (x) contains only even (uneven) 
powers of x. 

(3) The coefficients of A‘?(x) are integers. 

(4) For each fixed triple (n, k, i) and n = k the coefficients of AS? (x) have 
alternating signs and the coefficient of the highest power in A{) (x) is always 
positive and in A%? (x) always negative. Forn < k the coefficients of AS? (x) 
are all positive. 

(65) Ta2ek>O0, A (2) has a zero of nth order for x = 1. 

The formulae are as follows: Let 7 = 3 log (x + 1)/(a — 1). Then 


Qo'(x) = 


Q(z) = 2r—-1 

Q."(x) = 4 [(3a” — 1)r — 3a] 

Q;(a) = = [C152 — 9x)r — 152° + 4] 

Q(x) = a [(1052* — 9027 + 9)r — 1052 + 552] 
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Qe(z) = —— 
2-1 
Q(x) = FF (2 — 1)r — a] 
Q(x) = == ; [(3a*° — 32)r — 32” + 2] 
Q(x) = FT (15a — 182° + 3)r — 152 + 132] 
Q(x) = rae [(1052° — 1502° + 45x)r — 1052* + 1152* — 16] 
Qo (x) = =~ ; 
Qe) = =" 
Q(x) = + [(8a" — 6x" + 3)r — 30° + 5a] 
Q(x) = *~ [(150° — 302° + 150)r — 152" + 250" -8 
Q2(x) = Gp [(1052° — 2252 + 1352 — 15)r — 1052° + 1902* — 81] 
Q(z) = -gi 
Q%2) = -a 
Q(x) = — mate 
Q:°(x) = Ea [(15a° — 452* + 452° — 15)r — 152° + 402° — 332] 
Q(x) = : [(1052° — 3152° + 3152° — 105x)r 


(# — 1)" 


— 105° + 2802* — 2312" + 48] 
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Qo'(x) = me te 

ae) = Sp 

Q(x) = ayy 

Qs'(x) = at 

Q(x) = iy [(105a° — 4202° + 630z* — 4202” + 105)r 


— 1052’ + 4302° — 5liz* + 2312]. 


The above formulae have been calculated by use of the recurrence relations [2] 





Qe (2) 


(n — k)(n +k — 1)Q,*(x) — 2(k + 1) 


x 
? 


Wau i n (x) 





(n — k + 1)Qi41 (x) = (2n + 1)2Q,"(2) — (n + k)QLa(z) 
V2 —1 Q(x) = (n — k)aQ,' (x) — (n + k)QiAa(z) 
starting with 

a TEs 
V2-—1 , 


Mathematisches Institut der Technischen Hochschule, 
Miinchen, Germany 


Qo'(x) = 7, Qo'(x) s=— Q,°(x) =ar-—l. 





1. W. Magnus & F. OBERHETTINGER, Formeln und Sdize fiir die speziellen Funktionen der 
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2. A. Erpetyi, Editor, Higher Transcendental Functions, vol. 1, McGraw-Hill Book Co., 
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A Short Table of [  Jo(t)t" dt and [ ” Si(t)t* dt 


By I. M. Longman 


1. Introduction. In various physical applications, and in particular in geo- 
physics, there arises a need for the numerical evaluation of integrals of the type 


(1) tn [ Jo(t)g(t) at 
and 
(2) M - | Ji(tg(t) dt 


where x is a positive number, and where g(t) is, at least for sufficiently large t, a 
monotonically decreasing function of ¢. Here Jo(t), J:(t) denote the Bessel func- 
tions (of the first kind) of orders zero and one, respectively. 

A method of computation of (1) and (2) has been given by the author [1] for 
the case x = 0, and this method can also be applied to the evaluation of (1) and 
(2) for small values of x (such that the interval 0 — z is within the first few cycles 
of Jo and J; as the case may be), even if the integrals diverge for x = 0. The present 
paper presents an alternate method (which is easily extended to large values of x) 
for the case where g(t) is a function of the type that can be expanded as a series 
of inverse powers of ¢ 


(3) g(t) =a tat'+atf?+---, 
when ¢ is sufficiently large. 

2. Description of the Method. Suppose, then, that we wish to evaluate an 
integral of the type (1) or (2), where for simplicity we will suppose that an ex- 


pansion of the form (3) is valid when ¢ > x. Then our integrals can be expanded 
in the forms 


(4) L 


apf Jolt) atta f Jo(tyet att a [Ite a+ ---, 


(5) M ay | Ji(t) dt + a | Ji(t)t" dt + a. [ Ji(t)t* dt + ---, 


so that (1), (2) can be evaluated if we can evaluate integrals of the type 


(6) L, = / Jo(t)t™ dt, n = 0,1, 2,--- 
and 
(7) M, = / Ji(t)t™ dt, n=0@,1,2,-:-. 


We will see later that the series (4), (5) are quite rapidly convergent for x > 1, 
the rapidity of convergence increasing with x. However it should be noted that in 
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a number of applications (for example propagation in layered media, seismology 
problems) the integrals we need to evaluate have pt in place of ¢ in the arguments 
of the Bessel functions, where p is a positive number which may be as high as 20. 
After making the appropriate scale-factor transformation we obtain the integrals 
(6) and (7) with a factor p” in front of each integral. For this reason the inte- 
grals (6) and (7) are often required with high precision and at least ten place 
accuracy would be needed for such application. 

With regard to the integrals L, , M, , the only published tables known to the 
author are one of 


[ sett) at, x = 0(0.01)10 
0 
by Lowan and Abramowitz [2], from which we can obtain Lp as 
ly =1— | Jolt) at, 
() 


and one of L,, for x = 0 (0.1) 10 (1) 22, by Lowan, Blanch and Abramowitz, 
[3]. Watson [4] gives a table of the maxima and minima of 


[ Jo(t) dt 


and this table has been extended by the author [1], who has also pointed out an 
error in Watson’s table. Also Smith [5] and Lowan, Blanch and Abramowitz [3] 
give asymptotic expressions for L, . For large values of x beyond the range of the 
published tables we can use the series 


(8) [ Jolt) dt = W(x) + Jala) + Isla) + ++]. 


Also we have the obvious result 
(9) Mo == Jo(x). 


It is not the purpose of this paper to give extensive tables, but rather to present 
a general method for their computation, and to illustrate their use by way of an 
example. 


3. Method of Computation. Integration by parts show that the L, , M,, satisfy 
the following recurrence relations for a given value of z: 





(10) L, = —« "Ji(x) + (n + 1) Mays 

(11) M, = x “Jo(x) — Danas 

from which we deduce 

oa Lass = (n tes ws hae ~ in - 1)? 
and 

(13) es Jo(x) 4 Ji(a) M, 





n(n + 2)a"™ ~~ (n + 2)a"* ~ n(n + 2) ° 


Equation (12) has been given by Smith [5] where he uses it to obtain an asymp- 
totic expansion for L, . Thus starting from LZ» and L, , which are either obtained 
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from existing tables, or must be calculated, we can use (12) to compute successive 
L,’s. Furthermore the form (12) shows that for even moderately large x, L, tends 
to zero very rapidly as n increases. This makes the calculation of the series (4) 
a very rapid process. 

With regard to the M,’s, starting with the relations 


My = Jo(2) 
M, = Lo a J1(x) 
M, = (Ly + a 'Ji(x)) 


(which are obtainable from (10), (11)), successive M,’s can be computed by 
means of (13). Some L,’s and M,’s have been computed in this way, and these 
are given in tables 1 and 2 

With regard to Ly and L, , these are obtainable (by interpolation if necessary ) 
from existing tables ({1], [2], [3], [4], (5]), or, for large x, by means of asymptotic 
expansions. The interpolation can easily be carried out accurately using Taylor’s 
theorem. Examples of such interpolation are given in Longman [1]. Asymptotic 
series for L, are given by Smith [5] and by Lowan, Blanch and. Abramowitz [3], 
while for Io we can apply Smith’s method to obtain the result 





” Jolt) dt ~ Jo(x) [t-* wb By ESS 1334, 


x xz‘ 
rs V2s 
oh ra tel, 


and for large x this can be used to compute Ly using the-known asymptotic series 
for Jo(X), Ji(X). 


z= 


(14) 





4. Tables. Tables 1 and 2 were computed from equations (10), (11), (12) and 
(13) and existing tables [2] and [3] of Z» and LZ, by the use of an IBM 709 computer 
at the Western Data Processing Center of the University of California, Los Angeles. 


5. Example of the Use of the Tables. Suppose we wish to compute 


r= | wolt/(é + 1) a 
10 
Using the expansion 


2 


(@+1it]|=f-+r---., (t > 1) 
we have 

T=k,-K+he—[s+-:::, 
from which we obtain 


I = —0.001042. 


Institute of Geophysics, 
University of California, 
Los Angeles, California 
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1. I. M. Loneman, ‘‘Tables for the rapid and accurate numerical evaluation of certain 
infinite integrals involving Bessel functions,’’ MTAC, v. 11, 1957, p. 166 
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Numerical Integration of l "7D, () ds (2) a dx 


By Gerard P. Weeg 


1. Introduction. This integral arises in a certain thermo-elastic study, and 
although it is similar to many evaluated in Watson’s “Bessel Functions,” this 
integral is not listed. Numerical integration techniques were applied. By means of 
Simpson’s rule with the limits of integration set at 0 to 12 and making use of 400 
points, the above integral was evaluated to four decimal places for seven different 
combinations of 7 and é. 

The same integral was evaluated by the Gauss-Laguerre formula. Although the 
integrals in question can be expressed in terms of complicated elliptic integrals, the 
natural method of evaluation, to avoid lengthy analysis, is by quadrature. It is 
deemed that the experience with the Gauss-Laguerre formula is worth reporting, 
even though the integrals are of a special nature. Using the Gauss-Laguerre formula 
with only six points, the results agreed with Simpson’s rule entries to three decimal 
places, although the upper bound of the truncation error associated with the 
Gauss-Laguerre formula was greater than unity. 


2. Numerical Results. The integral under consideration is 


(1) I(n, &,n) = I e Jo (") Ji (2) xz” dz 


where 7 and é are integers and n is 0 or 1. A careful analysis of Simpson’s rule 
applied to I(n, —, 0) indicates accuracy to within a unit in the fifth decimal place 
in the values listed in Table 1. It is possible to show* that 


1 
(2) I(n, &,0) + I(n, & 1) = —; [Qualx) — 1Q(zx)] 
at 
where « = (1 + & + 7')/2n and Q,,(x) is the Legendre function of the second 
kind. Use of (2) then verifies that the values of J(n, —, 1) are also correct to the 
number of places given in Table 1. 

















TABLE 1 
gE 2 I@, &, 0) | T(n, €, 1) | Maximum Error in J(m, =, 0) 
1 1 17868 .33880 | 1.12 x 10-5 
1 3 .01741 . 15950 8.37 xX 10-8 
1 6 .00229 .08246 7.97 xX 10° 
3 1 . 13463 . 15486 2.48 X 10-5 
3 3 .05948 .11742 G.22 =< ip 
6 1 .07850 | .08167 La x we 
6 6 | .02954 | .05887 7.97 X 10-* 





Received May 12, 1958; revised April 13, 1959. 
* The author is indebted to Dr. Gertrude Blanch for this relationship. 
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I(n, &, 0) was also evaluated for § = 7 = 1 by means of the Gauss-Laguerre for- 
mula. Using four points the results agree to two decimal places with those in 
Table 1. Using six points the results agree to three places. In both cases, the trunca- 
tion error found exceeded unity. 


Computer Laboratory, 
Michigan State University, 
East Lansing, Michigan 


1. G. N. Watson, A Treatise on the Theory of Bessel Functions, The Macmillan Company, 
New York, 1945. 











Calculation of Gamma Functions to 
High Accuracy 


By M. E. Sherry and S. Fulda* 


The two constants, I (4) and T (4) have been calculated to 35 decimal places 
on the Cambridge Computer by means of an interpretive routine that treats float- 
ing-point numbers of 37 significant digits [1, 2]. 

Stirling’s asymptotic expansion for In T (x) can be written 


C, 


g?r-l P 





(1) InT(xz) ~ («2 — 3) nx —2+43l2r4+ >> 
r=l 


With the remainder of the series made very small, the accuracy of In I'(2) 
depends primarily on the accuracy of In 27, In x, and the C,. Uhler [3, 4] has pub- 
lished the In x for all primes through 101, as well as the In z, to more than 100 
significant figures. Uhler has also calculated the C, to over 100 significant figures 
(5). 


The recursion formula for the Gamma function 


(2) T(z +1) = 2I(z) 

can be extended to 

(3) l(c + 2) = (2),l(z) 

where (2x), is the Pochhammer-Barnes symbol 

(4) (t)n = (u@ +n —1)(a2 +n — 2)--: (4 1)z. 


If we take the logarithm of Eq. (3), 


n—l 
(5) InT(2 +n) = n(x) + Vln (a+ J). 

j=0 
Substituting Eq. (1) into Eq. (5) gives, after solving for In T(x), 
(6) ef (x) ~ Mz) + = 


rat (x + nn)?!’ 


where 


n—l 
(7) A(z) = (a@ +n — 4) In(ex+n) — (x4+n) + 4iln27 —- > In (a +). 
j=0 
The term A(x) is calculated quite easily by hand, and it is not difficult to calculate 
In I(x) from Eq. (6) with a digital computer. 
To calculate ['(2), let 


= T(x) = event@—v 


Received December 31, 1958; revised May 4, 1959. 
*§. Fulda, Bell Telephone Laboratories, Murray Hill, New Jersey, co-authored this 
paper while on temporary duty at the Air Force Cambridge Research Center. 
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where 
(9) vy = +0.1n 


is the exponent of e nearest in value to In T(x). Values of e*°'", n being any integer 
from 0 to 10, have been tabulated by Van Orstrand [6]. 


(10) Since |InT(rz) — ¥| S 0.05 


it will always be possible to calculate e”™"@ 


series expansion using very few terms. 
These techniques led to the following values: 


by means of the exponential power 


I'(4) = 2.6789 38534 70774 76336 55692 94097 46776 


(2) = 1.3541 17939 42640 04169 45288 02815 45138 
In P'(4) = .98542 06469 27767 06918 71740 36977 96139 
In P(2) = .30315 02751 47523 56867 58628 17372 01104. 


It was possible to examine round-off and truncation errors at each step in the 
calculation. The final relative error was less than +3 X 10° in each case. 
As a final independent check, the values of T'(4) and T'(2) were put into 


9 
r(3)r@) = —. 
(11) (3)T() V3 


The error analysis was consistent with this identity. 


Computer & Mathematical Sciences Laboratory, 
Air Force Cambridge Research Center, 
Bedford, Massachusetts 


1. B. Zonpex, “The values of P'(4) and (3) and their logarithms accurate to 28 decimals,’ 
MTAC, v. 9, 1955, p. 24-25. 
2. Air Force Cambridge Research Center parece, aw Computer Interpretive Routine 
for ‘eae: Precision Numbers, Series 3, TN-59-155, 1959. 
. H. S. Unuer, ‘Natural logarithms of small primary numbers,” Proc., Nat. Acad. Sci., 
v. 2 1943, p. 319-325. 
H. 8. Unter, ‘‘Log z and other basic constants,’’ Proc., Nat. Acad. Sci., v. 24, 1938, 
2330. 
m 5. H. S. Unter, ‘‘The coefficients of Stirling’s series for log Tf (x),’’ Proc., Nat. Acad. 
Sci., v. 28, 1942, p. 59-62. 
6. C. EB. Van Gasmnane, “Tables of the exponential function and of the circular sine and 
cosine to radian argument,’’ Mem., Nat. Acad. Sci., v. 14: 5, 1921 (Tables IV and V). 











REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


49 [F, Z|.—Martin Davis, Computability and Unsolvability, McGraw-Hill Book 
Company, Inc., New York, 1948, xxv + 210 p., 23 em. Price $7.50. 


The following is a quotation from the author’s preface. 

“This book is an introduction to the theory of computability and noncomputa- 
bility, usually referred to as the theory of recursive functions. This subject is con- 
cerned with the existence of purely mechanical procedures for solving problems. 
Although the theory is a branch of pure mathematics, it is, because of its relevance 
to certain philosophical questions and to the theory of digital computers, of poten- 
tial interest to nonmathematicians. The existence of absolutely unsolvable problems 
and the Gédel incompleteness theorem are among the results in the theory of com- 
putability which have philosophical significance. The existence of universal Turing 
machines, another result of the theory, confirms the belief of those working with 
digital computers that it is possible to construct a single ‘all-purpose’ digital computer 
on which can be programmed (subject of course to limitations of time and memory 
capacity) any problem that could be programmed for any conceivable deterministic 
digital computer. This assertion is sometimes heard in the strengthened form: any- 
thing that can be made completely precise can be programmed for an all-purpose 
digital computer. However, in this form, the assertion is false. In fact, one of the basic 
results of the theory of computability (namely, the existence of nonrecursively enu- 
merable sets) may be interpreted as asserting the possibility of programming a given 
computer in such a way that it is impossible to program a computer (either a copy 
of the given computer or another machine) so as to determine whether or not a 
given item will be a part of the output of the given computer. Another result (the 
unsolvability of the halting problem) may be interpreted as implying the impossi- 
bility of constructing a program for determining whether or not an arbitrary given 
program is free of ‘loops.’ 

“Because it was my aim to make the theory of computability accessible to per- 
sons of diverse backgrounds and interests, I have been careful (particularly in the 
first seven chapters) to assume no special mathematical training on the reader’s 
part.” 

The author succeeds admirably in presenting a readable, motivated, yet not 
verbose, exposition of the theory of recursive functions. The not mathematically 
mature reader will, however, find the reading slow, but this is probably intrinsic to 
the subject matter. The notion of a Turing machine is taken as basic, and the re- 
viewer is gratified to note that the author sharply distinguishes what a Turing ma- 
chine is from how it behaves. Up to the appearance of this book hardly any two 
articles on Turing machines employed exactly the same concepts of a Turing ma- 
chine and functions computed by such a machine. This made it necessary to check 
that results established for one concept held for another. Since the concepts (men- 
tioned above) employed by the author seem highly satisfactory it is to be hoped 
that the tendency to introduce even more variants will be attenuated. 

While the book appears suitable for the classroom, this usefulness would have 
been greatly enhanced by the inclusion of exercises which, in part, would indicate 
alternative lines of development. The reader familiar with programming digital 
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computers will note that construction of Turing machines is an activity similar to 
programming. One should be wary, however, in interpreting results concerning 
Turing machines as results concerning digital computers, for this implies, in partic- 
ular, a willingness to endow the (idealized) digital computer with (potentially) infi- 
nite memory. For some purposes, at least, it is more appropriate to take as a mathe- 
matical model of a digital computer the notion of a finite (memory) automaton. The 
significance of this change in mathematical model is indicated by the fact that the 
analogue of the holding problem for finite automata is recursively solvable. This 
may be interpreted as implying the possibility of constructing a program for a Uni- 
versal Turing Machine for determining whether or not an arbitrarily given program 
for a digital computer is free of “loops.” 

The chapter headings are: Computable Functions, Operations on Computable 
Functions, Recursive Functions, Turing Machines Self-Applied, Unsolvable Deci- 
sion Problems, Combinatorial Problems, Diophantine Equations, Mathematical 
Logic, The Kleene Hierarchy, Computable Functionals, and The Classification of 
Unsolvable Decision Problems. 


C. Eicor 
The University of Michigan, 
Ann Arbor, Michigan 


50 [G, H, W, X, Z].—MaruematiscuHes Lasor DER TECHNISCHEN HocHSCHULE 
Wien, MTW Mitteilungen, Vol. 5, 1958, 380 p., 24 em. Price DM 15. 


Of the 380 pages in this volume, two hundred are devoted to contributed articles, 
the remainder to various progress reports, reviews, news, etc. The relative propor- 
tions of these features are about the same as in earlier volumes [1]. The “literature 
reports” (references to new publications, arranged by subject matter) seem to 
have been discontinued after March, 1958. 

Of the contributed articles, some are accounts of new research results, some are 
expository surveys, and some are reports on meetings or talks. The following are 
likely to be of interest to the computer mathematician. Most of the articles are in 
German; the titles Uelow are translated. 

W. KnoGpet: A Transportation Problem. The usual algorithm for solving trans- 
portation problems is applied to the Austrian sugar industry, using standard 
punched-card machines. 

W. SPINDELBERGER: Determination of complex roots of algebraic equations on 
small electronic computers. 

W. GrésneErR: Computation of powers of matrices. The mth power of ann X n 
matrix A, where m is considerably larger than n, is expressed as a polynomial p,, (A) 
in A of degree at most n; more precisely, of degree at most d, where d < m is the 
degree of the minimal polynomial of A. The coefficients of p,, can be found by recur- 
rence with respect to m, if the coefficients of the minimal polynomial are known. 

H. RecHBERGER and H.SrqueEnz: Basic problems in simulation of transients in 
reactors. A discussion of analog and mixed digital-analog systems for this problem. 

K. Zuss: “Feldrechenmaschine.” This is the name given by the author to a pro- 
posed logical design for a digital computer with single commands causing arithmetic 
operations on “fields of data,” such as vectors or matrices. 

E. Ecrervary: Remarks on the transportation problem. A theorem of the author 
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has been used by H. W. Kuhn and M. Flood to develop an algorithm for the As- 
signment Problem of linear programming. In the present paper the author modifies 
the method of Kuhn and Flood for transportation problems with integer coefficients; 
these can be considered as special cases of assignment problems. 
F. G. Tricomr: Active problems in the theory of ordinary differential equations. 
This is a survey of the recent literature. 
Franz L. Aur 


National Bureau of Standards, 
Washington, D. C. 


1. MTAC, v. 13, 1959, p. 67-68 (Review 10). 


51 [G, S].—B. J. Sears & M. G. Raprke, Algebraic Tables of Clebsch-Gordan Coeffi- 
cients, Report AECL No. 746, 1954, II + 10 p., 27 em. Available from Scientific 
Document Distribution Office, Atomic Energy of Canada Limited, Chalk River, 
Ontario, Canada. Price $.50. 


Prior to the development of quantum mechanics, the Clebsch-Gordan series 
played an important role in the theory of invariants of binary forms and elsewhere 
in group theory. In quantum mechanics, the coefficients in these series give the 
matrix components needed for the coupling of angular-momentum vectors. If two 
angular momenta having quantum numbers 7; and j2 and components m and m2 
(all integral or half-integral) are coupled to give resultants 7 and m, the matrix 
component (71 J2 m1 M2 | j m) is a Clebsch-Gordan coefficient. 

The algebraic formula for these coefficients involves a complex sum; hence, to 
facilitate computation of the numerical values needed by theoretical physicists 
studying atomic and nuclear structure, algebraic tables were devised. The present 
volume contains six tables in which the values of 7; are successively fixed at 3, 1, %, 
2, $, and 3. The rows of the tables correspond to the permitted values of 7, from j2 + 
ji to J2 — ji (J2 = ji), while the columns correspond to the values of m; from +); 
to — ji. The entries -re comparatively simple algebraic expressions in j and m. 

Only the table for 7: = $ represents a new calculation from the definitions. The 
first four tables are rearrangements of those appearing in Condon and Shortley, 
Theory of Atomic Spectra (1935), while the last was computed by Falkoff, Collady, 
and Sells in 1952. The rearrangement by the present authors involves expressing 
the entries in terms of j and m instead of j2 and m (actually j; and m in previous 
tables where 72 was fixed), and leads to slightly simpler algebraic expressions but 
no essential reduction in time required for numerical substitution. 


GEORGE SHORTLEY 


Booz, Allen Applied Research, Inc., 
Bethesda, Maryland 


52 [G, X].—Franz E. Houn, Elementary Matrix Algebra, The Macmillan Com- 
pany, New York, 1958, xi + 305 p., 24 cm. Price $10.00. 


For several years the reviewer has taught beginning graduate courses in the com- 
putational methods of linear algebra, for which it is essential that the student have 
a strong knowledge of concrete linear algebra, of finite-dimensional linear analysis, 
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and of the analytic geometry of hyperplanes and quadric surfaces in n dimensions. 

Sven reasonably good students who have studied Birkhoff and MacLane [1] seem 
astonishingly unable to use the basic theorems of linear algebra in the concrete 
analysis of actual computations. As an example, students need a working familiarity 
with the Jordan canonical form of a linear transformation and its relation to the 
biorthogonal system of eigenvectors and principal vectors of the transformation. 
Without this knowledge it is quite impossible to understand, for example, the inner 
working of the Liebmann process or successive overrelaxation. Another important 
need is a familiarity with the principal vector norms of numerical analysis—the 
pth-power norms for p = 1, 2, ©—and the associated matrix norms. 

The reviewer is often asked to recommend suitable textbooks for remedial work 
in this area. Such books as Perlis [2] seem a little too elementary, and skip the analy- 
sis and geometry in favor of the algebra. Perhaps the best reference has been the 
first chapter of Faddeeva [3], until now available only in Russian or in an informal 
translation. And Faddeeva, good as it is in algebra and analysis, lacks the geometric 
point of view and has no exercises. 

The reviewer therefore welcomed the arrival of Hohn’s book as a new textbook 
in this field, not written by a pure algebraist (the author is an associate professor 
of mathematics at Illinois), and directed to future applied mathematicians who are 
at the junior-senior level in a university. The book is now reviewed from the stand- 
point of numerical analysis only, although it is directed to a much wider audience. 

The reviewer’s conclusion is that Hohn has written an unusually good exposition 
of the material covered, and that the exercises are especially good. Students will find 
it very readable, and the author presupposes only college algebra. On the other 
hand, this is basically another book on algebra, with very little analysis or geometry. 
Moreover, the algebra, while adequate on many topics (like Hermitian forms), 
seems deplorably lacking in any mention of the canonical forms under similarity. 
For the needs of numerical analysis as the reviewer sees them, the book would have 
to be supplemented by other sources of linear algebra and analysis. 

There are nine chapters and three appendices. Chapter 1 introduces the reader 
to matrix notation. Chapter 2 defines a determinant (as the sum of n! terms), and 
gives basic properties. Chapter 3 deals with solving linear systems and inverting 
matrices. Chapter 4 covers rank and equivalence of matrices, with a proof that the 
rank of AB cannot exceed the rank of A or B. The concept of a number field (subfield 
of the complex numbers) is used in most further work. 

Linear dependence is the topic of Chapter 5, and the fundamental existence 
theorem for a linear algebraic system is stated here in terms of rank. Chapter 6 is 
devoted to vector spaces and linear transformations. We get to Sylvester’s law on 
the nullity of AB in terms of the nullities of A and B. The Jordan canonical form is 
apparently nowhere even alluded to, nor is any other canonical form under simi- 
larity, and this disappoints the reviewer very much. In Chapter 7 we meet some 
analysis during the development of unitary and orthogonal transformations. We 
again meet the important existence theorem for a linear algebraic system AX = B, 
this time in terms of the orthogonality of B to the nullspace of A* (called the tranju- 
gate of A), but the nature of the general solution is not mentioned. 

In Chapter 8 the characteristic equation is discussed, with the Cayley-Hamilton 
theorem, and we do learn about diagonalizing a Hermitian matrix. In Chapter 9 are 
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discussed quadratic and Hermitian forms, definiteness, and even Cochran’s theorem 
on representing a semidefinite form as a sum of squares. 

The appendices are devoted to the >> and J] notations for sums and products, 
to the algebra of complex numbers, and to isomorphism as a concept. There follows 
an excellent 7-page bibliography of books on linear algebra and its applications in 
pure mathematics, numerical analysis, economics, psychology, electrical engineer- 
ing, chemistry, etc. 

The author consistently denotes matrices and column vectors by capital letters, 
but frequently uses component equations. He thus happily seems to make the nota- 
tion the reader’s ally instead of an enemy. 

One should mention that a little desk-machine numerical analysis is contained in 
the book. First, there is an elimination method for evaluating a determinant (called 
the “sweep-out” process), with a numerical example of order 3 and a reference to 
Dwyer [4] for more details. The other exposition is of Gaussian elimination for 
solving a linear algebraic system, and an abridged form thereof, with two more nu- 
merical examples of order 3. For the student interested in round-off errors, there is a 
reference to von Neumann and Goldstine [5], but not to its sequel [6]. There is no 
mention of automatic computers as a tool for solving linear systems! 


GeorGeE E. ForsytTHE 
Stanford University, 
Stanford, California 
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53 [H, I, X].—W. J. Cunnrneuam, Introduction te Nonlinear Analysis, Electrical 
and Electronic Engineering Series, McGraw-Hill Book Co., Inc., New York, 
1958, ix + 349 p., 23 cm. Price $9.50. 


This book is the well planned and well written outgrowth of a graduate course in 
electrical engineering. The subject matter is the practical analysis of nonlinear ordi- 
nary differential equations. A list of the chapter headings follows. 

1. Introduction 
. Numerical Methods 
. Graphical Methods 
. Equations with Known Exact Solutions 
. Analysis of Singular Points 
. Analytical Methods 
. Forced Oscillating Systems 
. Systems described by Differential-difference Equations 
. Linear Differential Equations with Varying Coefficients 
. Stability of Nonlinear Systems 
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The 47 pages of the second and third chapters give a short and practical account 
of methods for desk calculation and for graphical solution. The remaining chapters 
frequently refer to these methods. The book contains many excellent figures. The 
exposition is bolstered with carefully worked out examples for each new idea, as is 
appropriate for a good textbook. The homework problems are given in 16 pages at 
the end of the book and are arranged according to the appropriate chapter. A bibli- 
ography, with comments on relevance to the book, is given in seven pages. A five- 
page index completes the book. 

E. I. 


54 [H, S).—A. Jerrrey, “Tables of characteristic roots and functions arising in 
neutron transport problems,” 15 p. Deposited in the UMT File. 


Application of the spherical harmonic method to neutron transport problems 
[1], [2] requires the knowledge of the roots v of a certain determinantal equation 
(equation (1)), and of the functions G, (v) of those roots (given in equations (2) 
and (3) below). 
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(2) (n+ 1)Gyus(v) + G,(v) + nG,1(v) = 0,n = 1,2,3,---, N 
4s 9 
(3) Go(v) = 1, Gi(v) = (¢ . 





Tables have been constructed listing the positive roots v of equation (1) for N = 
9, 11, 13 and for c in the range 0(0.05)0.9(0.01)1.0. An abbreviated tabulation of 
the functions G,(v;) of equations (2) and (3) is also given for argument c in the 
range 0(0.1)0.8(0.05)0.95, 0.99, 1.0. These 7D data have been deposited in the 
MT AC repository for Unpublished Mathematical Tables. 
K.P. 


1. J.C. Mark, The Spherical a Method, Parts I and II, Atomic Energy of Canada 
Limited, Reports CRT-340 and CRT-338. 
2. B. Davison, Neutron Transport Theory, Oxford University Press, 1957. 


55 [I].—H. E. Sauzer, Tables of Osculatory Interpolation Coefficients, National Bu- 
reau of Standards Applied Math. Series, No. 56, U.S. Government Printing 
Office, Washington, D.C., 1959, xi + 25 p., 26 em., Price $.30. 

Let f(x) be a function defined over a region X, with at least 2n bounded deriva- 
tives in X. Let fi = f(x,); fi’ = (df(x)/dz),..,. The well-known Lagrange-Hermite 
formula for osculatory interpolation has the form 

[n/2] 


f(z) = flto+ ph) = De {AM Pdf + BPM P)f'} + Ren 


i=—[(n—1 











322 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


where 


 haned (w)h™ j=[n/2] 


Re a 
‘ (2n)! — jm—t(n—1)/21 


(p—j)’. 
In the above [s] denotes the greatest integer in s, and w is an inner point of X. 

An inspection of the form of the remainder shows that for comparable values 
of the derivative of order 2n, the formula given avove is more accurate than the 
corresponding ordinary Lagrangian formula of degree 2n (which involves the same 
amount of work). The above formula is therefore more advantageous in cases 
where f(x) is required, and where both f(x) and f’(x) are known or easily obtainable 
at points z;. The author tabulates the coefficients A; and B; over the following range: 
n = 2,3, p = —[(n — 1)/2] (.01) [(n/2)], exact; n = 4, 5, to 9D, same range and 
interval. 

The table should prove to be very useful in application. The author gives the in- 
exact values to within 1.5 units in the last place, but adds: ‘““This accuracy is not 
absolutely guaranteed, owing to the large amount of extra work that would be 
necessary to check it completely.”” Considering the author’s past performance, a 
table issued under his name has indeed a high probability of being accurate. It is, 
nevertheless, a sobering fact that rigorous methods of checking tables are nowadays 
considered too expensive to be worth while! 

GERTRUDE BLANCH 


Wright Air Development Center, 
Wright-Patterson Air Force Base, Ohio 


56 [I, X).—L. Fox, The Numerical Solution of Two-Point Boundary Problems in Ordi- 
nary Differential Equations, Clarendon Press, Oxford, 1957, xi + 371 p., 24 cm. 
Price $9.60. . 


This volume is the first in a series of monographs on numerical analysis planned 
by E. T. Goodwin and the late D. R. Hartree. In it, the numerical solution of one- 
dimensional boundary problems by finite difference methods receives the most com- 
prehensive treatment which has as yet appeared in print. The emphasis is on purely 
numerical methods of universal applicability; ‘semi-analytical’ techniques such as 
the methods of Ritz and Galerkin are not treated. 

The author modestly disclaims all pretensions to rigorous or highbrow mathe- 
matics. According to his own words, the book is intended for the practical computer 
and for the student seeking a career in computation. Yet the book can hardly be 
dismissed as mathematically trivial; it would take a competent analyst to prove, 
say, the convergence of some of the algerithnmas proposed in the book, or to establish 
rigorous error bounds. Students cf numerical analysis looking for topics for Ph.D. 
dissertations will find ample raw material here. 

The recurrent theme of the author’s treatment of boundary problems is the sys- 
tematic use of the difference correction. This is mathematically (although not nu- 
merically) equivalent to improving a crude first approximation to a solution by 
taking into account the error terms in the finite difference approximation to the 
differential operator. As the author puts his basic philosophy (p. 35): ““Most writers 
neglect the difference correction, using more accurate formulae for the derivatives, 
or a small enough interval, so that the error is either negligible or kept within reason- 
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able bounds. In this book we shall work with a reasonably large interval and take 
full account of the difference correction. A large interval involves a correspondingly 
small number of algebraic equations and a proportionately small amount of labour: 
the inclusion of the difference correction yields high accuracy at this large interval.” 
He also says in the preface: “Throughout I have taken the view that truncation 
errors in finite difference equations should not be tolerated.’”’ While it would be easy 
to take issue with these principles, if taken too literally, it cannot be denied that 
the author makes a strong case for a reasonable interpretation of them. 

After introductory chapters of a general nature on the algebra of finite difference 
operators and on the solution of algebraic equations, the author considers sec- 
ond-order problems in Chapter 4. There follows, somewhat surprisingly, a chapter 
on the solution of initial-value problems by boundary-value techniques. Further 
chapters concern equations of higher than second order, eigenvalue problems, accu- 
racy and precision, and miscellaneous methods. The discussion proceeds almost 
wholly by examples, all fully worked (on a desk computer). Many of the examples 
illustrate special difficulties and points of interest, such as the treatment of singu- 
larities, changes of the basic interval, and boundary conditions at infinity. Some of 
the famous tricks due to classical British numerical analysts, such as L. F. Richard- 
son’s deferred approach to the limit and A. C. Aitken’s 6*-method, are discussed in 
detail. 

The author writes a terse, but at the same time fluent and highly readable 
English. The Clarendon Press has with this volume once more demonstrated its 
tradition of superb mathematical printing. 

Perer HEnrRIcI 
Numerical Analysis Research, 
University of California, 
Los Angeles 24, California 


57 [K].—Forman 8. Acton, Analysis of Straight-Line Data, John Wiley & Sons, 
Inc., New York, 1959, xiii + 267 p., 24 em. Price $9.00. 


As the title states, this is a book fitting a straight line to a set of data. Every 
conceivable set of a priori situations is assumed, and then the techniques for making 
the fit are fully and clearly discussed based on these assumptions. All the neces- 
sary side problems are considered and illustrated by examples. 

The author’s preface states very well why the book was written, and I would 
like to quote a few lines. 

“This book was written from the conviction that we need more detailed exposi- 
tions of classical and modern statistical techniques than are now possible in a gen- 
eral text surveying the entire field. Accordingly, I have selected the topics most per- 
tinent to the engineer or physical scientist when he deals with data containing one 
or more lines entrapped within his experimental variability. He may seek to reveal 
the line, or he may desire to remove it in order better to examine the residual fluc- 
tuations. In either event, he must encounter the philosophy and arithmetic of the 
analysis of variance—a philosophy that all too often has been shoved off as a step- 
child because its mathematical theory is rather tedious and pedestrian. 

“Tn particular, I wanted to set down a strong plea for the use of these analytical 
techniques by the experimentalist himself, rather than by a professional statistician. 
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The structure of the experiment determines the proper method for analyzing its 
data, and the experimenter himself knows that structure best. Since his statistical 
background is apt to be shaky, he needs to be guided—at times in considerable 
detail—and he usually prefers concrete examples from which general principles may 
be inferred. I have therefore tried to introduce each subject with data from a real 
experiment—examining them at first broadly, and then again with an increase in 
detail that occasionally they do not deserve. When I have felt that they will not 
mislead, I have sometimes relied on explanations that are heuristic, but in general 
I have tried to emphasize the essential questions that are posed by the mathematical 
models our experimenter may choose in analyzing his data. It is with the influence 
of the model on the nature of the extractable information that I have been chiefly 
concerned.” 

The author also states that he did not intend the book for ‘‘classroom use,” but 
rather as a detailed reference book. It is quite clear that Professor Acton has done 
a very excellent job and carried out fully the goal he set for himself. 

After a brief introduction in Chapter 1, the author discusses in Chapters 2 and 3 
the standard straight-line-fit problem in every way possible. Some needed informa- 
tion about the bivariate normal is presented in Chapter 4. Chapter 5 then considers 
the non-standard straight line fit when both variables are subject to error. Chapter 
6 takes up the straight line fit in the analysis of variance and presents (for the first 
time in a text, to my knowledge) information of great value when ‘‘doing” an analy- 
sis of variance. Following Chapter 7, which isdevoted to non-linearfit and orthogonal 
polynomial fit, the last three chapters round out the text by treating questions which 
usually come up in connection with least-squares problems, namely, transforma- 
tions, rejection of extreme values, and cumulative frequency fitting. 

Having this book available will really be a great boon for any statistician who is 
in daily contact with non-statistical scientific people. The nonstatistician would 
profit from exposure to it. It is an excellent book. 

Harry WEINGARTEN 
Special Projects Office, 


Department of the Navy, 
Washington, D. C. 


58 [K, P].—SrerHen H. Cranpatu, Editor, Notes for the M. I. T. Special Summer 
Program on Random Vibration, The Technology Press of the Massachusetts In- 
stitute of Technology, Cambridge, Massachusetts, 1958, viii + 418 p., 196 Fig., 
22 cm. Price $10.00. 


The first four chapters of these notes cover mechanical vibrations with deter- 
ministic excitation (S. H. Crandall), the mathematical description of random proc- 
esses (W. M. Siebert), stochastic processes of mechanical origin (H. Poritsky), and 
statistical properties of response to random vibrations. The remaining eight chap- 
ters cover instrumentation for random vibration analysis (T. F. Rona), simulation 
equipment (K. J. Metzgar & D. E. Priest), and a very wide variety of discussions 
of special fields in which random vibration analysis is important: for example, struc- 
tural damping (T. H. H. Pian), the fatigue of metals (F. A. McClintock), the re- 
sponse of structures to random pressures and to jet noise in particular (A. Powell), 
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the estimation of sound-induced missile vibrations (I. Dyer), and mechanical design 
for random loading (R. M. Mains). 

The treatment of the basic material in the early chapters is good, but rather ab- 
breviated, especially if the students taking the course are not already reasonably 
familiar with the leading concepts and the mathematical methods that have been 
developed for the study of random vibrations. A person without such a background 
would probably not get much out of this brief and rapid treatment, and one well 
equipped would presumably find it superfluous. The chapters on applications cover 
a wide field, and in each case would appear to be of interest only to specialists. It 
seems to me that a far sounder course would have resulted if no attempt had been 
made to cover detailed applications, and if instead the course had been devoted to 
a more careful and detailed presentation of the basic theoretical material, followed 
by the chapters on instrumentation and on testing devices. 

C. V. L. 8. 


59 [L].—Lovis Rosin, Fonctions Spheriques de Legendre et Fonctions Spheroidales, 
Gauthier-Villars, Paris. Tome 1, 1957, xxxvi + 201 p.; Tome 2, 1958, vii + 
384 p., 24 em. Price $18.50. 


The mathematical community will be grateful to Dr. Robin for providing an up- 
to-date and virtually encyclopedic work on Legendre functions, which is certain to 
take its place among the most important books on its subject, and in many respects 
is likely to become the standard work of reference in this field. 

The book is written primarily from the point of view of applied mathematicians, 
having in mind the numerous classical applications of Legendre functions and 
spherical harmonics as well as some modern applications, such as in antenna theory. 
It gives a wealth of formulae and results. These are clearly stated, with the condi- 
tions of validity scrupulously described in each case. The majority of the results are 
proved in detail; in cases where proofs have been omitted or severely curtailed, ref- 
erences to the literature are included. Errors found in the literature have been cor- 
rected, so that some of the results appear here correctly for the first time, and others 
are new or due to the author, who has been engaged in research on Legendre func- 
tions for a number of years. 

The contents of the first two volumes of this treatise will now be summarized 
and briefly commented upon. 

Topics considered in Chapter I include the separation of the wave equation in 
curvilinear coordinates, Legendre’s differential equation, Legendre polynomials and 
their properties, and Legendre functions of the second kind and their properties. 

Chapter II deals with the associated Legendre functions P,,.” and Q,,” for integral 
values of m and non-negative integral values of n. A special feature of this chapter 
is the careful consideration of P,” for negative integers m, including m < —n. 

Chapter III treats spherical harmonics of (positive, negative, or zero) integral 
degree. Included therein is an account of Maxwell’s generation of spherical har- 
monics and transformation formulae to new coordinate axes, as well as the addition 
theorem for Legendre polynomials and expansions in spherical surface harmonics. 
Appendices contain explicit expansions and formulae for Legendre functions and 

associated Legendre functions with integral m and n. 
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The second volume begins with Chapter IV, which occupies 211 pages. It is de- 
voted to an investigation of the associated Legendre function P,,” (z) for general 
(complex) values of m, n, and z. Corresponding to such values, this function is de- 
fined by a double loop integral. The connection with the hypergeometric equation 
is established, and numerous hypergeometric series are given for P,,” and Q,. The 
analytic nature of the Legendre functions as functions of the three complex variables 
m, n, z is carefully described, and expansions are given also in the exceptional (loga- 
rithmic) cases as well as for the derivatives of the functions with respect to m and n. 
Chapter V is concerned with approximations, asymptotic expansions, and inequali- 
ties. This is the largest published collection of material of this kind, and some of it 
is new or has not previously appeared in a book. Chapter VI relates to expansions in 
Legendre polynomials, associated Legendre functions, and spherical surface har- 
monics. Both the (pointwise) convergence and various kinds of summability of these 
expansions are investigated, and some examples of such expansions are given. Ex- 
pansions involving P,,,” (cos @) are also discussed when the n; are roots of a trans- 
cendental equation. The Gibbs phenomenon is described. 

A forthcoming third volume will contain four chapters on addition theorems, 
zeros, applications, and related functions, including spheroidal functions, and an 
appendix on numerical tables of Legendre functions. 

The book is very lucidly written and can be recommended both for study and as 
a work of reference. References to the literature are given throughout the text at ap- 
propriate places, but there is no systematic bibliography in the two volumes under 
review. 

A. ERDELYI 
California Institute of Technology, 
Pasadena, California 


60 [P, S, T, X, Z].—L. H. Hersacu, Editor, Symposium on Digital Computing in 
the Chemical and Petrochemical Industries, New York University-International 
Business Machines Corporation Symposium, New York, 1958, vi + 189 p., 28 
em. Price $7.00. 


This report contains a series of papers given at the Symposium on Digital Com- 
puting in the Chemical and Petrochemical Industries held at New York University 
in May 1958. This was the second symposium in this series. The first symposium 
was held in early 1957 and was devoted to the application of computers in the air- 
craft industry. A list of papers follows. 

1. “A Numerical Solution to the Miscible Displacement Equation” by D. U. 
von Rosenberg. 
2. “Thermodynamic Properties of Neon” by J. B. Butt, J. C. Leutwyler, 

B. F. Dodge & R. W. Southworth. 

. “Numerical Experiments in Chemical Reactor Design” by L. Nemerever. 
. “Design of a Distributed Control System” by O. Bilous. 
. “Supervision of an Operating Plant by a Digital Computer” by W. D. 
Mohr. 
6. “Multicomponent Ion Exchange Column Calculations” by J. Dranoff 
& L. Lapidus. 
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7. ‘Multicomponent Distillation Calculation on the IBM 704” by J. Green- 
stadt, Y. Bard, & B. Morse. 
8. “Source-Sink Method for Oil-well Fracture Problem” by M. Friedman, 
E. Mehr, & C. Harris. 
9. ‘Solution of Diffusion Equations for Photochemical Initiation by Light that 
is Intermittent in Space” by R. M. Noyes. 
10. ‘““Molecular Orbital Calculations” by J. Auffray. 
11. “The Role of the Digital Computer in the Structure Determination of 
Complex Organic Molecules by X-ray Analysis” by G. A. Jeffry. 
12. ‘‘The Use of Digital Computers in X-ray Crystallography” by W. G. Sly. 
13. “The Engineer, the Computer, and the General Purpose Program” by 
E. P. Bartkus & W. M. Carlson. 
14. “The Use of Linear Programming Techniques in Crude Oil Production” 
by A. 8. Lee & J. Aronofsky. 
15. ‘“‘Product Allocation by Linear Programming” by K. H. Shaffir & T. Zang. 
6. “Library Programs for Chemical Engineering Applications” by G. Trimble. 


H. P. 


61 [P, X).—Rvex V. Cuurcui.1, Operational Mathematics, McGraw-Hill Book Co., 
Ine., New York, 1958, ix + 337 p., 23 em. Price $7.00. 


The first edition of this book, published in 1944, was entitled Modern Opera- 
tional Mathematics in Engineering. With the exception of the last chapter, which is 
concerned with Fourier transforms, both editions are devoted to an exposition of the 
theory of the Laplace transformation and its applications to physical problems, at 
about the junior-senior level of difficulty. The present edition is an enlargement by 
about ten per cent of the old, with many sections extensively rewritten and new 
problems added. Considerably rewritten are the portions of the book concerned 
with fundamental properties of analytic functions and with Sturm-Liouville systems. 
The arrangement of chapters is changed not at all, and of subject matter within 
chapters, very little from the earlier edition. The exposition is sound and lucid, and 
the book should have wide appeal. 

Tables are essentially as in the earlier edition, namely, a table of Laplace trans- 
forms containing 125 entries, and a short table of finite Fourier sine and cosine trans- 
forms. 

Joun W. GREEN 
Department of Mathematics, 
University of California, 
Los Angeles 24, California 


62 [S].—J. P. Scuirrer, Tables of Charged Particles Penetrabilities, ANL-5739, Ar- 
gonne National Laboratory, Lemont, Ill., 1957, 59 p., 28 em. Price $.40. Office 
of Technical Services, Washington 25, D. C. 
Charged particles penetrabilities defined as A,? (x, g) = Fx? (x, g) + G.? (2, g) 
where F, and G, are the regular and irregular Coulomb functions, have been cal- 
culated for L = 0 to 4, « = p/2 = 0.05(.01)0.70, g = ~/2pn = 1.6(.1)7.0. The 
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range covered corresponds roughly for protons on nuclei to z ~ 15 (—1)75 and 
bombarding energies up to 0.7 of the Coulomb barrier height. 

The calculations were carried out on an IBM 650, neglecting F',? and using the 
approximate formulae of Feshbach, Shapiro and Weisskopf [1] for Gp and G)’. 

The tables are reproduced from the machine listings in floating decimal format 
with all quantities given to four figures. The accuracy was checked where possible 
against the entries of Bloch, et al [2]; regions of agreement to better than two per 
cent and better than five per cent with the latter reference are indicated. Other re- 
gions are admittedly meaningless. 

M. MELKANOFF 
Numerical Analysis, 
University of California, 
Los Angeles, California 

1. H. Fesupacu, I. SHaprro & V. F. Weissxorr, Tables of Penetrabilities for Charged 
Particle Reactions, New York Operations Office of AEC 3077, 1953. 


2. I. BLocn, ET au., ““Coulomb functions for reactions of protons and alpha-particles with 
the lighter nuclei,’”’ Revs. Mod. Phys., v. 23, 1951, p. 147-182. 


63 [S, X].—B. Nose, Methods Based on the Wiener-Hopf Technique for the Solution 
of Partial Differential Equations, Pergamon Press, Inc., New York, 1959, 246 p., 
21 em. Price $10.00. 


This excellent book is a definite contribution to modern applied mathematics. 
It brings together, in a well-integrated manner, results which are widely dispersed 
in the scientific literature. 

Professor Noble has isolated the fundamental idea behind the Wiener-Hopf 
technique. He recognizes that it deals with certain factorization problems for func- 
tions of a complex variable and that it is only indirectly related to the Wiener-Hopf 
integral equation with which it is usually associated. : 

The classical Wiener-Hopf procedure starts from the Wiener-Hopf integral 
equation: 


(1) AA(z) + f(x) = I k(x — t)h(t) dt, Pe eaten, 


In the applications which are considered in Noble’s book, \ vanishes. Define h(x) 
and f(x) to be 0 for negative x, and let e(x) be the value of the integral in (1) for 
negative x. Introduction of Fourier transforms reduces equation (1) to: 


(2) F(a) + E_(a) = K(a) H,(a), 


where F, and H, are analytic in an upper half plane, E_ is analytic in a lower half 
plane, K is analytic in a strip. K(a) can be factored as K,/K_ and F,K_ can be 
split as C, + C_. This yields 


(3) E_(a)K_(a) + C_(a) = H,(a)K,(a) — C,(a), 


where E_ and H, are unknown. The remarkable fact is that both can be determined 
from the single equation (3). The left and right sides of (3) usually have a common 
strip of analyticity, so that either side is the analytic continuation of the other. Both 
sides are therefore equal to the same entire function M(a), which often turns out to 
be a constant or a polynomial. M(a) can be determined from the asymptotic be- 
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havior of the functions appearing in equation (3). Some of these asymptotic beha- 
viors are postulated on physical grounds, the others being found from equation (1). 

Once M(a) has been determined, H, and E_ are computed from (3), and inver- 
sion gives h(x). 

A typical physical application occurs in the diffraction of waves by a half-plane. 
The boundary value problem can be reduced to an integral equation of type (1), 
where k(x — t) is the Green’s function for the partial differential equation, h(x) is 
the current on the screen, and f(x) is the incoming wave on the screen. 

Noble points out that the functional equation (2) can be obtained directly from 
the boundary-value problem without reference to the integral equation. The advan- 
tages of this point of view are the following: 

a. K(a) can be obtained with ease by transforming the partial differential equa- 
tion directly without introducing the Green’s function k explicitly. 

b. We are led to possible further applications of the Wiener-Hopf technique to 
problems which cannot be formulated as an equation of type (1). The generalized 
two-part boundary-value problems of 8. Karp fall in this category. The Wiener- 
Hopf technique can be extended to cope with these problems, the appropriate trans- 
form being a Mellin, Hankel, or Kantorovich-Lebedev transform rather than the 
ordinary Fourier transform used for problem (1). 

There are many other interesting features of Noble’s book. A comparison of the 
Wiener-Hopf technique with Titchmarsch’s dual integral equation approach is 
given. One also finds a careful discussion of the relationship between the Wiener- 
Hopf problem and the Hilbert problem. An extensive collection of problems (taken 
from acoustics, electromagnetic theory, fluid flow, elasticity, and neutron diffusion) 
is worked out by the author. He has also included a wealth of non-trivial exercises, 
many of which are likely starting points for further research. 

There are a number of topics omitted by Professor Noble, which I would have 
liked to see included. These omissions do not detract appreciably from the over-all 
excellence of the volume, but they might easily be remedied in the next edition: 

a. The work on approximate methods in connection with Wiener-Hopf prob- 
lems. The powerful techniques developed by Harold Levine receive only cursory 
mention. 

b. The attempt to develop a two-dimensional Wiener-Hopf technique. J. Rad- 
low has made a start toward the solution of the quarter-plane problem. This leads 
to difficult questions in the theory of functions of two complex variables. 

c. There is a fairly extensive literature in probability theory on the Wiener-Hopf 
method. The interest is with solutions of integral equations of type (1) which are 
required to be cumulative probability distributions. Significant results have been 
obtained in this direction by F. Pollaczek, W. L. Smith, and F. Spitzer. 

In conclusion, Professor Noble is to be thanked for his valuable analysis of an 
important tool of mathematical physics. 

Ivar STAKGOLD 


Office of Naval Research, 
Washington 25, D. C. 


64 [W, X].—K. J. Arrow, L. Hurwicz, & H. Uzawa, Studies in Linear and Non- 
Linear Programming, Stanford Mathematical Studies in the Social Sciences, IT, 
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Stanford University Press, Stanford, California, 1958, iv + 229 p., 25 cm. Price 
$7.50. 


The general mathematical programming problem can be defined as: Determine 
the n-vector x which maximizes f(x) subject to the conditions 


gi (x) = 0, +=1,2,...m 
z.2-0. 


Linear programming is the special case where f(x) and the constraints g,(x) are 
linear. 

Since the development of the Simplex Method, linear programming has had many 
applications to management problems. Also, as a representation of the economic 
problem of optimum allocation of scarce resources, linear programming (with its 
duality theorem) has become an important concept in economic analysis. 

Refinements to linear programming have been extensive in recent years, but de- 
velopments in the more general aspects of mathematical programming have not 
been nearly as rapid. Computational procedures have been developed for the case 
where the functional f(x) is quadratic and where f(x) is convex and separable for 
each x;. All of these procedures have used the Simplex Method, or some variant, as 
the basic computational technique. Of more general nature is the basic paper by 
H. Kuhn and A. W. Tucker entitled ‘““Non-linear Programming,”’ which appeared 
in Proceedings of the Second Berkeley Symposium on Mathematical Statistics, 1951, 
p. 481-492. In this paper, Kuhn and Tucker give conditions under which the mathe- 
matical programming problem is equivalent to finding the saddle point of the La- 
grangian function, g(x, wu) = f(x) + u-g(z). 

The Kuhn-Tucker theorem on non-linear programming provides the basis of 
this book, which is in three parts. Part I, entitled Existepce Theorems, is mainly 
concerned with extensions of the Kuhn-Tucker theorem. Part II, entitled Gradient 
Method, is devoted to extensive studies of the gradient method applied to the 
saddle point of the Lagrangian function. Of special interest is the development of a 
modified gradient method which converges for all concave programming problems, 
(that is, where f(x) and g(x) are concave functions of x) without restriction. Exam- 
ples of this modified gradient method applied to two small linear programming 
problems are presented in detail. Part III, entitled Methods of Linear and Quad- 
ratic Programming, discusses the solution of three specific problems whose special 
structures are used to develop methods of solution. Also, a method for solving linear 
programs based on a theorem on extreme points of convex polyhedral cones is pre- 
sented. 

This book, which is a collection of related papers published for the first time, is 
the second in the Stanford Mathematical Studies in the Social Sciences series. Like 
the first, Studies in the Mathematical Theory of Inventory and Production, by Arrow, 
Karlin, and Scarf, it is a valuable addition to the literature. It should be of consider- 
able interest to mathematical economists and management scientists concerned 
with generalizations of linear programming. 

GEorRGE SuzvuKI 
Department of the Navy, 


Navy Management Office, 
Washington 25, D. C. 
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65 [W, X].—A. W. Tucker & R. D. Luce, Contributions to the Theory of Games, 
Annals of Mathematics Studies, v. 4, Princeton University Press, New Jersey, 
1959, 453 p., 24 em. Price $6.00. 


This is the fourth in a series of books presenting papers relating to the theory 
of games; the editors claim that it is the last. The whole series has been a most 
valuable collection of papers; if nothing else, they indicate the mathematical, and 
hence the numerical, difficulty which is to be encountered in attempts to treat 
competitive situations quantitatively. 

The present volume is perhaps further removed from computational method 
than any of the earlier ones. It still contributes to the description of the quantitative 
difficulties occurring in the description of competitive situations. 

One of the most valuable contributions of this volume is a bibliography by 
Dorothea M. Thompson and Gerald L. Thompson; this is a list of about a thousand 
of the most readily available contributions to game theory. 

A second valuable contribution is the introduction, presumably written by 
the editors, outlining what is in the volume and what the general significance of 
each paper is. This introduction constitutes what would be the most valuable 
review possible for this volume, but it is unfortunately too long to reproduce here. 

A third interesting item is a translation into English by Mrs. Sonya Bargmann 
of von Neumann’s famous early paper “Zur Theorie der Gesellschaftsspiele,” 
which appeared in 1928. 

While the present volume has only these fairly remote offerings to make to 
computation, the fact that so much present day computing is concerned with 
economic problems makes the material worth noting here. In addition, the reviewer 
would like to add his applause to that already bestowed on the editors of this series 
of volumes for their contributions in publishing research works in a single field 
grouped into a single volume with a good editorial introduction to guide the reader. 
The time the editors have expended on this task is certainly great, but to the reader 
the contribution is most valuable. An expansion of this idea would be most valuable 
if the editorial time could be made available. Specifically, it would seem appropriate 
to have journal editors relay some of the information they receive in the form of 
referees’ reports to aid the cursory readers of their journals. 

C. B. T. 


66 [X].—GLEeNN James & Ropert C. James, Editors, Mathematics Dictionary, 
2nd Edition, D. Van Nostrand, Princeton, 1959, 546 p., 24 em. Price $15.00. 


The new edition of this well known reference is about 25 percent longer than 
the 1949 edition. The new material includes terms from “modern algebra, number 
theory, topology, vector spaces, the theory of games and linear and dynamic pro- 
gramming, numerical analysis, and computing machines.” Entirely new are the 
four foreign language dictionaries of mathematical terms in the appendix: French— 
19 pages, German—16 pages, Russian—22 pages, and Spanish—15 pages. 

The dictionary contains a wealth of information and makes interesting browsing. 
While the book is a ‘‘must’’ for mathematical libraries, it is not difficult to find 
many minor defects. Since the editors, in their preface, invite comments, here are 
some of mine. 
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Sins of omission and commission: Where are “‘Class Number,” ‘Continuum 
Hypothesis”, “Diophantine Approximation’’, ‘‘Dirichlet Series’, and “Distribu- 
tion, (L. Schwartz)”? On the «ther hand must we have “Commercial Bank”’, 
“Horsepower”, ‘‘Mariner’s Compass’’, etc.? Sins of misplaced emphasis: After 
““MObius function” one would expect something about the Mébius inversion for- 
mula, but instead we find the more sophisticated, but less fundamental, relation 
to the Riemann Hypothesis. Improper definitions: ‘‘Unique Factorization Theorem” 
and ‘Factorial Series’ are not properly defined. Improper spelling: ‘“Fresnal”’, 
“Marriot”. While ‘“Tchebycheff” and even ‘“Tschebyscheff” are not “wrong,” the 
more modern “‘Chebyshev”’ could at least be given as a cross-reference. Typog- 
raphy: “Joukowski Transformation” is listed as a subheading of ‘Jordan.’’ Non 
sequitur: ‘Fibonacci Sequence” has a mysterious reference to ‘‘Farey Sequence”’, 
but this should cause the reader little inconvenience since it is only one page away. 

The appendices contain elementary tables of logarithms, trigonometric func- 
tions, compound interest, etc. There is also a table of ‘“Denominate Numbers” 
(2 pints = 1 quart, 7 days = 1 week, 20 grains = 1 scruple, 2352 pounds = 1 
Cornish mining ton, 24 sheets = 1 quire, 500 sheets = 1 long ream, etc.). British 
readers susceptible to acute nostalgia should be forewarned of the item which 
states that the British monetary pound is normally equal to 5 dollars. The editors 
seem especially fond of the fact that a gallon equals 231 cubic inches, since this is 
repeated three times. With due respect to them, could not the space be better 
used in the tables of integrals which follows? The type size there is almost un- 
believably small. 

The foreign-language dictionaries and the hundreds of new terms such as ‘“‘Haar 
measure”, ‘“Tychonoff space’’, etc., add considerably to the value of this useful 
reference book. 

‘ D.S. 


67 [Z|.A. D. Boor & K. H. V. Bootu, Automatic Digital Calculators, Second 
Edition, Academic Press Inc., New York, 1956, 261 p., 21 cm. Price $6.00. 


This book is a non-mathematical introduction to the principles and applications 
of computers employing tubes and other electronic devices. It is a very general 
coverage of the entire digital computing field. 

The first three chapters of the book are devoted to the evolution of calculating 
machines to the modern general-purpose electronic digital computer. 

The next four chapters discuss the overall design of a computing system: the 
memory, control, arithmetic unit, and input-output units. Each unit is described 
in terms of serial, parallel, and decimal operation, along with the peculiarities of 
each type of operation and number system being used. 

Four chapters are devoted to a general discussion of circuits and hardware which 
make up the logic of the different systems. The different types of storage, such as 
magnetic core, magnetic drum, magnetic disc, magnetic wire, mercury delay lines, 
magnetostrictive delay lines, and other types of storage, are described in terms of 
general operating principles. 

The remaining five chapters are devoted to coding for a digital computer and 
digital computer applications. The definition of a code and discussion of its form 
and contents are described. 
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The book gives the reader a very good introduction to the general principles of 
coding, operation and application of high-speed electronic computers and computer 
components. Some previous background in mathematics, or a previous knowledge 
of computer circuitry, is required for a thorough understanding of the contents. 
In clear and concise terms, the authors cover the subjects in seventeen well written 
chapters. The book is written primarily for new workers in the field. Engineers in 
other fields who wish to become acquainted with the subject of computers will find 
the text both interesting and informative. The book lays a firm foundation for 
subsequent reading in the digital computer field. It is not intended, however, that 
the book be of use to engineers who are well informed in computing techniques and 
equipment. The excellent illustrations of circuitry serve well as an aid to com- 
prehension. 

CueEstTerR H. WALLIN 
Computing Laboratory, BRL, 
Aberdeen Proving Ground, Maryland 


68 (Z).—Warrts 8. Humpnurey, Jr., Switching Circuits with Computer Applications, 
McGraw-Hill Book Company, Inc., New York, 1958, viii + 264 p., 23 cm. 
Price $8.50. 

After an introductory chapter mostly devoted to the binary number system 
then follows a fairly long chapter on Boolean algebra: the presentation is intuitive 
and practical, the basic ideas and results, for example, being set forth by means of 
Venn diagrams. This is immediately applied in Chapter 3 to the study of relay 
contact networks. Chapter 4 is devoted to codes; included are discussions of the 
various kinds of undetected errors using binary-coded decimal and excess-three 
codes, of modulo-nine checking, of parity checking, and of error-correcting codes. 
Chapter 5 is on minimizing aids, and is done almost entirely in terms of 
the Karnaugh map. Chapter 6, on circuit logic, covers two- and three-level diode 
logic and transistor logic. Chapter 7, on Boolean matrices, deals with connection 
matrices of networks, their modification and manipulation. This lays the founda- 
tion for the treatment, in Chapter 8, of bilateral networks, both symmetrical and 
antisymmetrical, and their reduction. Chapter 9, on cascaded networks, deals 
with complex networks built up by cascading simpler ones, and considers clocked 
and unclocked pulse logic and level logic. Finally, Chapter 10 is devoted to se- 
quential circuits, defined as circuits containing memory (or time-delay) elements. 

The discussion throughout is practical and didactic rather than theoretical. 
Many worked examples are included in the text; these are very instructive, in most 
cases are easy to follow, and supplement and elucidate the occasionally rather 
sketchy presentation of ideas. This should make the book particularly valuable for 
self-study, though it is also very suitable as a textbook. The emphasis throughout 
is on practical methods of minimizing switching networks, and no single method or 
procedure is presented as a cure-all. The very sensible conclusion of the introduction 
to the chapter on Boolean matrices could well be taken as a maxim: “It is, how- 
ever, no different from the other switching-circuit design techniques: it is an aid to 
the designer and does not replace skill and experience.” 


i Vi. Boo 











NOTES 


Announcement of New Title for 
Mathematical Tables and Other Aids to Computation 


In order to reflect the broadened scope of Mathematical Tables and other Aids to Computa- 
tion, and to satisfy an increasing need in this country for a journal devoted to numerical 
analysis and computation, the members of the editorial committee of MTAC voted unani- 
mously at their last meeting to recommend a change in the name of the journal. Mathematics 
of Computation has been selected as the new title. This change has been approved by the 
executive committee of the Division of Mathematics, National Academy of Sciences-Na- 
tional Research Council, and will become effective with the January 1960 issue. 

The new title for Mathematical Tables and other Aids to Computation in no way represents 
a diminished interest by the editorial committee in mathematical tables, which will continue 
to be a subject of prime emphasis, as in the past. It recognizes, however, an increased interest 
in other areas in the field of Mathematics of Computation, which have grown in importance and 
in which rapid advances are being made in the present era of technological progress. The sub- 
scribers and readers of Mathematical Tables and other Aids to Computation may expect future 
issues of the journal to be a continuation of MTAC, with similar style, format and character 
of contents, but with increased emphasis on modern advances in the theory and application 
of computational methods. It will appear under the abridged title, Mathematics of Computa- 
tion. 7 - 


New Journal 


Information and Control is a journal devoted to the publication of original papers in the 
theories of communication, computers, and automatic control. Papers on experimental results 
in this field are also accepted. It is published by Academic Press Incorporated, New York, 
New York. Editors: Leon Brillouin, Colin Cherry, and Peter Elias. Volume I covers the years 
1957-58. 

; E. I. 


Conference on Communication and Scientific Information 


The October 1958 issue, Volume 2, Number 4 of the IBM Journal of Research and De- 
velopment is devoted to ten papers (plus a related contributed article) presented at the Con- 
ference on Communication and Scientific Information held on the occasion of the dedication 
of the IBM Research Laboratory, San Jose, California, May 26, 27, 1958. New developments 
are presented in the use of machines for efficient information processing and problem solving. 

E. I. 


Funds for the Publication of Mathematical Tables 


Funds are available to subsidize the publication of suitable mathematical tables through 
the Division of Mathematics of the National Research Council of the National Academy of 
Sciences. Inquires should be sent to: 

Professor George E. Forsythe 

Chairman, Committee on Revolving Fund for Publication 
of Mathematical Tables 

Mathematics Department, 

Stanford University, 

Stanford, California 
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CORRIGENDA 


A. M. Ostrowsk1, ‘‘On Gauss’ Speeding Up Device in the Theory of Single Step Itera- 
tion,’”’ MTAC, v. 12, 1958, p. 116-132. 
p. 124, formula (36), add to the last right-side term the factor £. 
). 124, immediately preceding formula (36), add, ‘‘for a sufficiently small positive £.”’ 
p. 124, end of section 21, add, “‘replacing a; + a2 by a; + a2 + E(n — 2).” 
WaLTEeR GAUTSCHI 


— 


Notes, MTAC, v. 12, 1958, p. 319, item (4) announces the publication of the monograph 
series Acta Polytechnica, Applied Mathematics and Computing Machinery Series. The more 
accurate title of this series is Acta Polytechnica Scandinavica, Mathematics and Computing 
Machinery Series. This joint Scandinavian series is published under the auspices of the Scan- 
dinavian Council for Applied Research. 

Fe 
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